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SUMMARY 


Development of the conservation equations in the first part starts from 
the Boltzmann equation, A general transport equation for an arbitrary particle 
property is derived; from this the species mass, momentum, and energy conser- 
vation equations are obtained. Collision terms are treated in detail in the 
approximation of two independent, displaced Maxwellian distribution functions, 
which takes into account the effects of different flow velocities and tempera- 
tures. A new formalism is developed for these terms, simplifying analytic 
evaluation; this is illustrated by analytic evaluation for a sufficient variety 
of interaction potentials to meet most ionospheric requirements. The three- 
fluid approximation is applied to the species conservation equations, resulting 
in separate sets of equations for electrons, ions, and neutral particles. 
Order-of-magnitude estimates, based on extreme values for ionospheric proper- 
ties in the altitude region 90 km to 800 km, are used to delete terms which 
are not important for the specified ionospheric conditions. Closure of the 
conservation equations is accomplished through transport tensors and coeffi- 
cients. A condensation of Shkarofsky's method for calculating electron trans- 
port tensors is presented. Use of this technique is facilitated by the deri- 
vation of a new method for representing the electron-neutral collision frequency 
velocity dependence as a power law. Transport coefficients for ions and 
neutral particles are also presented. The resulting set of conservation 
equations should provide appropriate starting points for a wide variety of 
studies, both theoretical calculations and data analysis. 


In the second part, these conservation equations are applied to the 
detailed examination of auroral E region neutral winds observed by the inco- 
herent scatter radar at Chatanika, Alaska, during geomagnetic disturbances (15 
May 1974). Two primary objectives of this study are: (1) to determine the 

effects of altitude structure in ion and neutral velocity fields on neutral 
velocities derived from the radar observations; and (2) to determine the 
relative importance of ion drag and auroral heating in generating these winds. 

Ion and neutral momentum equations for the ion drag model are solved numeri- 
cally, using observed electric fields and electron number densities. The 
coriolis force is included; its importance is demonstrated in model calculations. 

The following results are obtained from these velocity calculations. 

Large vertical gradients are found in the calculated velocities for altitudes 
below about 130 km. As a consequence of this structure and fluctuations in 
the electron density profiles, the data analysis procedure of Brekke, et al. 
(1973) for obtaining neutral winds from radar data are found to underestimate 
the wind speed by up to 40 percent, but it represents the direction and tem- 
poral structure reasonably well. Comparison of observed neutral velocities 
with calculated values shows that ion drag alone cannot account for the obser- 
vations. An equation is derived to estimate the pressure gradients required 
to resolve the discrepancy between calculated and observed neutral winds. 
Accelerations due to these pressure gradients are of the same order as those 
due to ion drag, but at least an order of magnitude larger than those due to 
solar heating. Directions of the pressure gradients are consistent with 
expected locations of auroral heating. During geomagnetic disturbances, ion 
drag and auroral heating both appear to play important roles in the generation 
and modification of neutral winds. 
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Section I 


INTRODUCTION 


Recent interest in upper atmospheric dynamics has been directed toward 
high latitudes, both because of increased observations in this region and 
because of increased awareness and understanding of the importance of the 
coupling between the magnetosphere and ionosphere (Folkestad, 1972; Chappel, 
1974) , which is most prominently displayed in the vicinity of the auroral oval 
(Akasofu, 1968). Observations of ionospheric properties in auroral latitudes 
(-60° to 80° geomagnetic latitude), made by a variety of remote and in situ 
techniques, have found electric fields, ionization number densities, flow 
velocities , and temperatures to be highly variable in space and time relative 
to conditions at lower latitudes. Time scales for significant variations of 
these properties are typically on the order of minutes at high latitudes 
(e.g.. Banks, et al. , 1974), compared with hours at middle and low latitudes 
(cf. Rishbeth and Garriott, 1969). Similarly, spatial scales for significant 
variations , on the order of several thousand kilometers at low and middle 
latitudes , are an order of magnitude smaller at high latitudes and even smaller 
near auroral arcs (e.g., Chan and Colin, 1969; Jelly and Petrie, 1969). These 
relatively small spatial and temporal scales complicate the comparison and 
interpretation of observational data (Maynard, 1972). They also increase the 
difficulties of theoretical studies because appropriate boundary conditions 
are uncertain, and approximations which permit physical insight and mathemati- 
cal tractability in lower latitude problems are not valid for high latitude 
conditions. In addition, the magnitudes of electric fields and associated 
plasma drift velocities are an order of magnitude larger at high latitudes 
(Maynard, 1972; Banks, et al. , 1974). As a consequence, different and some- 
times new physical considerations enter into the formulation of the governing 
equations. The primary objectives of this study are to develop a theoretical 
framework adequate for treating problems of the dynamics of magnetosphere- 
ionosphere interactions at high latitudes, and to use this framework in the 
investigation of such a problem. 
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Two approaches, continuum and particle, may be taken in formulating the 
governing equations for upper atmospheric gases. The goal of both approaches 
is to determine the relations among and the behavior of the macroscopic 
observables, those measured by experiments which detect the average properties 
of a small but finite volume of the gas, specifically: number density n, flow 

velocity v, and translational temperature T. In the continuum approach, 
equations governing the macroscopic observables are formulated directly from 
physical considerations of the bulk behavior of the* gases. This involves the 
introduction of transport coefficients, which indicate how readily spatial 
inhomogeneities in the density, momentum, and energy of the gas approach a 
homogeneous state. These coefficients must be determined experimentally in 
the continuum approach. Proceeding from a particle viewpoint involves formu- 
lating equations which incorporate interactions of individual particles with 
fields and other particles, using statistical methods to determine the macro- 
scopic behavior of the aggregate. To the order of approximation appropriate 
to the upper atmosphere, both methods lead to essentially the same macroscopic 
equations. Important advantages of the particle formulation are (1) the 
enhanced insight gained into the physical processes which affect the macro- 
scopic observables, (2) the possibility of direct calculation of the transport 
coefficients, and (3) the assurance of a consistent treatment, especially when 
new physical processes must be considered. In addition, the particle approach 
is particularly flexible in treating multicomponent gases, which recommends 
its use in upper atmospheric problems. This, then, is the approach to be used 
in this study. 

Relating the macroscopic behavior of a nonuniform, multicomponent gas to 
the microscopic interactions of its constituent particles is the task of 
kinetic theory. Methods used for accomplishing this end are discussed in 
detail by Chapman and Cowling (1970). The basic procedure begins with a 
Boltzmann equation for each particle species (e.g. , 'S 9 0^, e, 0+, 0 + , etc.); 
this equation governs the velocity distribution of each species. A general 
transport equation is obtained by multiplying both sides of the Boltzmann 
equation by an arbitrary function of particle velocity, and integrating over 
velocity space. When those particular particle properties which are conserved 
in (elastic) collisions (mass, momentum, energy) are used in this transport 
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equation, the conservation equations for number density, momentum, and energy 
result. In general, a set of conservation equations for each particle species 
is required to determine properly the aggregate behavior of the system. 
However, if several species have in common those properties which are of 
primary interest (e.g. , temperature and flow velocity), the appropriate equa- 
tions can be easily combined, resulting in a considerable reduction in compu- 
tational effort when the system of equations must be solved. This is the 
basis of the three-fluid approximation. 

As applied to a partially ionized plasma, such as the ionosphere, the 
three-fluid approximation implies that electrons, ions, and neutral particles 
have distinct properties, but that these properties are undifferentiated by 
species within these gases. Applicability to problems of ionospheric dynamics 
can be better understood from a brief review of typical ionospheric conditions 
which affect the dynamics. Basically, four factors play primary roles in 
determining and differentiating the motions of upper atmospheric gases: 

(1) relative number densities of different gases, (2) types of interaction 
between particles, (3) relative masses of particles, and (4) forces which act 
differently on different gases. These are examined in turn below. 

Although particle populations vary considerably in both space and time, 
Figures 1-la and 1-lb show typical daytime, midlatitude altitude profiles of 
constituent number densities. Both neutral and ion compositions are seen to 
shift from diatomic to monatomic species and toward the lighter species with 
increasing altitude. This is due to diffusive separation in the Earth’s 
gravitational field and results in a decreasing mean mass per particle with 
increasing altitude. Ions are only singly ionized (negative ions are negligi- 
ble above 90 km (Banks and Kockarts, 1973)), and charge neutrality (electron 
number density is equal to total ion number density) is assumed to hold every- 
where on the macroscopic scale of interest here. 

Strictly in terms of numbers, the upper atmosphere is seen to be a weakly 
ionized plasma at all altitudes shown. However, the collisional coupling 
between the gases also depends on the interactions between the particles and 
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on their relative masses. These considerations are taken into account to some 
extent in the momentum transfer collision frequencies, indicate the 

degree of momentum coupling among the gases. Altitude profiles of collision 
frequencies, based on the number densities in Figure 1-1, are shown in Figure 
1-2. Most of the altitude variation in Figure 1-2 is due to variation in 
total number densities, and to a lesser extent to the changing compositions 
(hence changing interactions) and changing temperatures. Temperature profiles 
used in computing these collision frequencies are shown in- Figure 1-3; these 
profiles are representative of daytime midlatitude conditions. 

In addition to collisional coupling among the species, electromagnetic 
forces, which selectively affect only charged particles, have a pronounced 
effect in differentiating the motions of ionospheric gases. For purposes of 
this study, detailed consideration is limited to the effects of the static 
geomagnetic field and electrostatic (curl free) fields, in addition to the 
Coulomb interaction between charged particles. The intensity of the geomag- 
netic field varies geographically, ranging from about 0.2 gauss to about 0.6 
gauss at ionospheric altitudes (Matsushita and Campbell, 1967). Since this 
range of values is relatively small, the more important variation is in the 
geometry. A somewhat schematic view of the magnetosphere is presented in 
Figure 1-4, showing the dipolar configuration of the geomagnetic field at 
ionospheric altitudes. Not indicated in this diagram is the fact that the 
dipole axis is tilted about 11 degrees with respect to the Earth’s rotational 
axis • 


The significance of this geometry derives from the Lorentz force (qv x B/c) 
which causes charged particles moving in a magnetic field to spiral about the 
field lines. This inhibits motion transverse to the field, so that charged 
particles are, to some extent, tied to the field lines. The strength of this 
coupling is indicated by the particle gyrof requency , 

o> s = |q s |B/m s c (Gaussian Units) (1-1) 

where q g is the particle charge, B is the magnetic field strength, m g is the 
mass, and c the speed of light. For a magnetic field strength of 0.5 gauss 
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Figure 1-2. REPRESENTATIVE COLLISION FREQUENCIES (v) AND GYROFREQUENCIES 
(w) OF UPPER ATMOSPHERIC GASES FOR DAYTIME, MIDLATITUDE 
CONDITIONS 



TEMPERATURE <°K» 


Figure 1-3. THERMAL STRUCTURE OF THE UPPER ATMOSPHERE, REPRESENTATIVE OF 
DAYTIME, MIDLATITUDE CONDITIONS (JACCHIA, 1971; NASA, 1971) 
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Figure 1-4. SCHEMATIC STRUCTURE OF THE EARTH'S MAGNETOSPHERE (BANKS AND 
KOCKARTS, 1973, p. B-305) 
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(representative of the continental United States), electron and mass-averaged 
ion gyrofrequencies are shown along with the collision frequencies in Figure 
1-2. Since collision frequencies and gyrofrequencies play similar roles in 
the charged particle equations of motion, they may be compared directly to 
determine whether the magnetic field or the other gases will dominate the 
motion of a given gas. 


From Figure 1-2 it is evident that the electrons are firmly coupled to 
the geomagnetic field at all altitudes shown, since they gyrate about a field 
line many times between collisions with either ions or neutrals. Ions, however, 
undergo a transition from strong coupling to neutral particles at low altitudes 
to strong coupling to the geomagnetic field at high altitudes. The character- 
istic distance associated with coupling to a magnetic field is the gyroradius , 
given by 


r 

s 



0) 


( 1 - 2 ) 


s 

where Wj_ is the particle velocity component perpendicular to the magnetic 
field. If thermal velocities are used for order-of-magnitude estimates, 
representative gyroradii are about 3 cm for electrons and 4 ra for ions, corre- 
sponding to T e = 2400°K (w^ - 3 x 10^ cm s ^) and = 1200°K (w^ - 8 x 10^ cm 
s 1 ). As a result of this restriction in the transverse direction, charged 
particles move most freely parallel to the geomagnetic field lines, that is, 
in a north-south direction over the magnetic equator and vertically over the 
magnetic poles. 


At this point, several questions can be examined with respect to appli- 
cability of the three-fluid approximation in studies of ionosphere dynamics, 
particularly at high latitudes where magnetosphere-ionosphere interactions are 
so important. Is a three-fluid approximation necessary, or would a two-fluid 
or even a single fluid approximation be sufficient? Is a three-fluid approxi- 
mation adequate, or must the different species of the ion and neutral gases be 
treated separately? If the answers to these are affirmative, are present 
three-fluid formulations adequate, or are further developments required? The 
answers to these questions essentially determine the scope of the theoretical 
framework developed in this study. 
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The first question inquires whether the macroscopic properties n, v, and 
T of the electrons, ions, and neutrals are sufficiently different to warrant 
consideration of three separate gases. Experimental and theoretical work on 
ionospheric temperatures have been reviewed by Banks (1969) and Willmore 
(1970). The prevalence of nonthermal equilibrium among the gases of at least 
the ionospheric F region is well established observationally and reasonably 
well explained theoretically. The basic mechanism for this condition is 
preferential heating of thermal electrons (energies £ 0. 3e V) by photoelectrons 
in the process of being collisionally thermalized. Preferential heating of 
electrons results from their small masses and long-range Coulomb interaction. 

In turn, electrons transfer energy collisionally to ions and neutrals, but 
preferentially to ions, again because of the long-range Coulomb interaction. 

At higher altitudes (above ~300 km), ions lose close thermal contact with the 
neutrals and, through a balance of energy gained from electrons and loss to 
neutrals, attain a temperature intermediate to those of the electrons and 
neutrals (e.g., see Figure 1-3). Theoretical calculations (Rees and Walker, 

1968) have suggested that heating due to large auroral electric fields may 
cause ion temperatures to exceed electron temperatures, although no direct 
observation verifying this is presently known. Thermal nonequilibrium in the 
E region has been a subject of some controversy for several years because of 
apparently conflicting results from different observational techniques (D’Arcy 
and Sayers, 1974). However, this is immaterial to the present question since 
conditions in the F Region clearly require consideration of distinct electron, 
ion, and neutral temperatures. 

For flow velocities, fewer observational data are available for several 
reasons. Flow velocity measurements are more difficult, particularly those 
which measure velocities for different gases simultaneously. Also, conditions 
in which substantial differences in flow velocities of the different gases are 
likely to occur are less prevalent in space and time than for thermal nonequi- 
librium. Nevertheless, there is sufficient evidence to establish that large 
differences in flow velocity, on the order of hundreds of meters per second, 
can and do occur at high latitudes. In this region, large electric fields, 
sometimes in excess of 100 mV m \ can occur due to sources far out in the magnet- 
osphere (Maynard, 1972). These fields are mapped from the magnetosphere into 
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the ionosphere along geomagnetic field lines (see Figure 1-4) , which act 
approximately as equipotentials (Farley, 1959; Reid, 1965). This results from 
the high conductivity along field lines relative to that transverse to field 
lines. (A qualitative grasp of this idea may be gained from the comparison of 
collision frequencies and gyrof requencies in Figure 1-2, recognizing that for 
, collisions inhibit motion along the field line while gyrations 
inhibit motion transverse to the field line.) Implicit in this conception is 
that electrostatic fields are transverse to magnetic field lines. The degree 
to which parallel electric fields are actually negligible in the ionsophere is 
presently the subject of considerable controversy (Maynard, 1972; Zmuda, et 
al. , 1974). However, the existence of large transverse ionospheric electric 
fields at auroral latitudes is an observational fact (Maynard, 1972; Banks, et 
al. , 1974) ; and for present purposes that is sufficient. 

Large transverse electric fields cause large electron and ion flow veloc- 
ities of several hundred meters per second. In the absence of collisions, 
this velocity is given by 

V = E * B c (Hall Drift). (1-3) 

B 

For B = 0.5 gauss, an electric field of 30 mV m ^ results in a drift velocity 

of about 600 m s At high altitudes where ion-neutral collision frequencies 

are much smaller than ion gyrof requencies (Figure 1-2) , electrons and ions 

drift with the same velocity, given by equation (1-3) (which is independent of 

charge sign). Thus, when an electric field suddenly becomes large, the charged 

particles respond rapidly, resulting in a large velocity difference between 

the charged and the neutral particles at high altitudes. Under the same 

circumstances at low altitudes, where the ions are collisiorally coupled to 

the motion of the neutrals, the electrons will again respond rapidly to the 

electric field (since w >> v ) while ions are constrained to move with the 

s e 

neutrals. So at low altitudes the large velocity difference is between the 
electrons and the heavy particles. Numerical calculations, by Fedder and 
Banks (1972), among others, demonstrate these effects theoretically. Barium 
cloud release experiments, for example, those analyzed by Meriwether, et al. 
(1973) , provide direct observational evidence of ion-neutral velocity differ- 
ences of several hundred meters per second. Comparison of high altitude and 
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low altitude incoherent scatter radar observations allows similar inferences 
to be drawn with respect to electrons and ions, if transverse electric fields 
are assumed to be mapped along magnetic field lines unchanged from F region to 
E region altitudes (Brekke et al. , 1973, 1974a). From these considerations, 
it is concluded that different flow velocities for electrons, ions, and neu- 
trals must be taken into account in any theoretical framework used for treating 
problems of dynamics at high latitudes. 

For the final macroscopic observable, number densities, different reason- 
ing must be used. If composition of either the ion or neutral gas is required 
as part of the solution to a problem under study, then a continuity equation 
for every particle species is required. If composition can be assumed and 
only total numbers of particles of each gas (electrons, ions, and neutrals) 
are of interest, as in some dynamical problems, then only a two-fluid approxi- 
mation is required, since the charge neutrality assumption eliminates one 
charged particle continuity equation. For present purposes, consideration 
will be restricted to the latter case, that in which composition is not required. 

With the necessity of considering separate properties of at least three 
gases established, the next question must be examined: Are three gases enough? 

For number densities, this question has been somewhat sidestepped by limiting 
the scope of applications to problems in which more than three gases need not 
be treated. For temperatures and flow velocities the question reduces to 
asking whether or not all ion species share a common temperature and flow 
velocity, and likewise for neutrals. Experimental techniques are not suffi- 
ciently selective to be able to measure the flow velocity or temperature of a 
single species of ion or neutral gas. However, theoretical calculations 
indicate that temperature differences between 0 + and H + in excess of 100°K 
could be expected to occur (Banks, 1967). Since H + is a minor ion below about 
800 km (see Figure 1-lb) , this will cause no difficulty for studies below that 
altitude; and even above that altitude it may be of little consequence. Since 
0 + is predominant between 200 km and 800 km, and below 200 km the ion tempera- 
ture is equal to the neutral temperature, thermal nonequilibrium among ion 
species is negligible for altitudes below about 800 km (and perhaps above). 
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No theoretical work has suggested that neutral species in the ionosphere may 
have differing temperatures. Similarly, no information is available to indi- 
cate that flow velocities of either ion or neutral gases might be signifi- 
cantly differentiated by species. Therefore, it appears that the three-fluid 
approximation is adequate for a wide variety of problems in the altitude range 
90 km to 800 km for all geographic locations. 

Finally, the question of the present state of development of multifluid 
formalism must be considered. Since gases interact by collisions (leaving 
aside collective phenomena of plasmas for the present) , the important multi- 
fluid effects are contained in the collision terms. Different temperatures 
mean different thermal velocities; and since collision cross sections are 
functions of the relative speed of the colliding particles in general, average 
collision frequencies depend on the temperatures of both gases. Intuitively, 
it is apparent that if the relative flow velocities are comparable to or 
larger than the thermal velocities, the collision frequencies will also be 

2 

affected. In the ionosphere, ion thermal velocities are of the order of 10 

3 

to 10 meters per second. Ion drift velocities of this order occur frequently 
at high latitudes due to relatively large electrostatic fields, as discussed 
previously. This is why proper treatment of flow velocities in the collision 
terms is important in problems of high latitude ionospheric dynamics. 

In the context of these ideas, various multifluid formalisms can be 
considered briefly; a more detailed examination is deferred to Section III. 
Chapman and Cowling (1970) present an extensive treatment of the kinetic 
theory approach to gas dynamics in the single- fluid approximation. Their 
treatment of gas mixtures is based on deviations of the individual species 
properties from the (global) properties of the mixture as a whole. S. T. Wu 
(1969, 1970) treats gas mixtures in a similar manner, but includes detailed 
consideration of inelastic interactions between the gases and with the radia- 
tion field; however, he includes no external forces (e.g. , electrostatic and 
magnetostatic fields). Because of the relatively loose coupling of the three 
ionospheric gases , more physical insight is to be gained by treating these 
gases individually, rather than in terms of global mixture properties. More- 
over, since deviations of the properties of the electron and ion gases from 
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those of the neutral gas (which dominates the mixture, due to much larger 
number densities) may be quite large, the global approach is likely to be 
either inefficient or mathematically troublesome. 

More appropriate for the ionosphere is the work by Burgers (1969). He 
considers explicitly mixtures of gases with arbitrarily different temperatures 
and flow velocities. However, his treatment of collision terms for the case 
of arbitrary flow velocities is limited to a few examples (hard sphere. Coulomb, 
and Maxwellian-molecule interactions) which only partially satisfy the require- 
ments of ionosphere studies. Moreover, since Burgers 1 (1969) study treats 
higher order approximations extensively, his framework (and notation) is more 
general and complicated than is necessary for ionospheric applications. This 
makes it cumbersome to use and difficult to extend to other interactions of 
interest. This latter objection has been remedied somewhat by Banks and 
Kockarts (1973), who take appropriate results (conservation equations) from 
Burgers (1969) and present them in a more transparent notation. The difficulty 
here is that only the results are presented, so that the mathematical origin 
of parts of the equations remains somewhat obscure. In particular, the colli- 
sion terms are presented with insufficient formalism to permit extension of 
the results to other interactions of ionospheric interest. In addition to 
these, separate treatments of the collision terms alone have appeared in the 
literature, to be discussed in detail in Section III; but none have extended 

results beyond those presented by Burgers (1969). Thus, there appears to be a 

need for a more thorough development of the collision terms. for the multifluid 
case, with extension to additional interactions of ionospheric interest. 

Conclusions which may be drawn from the considerations above are as 
follows: the three-fluid approximation is needed for and applicable to iono- 

spheric problems, particularly at high latitudes; and present formulations 
need further work on the collision terms. The first objective of this study 
is then to develop the formal framework of the conservation equations in the 
three-fluid approximation, suitable for application to ionospheric problems at 
all latitudes and 90 km to 800 km in altitude (E region and lower and middle F 

region). The second objective is to use this framework in the investigation 
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of a problem of current interest in ionospheric dynamics at high latitudes. 

In the formal development of the conservation equations the treatment is 
necessarily uneven, since parts are covered in many texts, while others have 
been little discussed. To ensure completeness and consistency, steps of the 
more familiar parts are at least outlined; detailed mathematics is reserved 
for those areas involving original formulations or results. 

The approach taken is to start with the Boltzmann equation for each 
species and, following procedures noted previously, to develop the general 
transport equation. This is then used to obtain species conservation equations 
for number density, flow velocity, and energy. These steps are presented 
briefly in Section II, along with an introduction of the notation and physical 
ideas associated with the mathematics. 

Detailed treatment of the collision terms is reserved for Section III. 
Emphasis is placed on this treatment because an original formalism is developed 
which allows more physical insight into the energy transfer process when gases 
move through one another with a large relative velocity, a condition frequently 
occurring at high latitudes. The Appendix is a logical extension of Section 
III. Included there are proofs that the collision terms reduce to accepted 
expressions in the limit of zero relative flow velocity, and analytic calcula- 
tions of collision frequencies for a sufficient variety of interactions that 
most ionospheric requirements should be met. Because these calculations are 
detailed, extensive, and tend to divert attention from the principal objective 
(formulation of the conservation equations) , they are relegated to the Appendix. 

Formal development of the species conservation equations is completed 
with the collision term treatment in Section III. It is the task of Section 
IV to turn these formal equations into working equations, applicable to the E 
and F regions of the ionosphere. Here the three-fluid approximation is applied 
to obtain three sets of conservation equations. Through conservative order- 
of-magnitude estimates, a system of working equations of general applicability 
to ionospheric problems is obtained. Nonelastic and external source-loss 
processes, which are excluded from the formalism, but which nevertheless 
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significantly affect the conservation equations, are also examined briefly. 

The results of Section IV represent the starting points for practical problems 
in the sense that all approximations of general validity have been made. 

Further modifications will depend on the conditions of any particular problem. 

In Section V, transport properties of the gases in the ionosphere are 
discussed briefly in order to provide a means for remedying the closure problem 
of the set of conservation equations. It is characteristic of any set of 
equations formed by taking velocity moments of the Boltzmann equation that a 
higher velocity moment is introduced that is not determined by the (finite) 
set of equations (T.-Y. Wu , 1966). Through the use of transport coefficients 
(calculated or experimental), these higher order moments can be related to 
lower order moments, permitting closure of the set of equations. They can 
also be used in simplified model studies to examine some effects of transport 
without solving the entire set of velocity moment equations (e.g. , Schunk and 
Walker, 1971). Emphasis here is placed on electron transport properties, both 
because electrons are much more mobile than other gas particles due to their 
small mass and because an adequate treatment is readily available. Mathemati- 
cal detail is restricted primarily to that necessary to use results from the 
literature. The main exception is the derivation of an original method for 
approximating the energy dependence of collision frequencies by a power law in 
velocity, thus making available for arbitrary velocity dependence the rela- 
tively simple transport coefficient results of the power law case. Explicit 
analytic expressions for collision frequencies, some obtained from the litera- 
ture and some from the present investigation, are introduced at this point 
since they are basic to the transport properties. 

A problem of contemporary interest is examined in Section VI, making 
partial use of the mathematical apparatus developed in earlier sections. This 
is an investigation of the altitude structure of auroral E region neutral 
winds and some of the consequences of this structure for incoherent scatter 
radar observations. Since the primary calculation is based on observational 
data at a single location as a function of time, severe limitations are placed 
on the spatial information available. This results in a considerably simpli- 
fied calculation for a problem which is not necessarily susceptible to such 
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simplification. The use of theory, observation, and calculation together, 
however, makes it possible to obtain ^estimates of some elements left out of 
the calculation by necessity. A summation is presented in Section VII. 


1-16 


Section II 


KINETIC DESCRIPTION OF GOVERNING EQUATIONS 

2.1 GENERAL FORM OF THE TRANSPORT EQUATION 

Investigations of transport phenomena based on a particle description of 
gases generally begin with the Boltzmann equation. Since the derivation of 
the Boltzmann equation is discussed in detail in many texts on kinetic theory, 
it is sketched only briefly here. At time t, the number of particles in a 

-y ■+ —y —y 

volume element drdw of six dimensional phase space located at r, w is given by 

“V —y -)■ 

f(r, w, t)drdw, where f(r, w, t) is the distribution function of the gas. 
Suppose each gas particle is subject to an external force F and no collisions 
between particles occur. Then at time t + fit, particles previously at r, w 
at time t will be located at r + w fit, w + (F/m) fit, where m is the particle 
mass. Due to collisions, however, not all particles at r, w at time t will 
arrive at r + w fit, w + (F/m) fit at time t + fit. Likewise, due to collisions, 
not all particles at r + w fit, and w + (F/m) fit at 'time t + fit originated at 

—y -y 

r, w at time t. This collisional effect is expressed in the following 
equation : 


[f(r + w fit, w + ^ fit, t + fit) - f(r, w, t)]drdw 




drdw fit. 


( 2 - 1 ) 


Dividing through by drdw fit, expanding the first term to first order in fit, 
and taking the limit as fit 0 results in the Boltzmann equation: 


ff + w • Vf + - • V f = (|f) 
3t m w \^t/ c 


(2-2) 


The right side will be treated explicitly later; for the present it simply 
denotes the time rate of change of the distribution function due to collisions. 


If there is a mixture of gases, each species has a distribution function 
governed by equation (2-2). For the species r this is written 
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(2-3) 
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+ w 


Vf + — 
r m 

r 


/ 9f 

V w f r = U r 
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where the collision term now includes collisions with other species in -addi- 
tion to self collisions. The distribution function is taken to be normalized 
such that the number density of species r at the location r is given by 


» 

n r (r, t) = f r (^> w, t)dw r> 


(2-4) 


where the integration extends over all velocity space. In this definition, it 
is assumed that f^_ vanishes for infinite velocities. Conventionally, equatior. 
(2-4) is called the zeroth (velocity) moment of the distribution function. 


The velocity w^ denotes the total velocity of an r-particle with respect 
to some fixed (laboratory) reference frame. It is convenient to define also 
the species random (thermal) velocity by 


c e w - v 


(2-5) 


where v^ is the average velocity (also referred to as the drift or flow vel- 
ocity) of species r. In the literature random or thermal velocities are 
frequently defined with respect to the mass average velocity of the entire 
fluid. That is convenient when the species flow velocities of all species are 
close to one another. However, in the ionosphere species flow velocities may 
differ significantly from one another. A more consistent and physically more 
transparent treatment is possible if the parameters for each species are 
defined in their own reference frames. 


The species average value of some particle property 0^ is defined to be 


if , 

<e> e — 0 f dw 

r n r r r 

r 


(2-6) 


In particular 


v = <w> = — 

r r n 


w f dw 
r r r 


(2-7) 


It then follows from equations (2-5) and (2-7) that 
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= 0 . 


(2-8) 


<c> 


<w> 


V 

r 


A transport equation for the property 0^ is obtained by multiplying both sides 
of equation (2-3) by 0^ and integrating over all velocity space: 


r dr r r r _ > 

0 dw + 0 w • Vf dw + 0 — • V f dw 

J r 3t rjrr rrjrm wrr 


(2-9) 


It is assumed that the integrands are well-behaved and vanish as velocities 
become infinite. 


For macroscopic analysis, a more convenient form for equation (2-9) can 
be obtained by considering each term individually. If 0^ is assumed to have 
no explicit dependence on position or time, then since time and velocity are 
independent coordinates, the first term in equation (2-9) may be written 


[ 0 dw = [ 0 f dw E 7— (n < 0 > ) 

J r at r at J r r r at r r 


(2-10) 


from equation (2-6). The second term can be treated in a similar manner: 


0 w • Vf dw = V • 
r r r r 


0 w f dw e V • (n <w 9> ) 
r r r r r r 


( 2 - 11 ) 


The third term in (2-9) can be integrated by parts to obtain 


0 — * V f dw = - 

r m wrr 

r 


— v • (e ? )dw 
~ ” r r r 


m w 
r 


n r 

L < V -0 F> , 

m wrr 
r 

where the first equality follows from the assumption that f^ vanishes on the 
infinite velocity surface. For velocity- independent forces 


<v • 0 ^> =<? • V e> . (2-12) 

w r w r 

But this equation also holds for velocity dependent forces which are perpen- 
dicular to the velocity (e.g. qw x B/c). Since all velocity-dependent external 
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forces to be treated in this work are of this nature, equation (2-12) may be 
used, with the result 


f F 

f 0 — 

I r m 

* T 


V f dw = 
w r r 


Substituting equations (2-10) , 



V 0> 

w r 


(2-11), and (2-13) into (2-9) gives 


(2-13) 


— (n <0> ) + V • (n <0 w> ) 
3 1 r r r r 




v e> 

W I 


S r (ir) di5 r 


(2-14) 


This equation (or ones very similar) is variously referred to in the litera- 
ture as Maxwell* s transport equation (Sutton and Sherman, 1965), equation of 
change of molecular properties (Chapman and Cowling, 1970), or transfer equa- 
tion (Burgers, 1969). 


In its present form the transport equation is quite general (except for 
inelastic effects) and no improvement over the Boltzmann equation in terms of 
computational requirements, since the distribution function must first be 
known before the average quantities in equation (2-14) can be determined. The 
purpose of this formulation is to establish a bridge between the microscopic 
properties of the gas, which can be inferred from theory and specialized experi 
ments, and the macroscopic properties, which are observed under more general 
conditions. 


2.2 SPECIES CONSERVATION EQUATIONS FOR MULTICOMPONENT GASES 

Conservation equations are obtained by selecting for use in the general 
transport equation, ' those particle properties which are conserved in collisions 
In particular, mass, momentum, and energy are the conserved quantities of 
interest here. The object is to obtain equations governing the macroscopic 
observables of the gas, i.e., number density, flow velocity and temperature. 

In this section only terms on the left side of the transport equation (2-14) 
are treated in detail; treatment of the right side (collision terras) is defer- 
red to the next section. 
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2.2.1 Conservation of Mass (6 f = m f ) 

Since is simply a constant quantity, equation (2-14) can be written 
down immediately for 0^ = m^: 


3n + r / 3f \ 

m — ~r~ + mV%nv=m ( -r-—) dw , (2-15 

r 3t r r r r J \3t J ^ r 

where equation (2-7) was used. The common factor m^ can be cancelled out 
unless it is desired to work in terms of mass density. If the right side of 
equation (2-15) is denoted by p , this equation can be rewritten 


3n 
i 

at 


+ v 



(2-16) 


2.2.2 Conservation of Momentum (6 f - m r v^ f ) 

With 6 r = m r w r , equation (2-14) becomes 


3n 

r -* -► -± -+ 

n <w> + m V • n <w w> - n <F • V w> 

rat rr r rr wr 


= m 


s r (ir) 


dw = ? 


(2-17) 


where represents the total momentum transferred to the species r through 
collisions. In the second term, the tensor <w w>^ can be expanded, using 


equations (2-5) , (2-7) , and (2-8) : 


<w w> = <c c> + <w> <w> . (2-18) 

r r r r 

In keeping with the method of defining all macroscopic parameters for each 
species in the species frame of reference, the species pressure tensor is 
defined by 


p e n m <c c> (2-19) 

*r r r r 

and the species temperature by 

n k T e — n m <c^> = p , (2-20) 

r r 3 r r r r 

where p^ is the scalar pressure of species r. Equation (2-20) is seen to be 

the perfect gas law (for species r) and comparison with (2-19) shows that the 
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scalar pressure is the mean value of the diagonal elements of the pressure 
tensor. These definitions correspond to values which would be measured by 
appropriate gauges sensitive only to the species r and drifting with the 
species r at a velocity v f . Equation (2-18) can now be written 


<w w> 



(2-21) 


where equation (2-7) is used. The second term in (2-17) may then be written 


m V • n <w w> = V • p + m [n (v • V) v + v V • (n v ) ] , (2-22) 
r r r r v r r r r r r r 

where the last two terms follow from Chapman and Cowling (1970, p. 19). 


The third term of equation (2-17) follows directly from a standard vector 
identity to be 

n <? • V w> = n <?> 
r w r r r 

With these results equation (2-17) may be rewritten as 


m tt (n v ) + V • p + n m(v • V)v + m v V • (n v) 
r 3t r r r r r r r r r r r 

- n <?> = ? 

r r r 


(2-23) 


2.2.3 Conservation of Energy (0 r = % m r w f 2 ) 

2 

In the literature the thermal energy moment (1/2 m^ c ) is more fre- 
quently used in obtaining an energy equation than is the total kinetic energy 
2 

moment (1/2 m^_ w f ) used here. Since in those treatments, drift velocities 
are generally small compared with thermal velocities, the practical conse- 
quences of this distinction are negligible. In the ionosphere, however, flow 
velocities of ions and neutral particles can attain magnitudes of the same 
order as thermal velocities. It is therefore important to treat the energy 
contribution of both flow and thermal components. 

2 

For 0 r = 1/2 w^ , equation (2-14) becomes 
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(2-24) 


1 9 

2 m r 9t 


2 1 

(n <w > ) + — m V 
r r 2 r 


(n <w 2 w> ) - -r n <f • V w 2 > 
r r 2 r w r 



where E^ represents the rate at which energy is transferred to the species r 
by collisions with other species. Again, a more usable form is obtained from 
a term by term examination. 


In the first term the time derivative can be expanded: 

1 9 , . 2 . 1 2 8n r 1 3<w> r 

2 m r H ( V W V = 2 V w "r 9T + 2 m r n r 9t 


(2-25) 


Although, as an independent coordinate, w (and hence w ) is independent of t, 

2 r r 

this is not true for <w > , due to the time dependence of the distribution 

2 r 

function over which w is averaged. From equations (2-5), (2-7), and (2-8) 


2 2,2 

<w> =<c>+v 

r r 


(2-26) 


since the cross term vanishes. With the definition of species temperature, 
equation (2-20) , equation (2-26) can be written 


2 

<w > 

r 


3 k T 




(2-27) 


r 

This is equivalent to the statement that average total kinetic energy is the 
sum of the average thermal energy and the average flow energy. Equation 
(2-25) now becomes 


1 

2 


3 , 2 N 

m — (n <w > ) 
r 3t r r 


k T +ym v*l — ~ 
12 r 2 r r I 9t 


3 8T r 

+ ! n r k IF 


+ n m v 
r r r 


3v 

t 

at 


(2-28) 


In treating tlie second term of equation (2-24), it is convenient to expand 
first the averaged velocity expression. Making use of equations (2-5), (2-7), 
and (2-8) yields 
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(2-29) 


2 ~y 9 -> ->* -)■ 2v ■+■ . 

<w w> = <(c + 2c • <w> + <w>^ ) (c + <w> r ^ > r 

2 -*■ _“>■ *> ^ 2 -> 2 

= <c c> + 2v * <c c> + <c > v + v v 
r r r r r r r 

For species r the heat flux vector q r is defined as 

q = —■ n m <c^ c> , (2-30) 

^r 2 r r r 

representing the flow of thermal energy. With this and previously defined 
quantities, equation (2-29) can be substituted into the second term of 
equation (2-24), giving 

-j m r V • (n r <w 2 w> r ) = V • + V • (v^ • p f ) + -| k V • (n r T r v r ) 

+ \ m r V • (n^ v r 2 v r ) . (2-31) 

Most of these terms can be expanded readily in vector notation; however, 
the second term on the right is more easily manipulated in Cartesian tensor 
notation. Making use of the convention of summing over repeated indexes allows 
this term to be written 


(v • P ) - 


3x. 


(v j v 


O I u 

- V. p . . + P . . V. 

3 dx ± Zx ± 3 

After Chapman and Cowling (1970, p. 16) , the double product of two 
tensors is defined as 


(2-32) 

second rank 


A :B = A. . B. . . (2-33) 

ij Ji 

Then, since p^ is a symmetric tensor (see equation (2-19)), equation (2-32) can 
be written in the form 

V * (v % p ) = v *(V*p)+p:Vv . (2-34) 

r r r r r *r r 


Using this result and expanding other terms in equation (2-31) gives 


2-8 



~ m V * (n <w 2 w>)=V*q + v *(V*p)+p:Vv 

2 r r r n r r r r r r 

+ -| n kT V • v + -| n k v • VT + 4 k T v • Vn 
2r r r 2 r r r2 rr r 

1 2 -> -* _ N -> n 
+ — m v v * Vn + n mv •* [ (v * v)v ] 

2 r r r r rrrr r J 


1 2 + 
+ — n m v V • v 
2 r r r r 


(2-35) 


The third term in equation (2-24) is straightforward: 

F • V w 2 - 2F • w . (2-36) 

w 

It was previously noted that the only velocity dependence of F to be considered 
has the form of a cross product in velocity, for which the scalar product in 
equation (2-36) vanishes. Therefore, the third term in equation (2-24) may be 
written, using equation (2-7), 

- ^ n <F • V w 2 > = - n <F> • v . (2-37) 

2 r w r r r 

Combining results from equations (2-28), (2-35), and (2-37) permits equa- 
tion (2-24) to be placed in the following form: 


~ 3T 9v To i 

3 , r, -* r , 3 . _ , 1 2 r 

— n k ■— — I- n m v • — K — kT + — m v —— 

2 r 8t r r r at \_2 r 2 r r J at 

+ V • q + v • (V * p ) + p :V v +-^n kT V * v 

M r r *r v r 2 r r r 

+ n kv • VT + 4 k T v - Vn +^m v 2 v • Vn 
2 r r r 2 rr r2rrr r 

+ m r n r v r * [(v r * V) v^_ ] + \ n r m r v % v r 


- n <?> • v - E 

rrrr 


(2-38) 


Tlie energy equation in this form is rather cumbersome; and since further 
modifications will be made subsequently, a detailed description of the physical 
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significance is deferred. However, for future reference, attention is directed 
to the final term on the left side of equation (2-38). This represents work 
done on species r particles by external forces. In the ionosphere, if F repre- 
sents an electrostatic field, this term is called the "Joule heating" contri- 
bution (Cole, 1962). A generalization of this concept will be discussed when 
the conservation equations are placed in final form. 

2.2.4 Formal Modifications of the Conservation Equations 

Equations (2-16), (2-23), and (2-38), as a coupled set of time-dependent 
partial differential equations are inconvenient to use in calculations because 
the quantities of interest are coupled in the time derivatives of the equa- 
tions. A more workable form is to have the time derivative of only one of the 
macroscopic variables n^, v^, T^_ in each equation. Each equation then describes 
the time evolution of a single observable, although coupling remains present 
in other terms. Fortunately, obtaining this result is a straightforward oper- 
ation, and the results are well worth the effort. 


Equation (2-16), the continuity equation, already has the desired form. 
It is used to eliminate the time derivative of n^_ from the momentum equation 
(2-23). That time derivative can then be written 


9 -* 9 \ -+ pr ~1 

m r- (n v ) = m n -r~— + m v I V • (n v ) 

r 9t r r r r 9t r r I m r r r \ 


(2-39) 


Substitution of equation (2-39) into (2-23) gives 
9v 


, .. „ — + m n (v • V)v + V • p - n <F> = P - p v , (2-40) 

rr9t rrr r r r rrr 


m n 


where all collision terms are collected on the right side. 


For the energy equation both equations (2-16) and (2-40) are used to 
eliminate time derivatives of n^ and v . After considerable algebra, equation 
(2-38) becomes 
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(2-41) 


3T 

2 n r k dT + 2 n r k \ ' VT r + - V * *r + V V \ 

= E - * • ? + i \l n v 2 - 1 k T 1 . 
r r rm!2rr 2 rj 
r 

This represents a substantial simplification in the form of the energy equa- 
tion. One consequence of this simplification is the absence of an external 
force term, which for an electrostatic force was termed the "Joule heating" 
term in the previous section. It will be seen in the following section that 
this contribution is contained implicitly in the collision terms, where it 
physically belongs. 


Since the left sides of the conservation equations are now in final 
(and simplest) form, the physical meaning of the terms can be examined briefly. 
Examination of the collision terms is deferred until the next section. In the 
continuity equation (2-16) , the first term represents the explicit time rate 
of change of n^ , while the second represents the net change in n^ due to more 
or less r-particles flowing into a small volume than out. 

In the momentum equation (2-40) , the first term gives the explicit time 
rate of change of momentum of species r particles. The second term is the 
nonlinear term which gives the change of momentum of particles in a small cell 
flowing with the species r particles at velocity v^ due to the variation of v^ 
in space. The third term gives change in momentum due to internal forces 
within species r, that is, pressure and viscous forces. Finally, the fourth 
term gives momentum changes due to external forces, excluding collision (fric- 
tional) forces which are treated explicitly. 

The energy equation (2-41) is seen now to involve only thermal or random 
energy. All directed (mechanical) energy contributions may be recovered from 
the momentum equation which governs directed motion. The first term in equa- 
tion (2-41) gives the explicit time rate of change of temperature while the 
second gives the temperature change of a cell moving at velocity due to 
spatial differences in temperature. Net changes in thermal energy due to more 
or less thermal energy flowing into a small volume than out are given by the 


\ 
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third term. The fourth term gives the thermal energy of viscous dissipation, 
that is, the randomizing effect of collisions of r particles among themselves 
when the drift velocity varies in space. 

The net effect of the mathematical manipulations in this section has been 
the transformation from microscopic to macroscopic variables. The resulting 
equations separately describe the time evolution of number density, directed 
motion, and random motion, while remaining closely coupled. Although the 
physical situation has been improved in the sense that the equations can be 
interpreted in terms of physically measurable quantities, the mathematical 
problem is no less complicated. The terms treated to this point involve only 
the averaged interactions of species r particles with themselves and with 
external force fields. Interactions with other species is the object of the 
next section. 



Section ill 


formal Transfer integrals 

3.1 GENERAL FORM OF THE TRANSFER INTEGRAL 

The transfer integrals to be treated are those obtained from the right 
side of equation (2-14) when mass, momentum, and kinetic energy are succes- 
sively substituted for 0^. Initially, however, the form of the transfer 
integral for arbitrary 0^ is determined. For this an expression for (3f r /3t) c> 
the collision integral, is required. In the standard formulation of the 
collision integral, only elastic collisions are considered, and particle 
potentials are assumed to be spherically symmetric and of sufficiently short 
range that only binary collisions need be considered. This latter assumption 
is too restrictive to include Coulomb collisions. However, for electron-ion 
collisions use of the Boltzmann equation has been found to give the same 
results as the Fokker-Planck equation, which treats many simultaneous but 
independent collisions (cf Shkarofsky, et al. , 1966 or Burgers, 1969). So for 
convenience, the Boltzmann collision integral will also be used for electron- 
ion collisions. 

With these assumptions, the collision term in the Boltzmann equation (2- 
2) for species r is given by 

= 1 f [ f r C«;) f s (P) - f r (w r ) f s (w s )]g bdb -d* dw s , (3-1) 

' ' c s J 

where f = distribution function for species r,s 
r ,s 

w = total velocity of species r,s particle before collision 
r ,s 

w 1 = total velocity of species r,s particle after collision 

r ,s 

g e |w s - w | = relative speed 

b = impact parameter 
cf) = aximuthal angle. 

The summation extends over all species and the range of integration over all 
w^-space. A schematic diagram of the collision process is given in Figure 3-1. 
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The derivation of equation (3-1) is presented in many books (eg. Chapman and 
Cowling, 1970; Sutton and Sherman, 1965; Huang, 1963) and is not repeated 
here. For notational convenience, only collisions between species r and 
species s are considered for most of this section; appropriate summations are 
reinstated as needed. 


When equation (3-1) is substituted in the right side of equation (2-14) , 
the transfer integral has the form 




[ f * £* - f f ]g bdbd<f> dw 
r s r s s 



(3-2) 


where f 1 = f(w'). By means of symmetry arguments. Chapman and Cowling (1970, 
Chapter 3) show that an equivalent form of equation (3-2) is 



(0* - 0 ) f f g bdb d<j>dw dw , 
r r r s b r r s’ 


(3-3) 


where 0^ corresponds to post-collision values of the property 0^. This is the 
form of the transfer integral on which subsequent calculations are based. 


3.2 SPECIFIC TRANSFER INTEGRALS 

At this point, it is helpful to define carefully the plasma system to be 
encompassed by the model being developed. The ionospheric plasma of interest 
is the thermal plasma, i.e. the low energy plasma with generally random 
velocities. Plasma particles in this category have energies <2eV, with pre- 
dominantly Maxwellian distributions characterized by temperatures of a few 
hundred to a few thousand degrees Kelvin. This distinction is important 
because occupying the same space may be charged particles with energies of 
order lOeV (photoelectrons) , of order 1-10 KeV (auroral particles) and of 
order 100 KeV - 10 MeV (trapped radiation). In general, number densities of 
these higher energy particles are orders of magnitude smaller than that of the 
thermal particles in the ionosphere; and they do not have the approximately 
Maxwellian velocity distribution associated with the thermal particles. 

Because of these substant ia L Ly different properties, it is convenient and 
plausible to l reat the higher energy particles as separate entities and, in 
fact, to consider them as external to the system being treated. Although a 
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strictly consistent, complete treatment of the problem would require inclusion 
of all such elements in the system (e.g. S.T.Wu, 1970), considerable success 
in ionospheric investigations has been achieved by means of this distinction 
because it has a sound physical basis and it greatly simplifies a number of 
bookkeeping and computational problems. 


In addition to practical considerations, there is some theoretical justi- 
fication for taking this approach. From numerical solutions of Boltzmann- 
Fokker-Planck equations, Carleton (1968) and Ashihara and Takayanagi (1973, 
1974) have found that for electron energies less than about 1 eV, the energy 
distribution is very nearly Maxwellian; above those energies significant 
deviations from a Maxwellian distribution occur. Solving a similar equation 
for the atomic oxygen ion distribution function, Banks and Lewak (1968) deter- 
mined that for ionospheric conditions the ion distribution is only slightly 
distorted from Maxwellian in the steady-state nonequilibrium (T^T^T^) case. 


With these ideas in mind, discussion of the collision terms can proceed. 
For the case 0 r = m r > the ri ght side of (3-3) clearly vanishes, since at the 
low energies considered here mass is unaffected by an elastic collision. 

From equation (2-16), this implies p = 0. Here, however, recognition of the 
physical situation described above must be taken into account. This result 
is a direct consequence of ignoring inelastic collisions, chemical reactions 
and higher energy particles. Although these effects cannot be included in 
the present formal treatment of collisions, their place in the equations can 
be retained by simply retaining the p notation until details can be filled in 
appropriate to the specific conditions of a problem. Hereafter, p^/m^ is 
viewed as a source/sink term resulting from processes external to the system, 
in the sense described above. 


For 6 r = m r w r > equation (3-3) takes the form 




f f gbdb dd> dw dw 
r s & Y r s 


P 

rs 


(3-4) 


3-4 



Because of the presence of the relative speed g in the integrand, it is 
necessary to transform the velocities to the center-of-mass system. New 
velocity coordinates are the center-of-mass velocities. 


m w + m w 
r r s s 

m + m 
r s 




(3-5a) 


(right equality follows from conservation of momentum) and. the relative 
velocities 


g = w - w (3- 5b) 

r s 

g f = w* - w* . (3-5c) 

° r s 

Equations (3-5a, b, c) can be solved for and w & to obtain the inverse rela- 
tions 


w 




m 


* rs -> 
w = G - - — g, 
s m 

s 


(3-6) 


where u = m m / (m + m ) is the reduced mass, and similar relations hold for 
M rs r s r s 

post-collision velocities. From conservation of momentum and energy it can 
be shown that 


g = g ? , (3-7) 

so that only the direction of relative motion is changed by the collision. 

From these relations, the change in velocity of particle r due to the 
collision can be expressed as 

(*; - » r ) = (g* - g). 0-8) 

r 

The scattering angle X (see Figure 3-1) is given by 



g 


(3-9) 
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where equation (3-7) has been used. It follows that in the direction of g 


g* - g = -g (1 - COS X) . 


(3-10) 


Because of azimuthal symmetry in the scattering, other components of the vector 
(g 1 * g) » when averaged over many collisions, will not contribute to the inte- 
gral in equation (3-4) . That integral may therefore be written 

? = - u I gg (1 - cos x) ff b db d(j> dw dw . (3-11) • 

rs rs j r s r s 

The cross section for momentum transfer is defined to be 


* 

Q- C (g) = (1 - cos x) b db dij> , 

rs J 


(3-12) 


where the g dependence follows from the fact that x = x(g» b) in general. From 
the transformation equations it can be shown that 

->■ ->• 

dw dw = dg dG, 
r s ° * 


so that equation (3-11) can be written 

? rs “ ~ y rs \ 8 g Q(g) f r f s d & d5 » 


(3-13) 


where now f = f (g, 2) from equations (3-6). This is the form which is 
r , s r , s 

used for calculations with specified distribution functions. 


For 0^ = 1/2 m r w r > equation (3-3) becomes 


4- m f w 2 dw = 4 m [ (w* 2 - w 2 )g f f b db d(f dw 

2 r J r \3t / r 2 r J r r ° r s T r 


dw 


= E 


rs 


(3-14) 


From equations (3-6) , it follows that 
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(3-15) 


2 m r (w r 2 " w r* = p r s & “ g) • 

As previously, averaged over many collisions, the vector (g* - g) contributes 
to the integral only in the component parallel to g. Making use of equations 
(3-10), (3-12), and (3-15) permits equation (3-14) to be written 

E rs * " y rs f g ® ^ Q(g) f r f s dg * <3-16) 

In this form the energy transfer integral is evaluated for specific distribu- 
tion functions. 

3.3 TRANSFER INTEGRALS FOR MAXWELLIAN DISTRIBUTION FUNCTIONS 
3.3.1 Background 

Before a lengthy calculation of this nature is begun, it is appropriate to 
question the need and the physical and mathematical justification for it and 
to review previous efforts in this direction. In order to evaluate the momentum 
and energy transfer integrals, and E rs > some form for the distribution 

functions f^ and f g must be used. The simplest, physically reasonable choices 
are two Maxwellian distributions with a common temperature and zero average 
velocity; reasons for choosing a Maxwellian distribution were discussed pre- 
viously. In this case, however, application of energy conservation in equa- 
tions (3-1) or (3-2) shows that the collision terms vanish. Moreover, such 
conditions do not in general hold in the E and F regions of the inosphere. It 
has been well established from both ground-based and in situ measurements that 
electron and neutral temperatures differ up to several thousand degrees Kelvin 

with ion temperatures intermediate to these. Flow velocity differences, while 

2-1 5-1 

generally small (< 10 cm s ) , may reach magnitudes of ~10 cm s under 

certain conditions; this is of the order of ion and neutral thermal velocities. 

Burgers (1969) has treated the transfer equations in considerable detail 
and considered these conditions specifically. He refers to the situation of 
small relative differences in temperatures and flow velocities as | the fully 
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linear case. Large temperature differences but small differences in flow 
velocities constitute the semilinear case. By extension, unrestricted differ- 
ences in both temperatures and flow velocities are termed here the nonlinear 
case. Although this terminology is not widely used, partially because the non- 
linear case has received little attention, it aptly describes the range of 
approximations for this aspect of ionosphere modeling. 

A completely rigorous treatment should also consider distortions of the 
Maxwellian distribution, for example, by use of a Chapman-Enskog method (cf 
Burgers, 1969). This is not attempted here because the calculational require- 
ments are large and the degree of accuracy to which basic data, such as momen- 
tum transfer cross sections, are known does not presently warrant such an 
effort. Moreover, as noted previously, theoretical determinations of distri- 
bution functions for ionospheric conditions have found deviations from the 
Maxwellian distributions to be small for the thermal energies considered here. 
The most important contributions to the transfer integrals are expected to be 
due to temperature and flow velocity differences, rather than distortions of 
the distribution functions. 

A variation of the semilinear approximation considers arbitrary tempera- 
ture differences, but assumes static gases. This problem has been treated in 
detail for Maxwellian distribution functions by Banks (1966 a, c) among others. 
It can be considered a limiting case for the nonlinear problem and is so used 
later in examining results. Treatments of flow velocities (i.e. displaced 
Maxwellian distributions) have not been so complete in either the semilinear 
or nonlinear approximations. 

Morse (1963) has made analytic calculations in the semilinear approxima- 
tion, presenting both general expressions (arbitrary interaction potential) 
and specific results for Maxwell, Coulomb, and hard sphere interactions for 

both momentum and energy transfer. However, his energy transfer integral is 

2 ' 

based on the thermal (random) energy moment (1/2 me ) rather than on the total 

2 

energy momeint (1/2 mw - see equation (3-14)) considered here. Noting this 
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distinction. Ingold (1968) has presented specific analytic results in the 
nonlinear approximation for Maxwell and Coulomb interactions. ^ However, he has 
presented no general equations, except in the semilinear approximation. 

Apparently without knowledge of these other studies, Stubbe (1968) has 
also investigated this problem. He has treated only the momentum transfer 
integral, obtaining a general result in series form, which can readily be shown 
to reduce Morse’s (1963) closed form integral. His specific results are 
obtained numerically and apply to charge transfer interactions. Burgers (1969) 
has presented a comprehensive, systematic development of the fluid dynamic 
equations from the Boltzmann equation. However, although he has laid out the 
Chapman-Enskog procedure for treating the nonlinear problem, he has not 
attempted to solve it. Rather his treatment of this problem extends only to 
hard sphere interactions, using displaced Maxwellian distribution functions. 

Recently Horwitz and Banks (1973), apparently independently of the others 
cited, have applied the consideration of large flow velocity differences to 
the ionosphere for both momentum and energy transfer in the nonlinear approxi- 
mation. However , they provide no general expressions ; and their results , 
although pertaining to charge transfer collisions, are in an approximation 
corresponding to a hard sphere interaction. As with Burgers’ (1969) similar 
results, the energy transfer calculation is based on the thermal energy moment 
rather than the total energy moment. 

A notable feature of this brief review is the number of times this type 
of problem has been investigated, generally independently, in a piecemeal 
manner. It indicates an interest in the problem for a variety of reasons. For 
ionospheric research in particular, there is good reason for interest. Much 
attention is currently being directed towards high latitudes where large elec- 
trostatic fields are observed driving the ionospheric plasma at high velocities. 
The previous partial treatments have been found to be inadequate for the pre- 
sent study, particularly with respect to the energy transfer integral. So 
one of the objectives of this investigation is to present a complete formal 
treatment of the nonlinear problem, in the approximation of two displaced 
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Maxwellian distribution functions, and to provide a set of specific results 
which should be adequate for most requirements. Below, general formal expres- 
sions are developed for arbitrary interactions. In the Appendix, the formal 
results are applied to several specific interactions to obtain analytic expres- 
sions suitable for numerical evaluation. 


3.3.2 Momentum Transfer 

Momentum transferred to species r by species s through elastic collisions 
is calculated from equation (3-13). Both distribution functions are assumed to 
be displaced Maxwellians, given by 


/ m \ 3/2 

"'W jj 20H S r • V 2 J 


(3-17) 


( " e \ 3/2 r " 


2 M - V 

s 




where parameters v , v , T , T have arbitrary values. Placing these expres- 

r S IT 3 

slons in equation (3-13) gives 

,3/2 


f - - y 
rs rs 


/mm \ 1 f 
” 1 8 8 Q ~ <e) 

£ 2kf; <S r - ? r )2 - 2kf; <3 s - *. )2 J d * iZ (3 ‘ 18 > 


The exponent of the exponential function must now be transformed to the center- 
of-mass system, using equations (3-6). When this is accomplished, the expo- 
nent, denoted by X(g, 5) , can be written in the form 


X(g, <$) - -(Ag 2 + BG 2 + Cg • £ + Da*g + E?"£+F) ( 


where 


A = 


rs 

2k 


[ m T + m T ~| 

r r ss I 

m m T T 
r s r s J 


(3-19) 

(3-20a) 
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(3-20b) 


B = 


m T + m T 
r s s r 

2k T T 
r s 


D = - 


E = - 


rs 

rv 

k 

L T r 

V 

'rs 

kT T 
r s 


1 


kT T 
r s 


, a = ( v T - v T ) 
r s s r . 


r r s s s r 


2 2 
m v m v 

F = r - r - + -g-§- 
“ 2kT 2kT 

r s 


(3-20c) 


(3-20d) 


(3-20e) 


(3-20f ) 


Since the only unknown in equation (3-18) is the momentum transfer cross 
section which is a function of g only, the integration over 2-space can be 
carried out immediately. For notational convenience, the integral in equation 
(3-18) is denoted by t p , which can be rewritten as 

l p = j dg g g Q rg (g) exp[-(Ag 2 + D ct • g + F)] 

• | exp [-BG 2 - 5 • (Cg + Eg)] d£ . 

With the notation 

t = eg + e| , 


(3-21) 


(3-22) 


the integration over 2-space can be performed by completing the square, giving 

3/2 

| expt-BG 2 - 2 . £]d2 = exp(S 2 /4B). (3-23) 

When this result is substituted back into equation (3-21), together with 
definition (3-22), may be expressed as 
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V 


- (t) 


3/2 


g g Q(g) exp[-(A' g 2 + y ' • g + C' ) ]dg 


(3-24) 


where 


A E 


A - C 2 p rs 


4B 


2kT 


(3-25a) 


y * = Dot - 


CEB 

2B 


rs v 

kT~ (v r - V 


(3-25b) 


C' = F — 


2 2 u „ 

E 3 rs /->■ ^ 2 

“4B~ = kT (v r - v s } 


(3-25c) 


The right sides of equations ( 3—25) follow from equations (3-20) and the defi- 
nition of a center- of -mass temperature 


“T T "| 

T E y — + — I . 
r rs m m I 

L r s J 


(3-26a) 


For convenience the following notation is also introduced: 


v E v - v 
o r s 


K E y /2kT 
M rs 


(3-26b) 

(3-26c) 


Equation (3-24) may then be expressed concisely as 


f p - (t )' 


3/2 


g I (g) exp[-K(g - v ) 2 ] dg . 
r s o 


(3-27) 


If a spherical coordinate system is selected for g-space, the angular 
integration can be performed immediately. This can be accomplished by means 
of symmetry arguments (e.g. Mathews and Walker, 1965 ^ Chapter 3). Since v q is 
the only direction specified in the integral, must be parallel or anti- 
parallel to it; this is specified mathematically by 


f p = A ^ , 

P o o 


(3-28) 
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where A q is a proportionality constant to be determined. Let v q designate the 
polar direction in a spherical coordinate system in g-space. Then from equa- 
tions (3-27) and (3-28) 

v • I p -j . .3/2 f 9 2 

A Q = = — 2 (g) v Q g COS0 Q rs (g) exp[-K(g + v q - 2v Q g cos0)] 

o o 

• g 2 dg sinBd0d<j> . (3-29) 


Because of azimuthal symmetry, the variable 4> may be integrated out directly 
and the 0 integral evaluated. For that part of the integral involving only 
angles, the result is 
1 

2tt | cose exp (2KgV Q cose) d cos 0 
-1 


4tt 

2K g V o 


cosh(2KgV Q ) - 2 Kg V ~‘ s inh(2Kgv o ) 


Placing this result in equation (3-29) gives 


A = 
o y 


/ . T / \3/2 exp (-Kv ) 

47TkT / 7T \ . o 

U/ 2 

rs v 


Q(g)exp(-Kg ) 


[ 


sinh(2Kgv 


cosh(2Kgv o ) - - 2Kgv 


■'] 


° i 

g dg . 


(3-30) 


With the relation 


d 

dg 


sinh(2Kgv ) 

, 2K S V o ° - 


1 

g 


cosh(2KgV Q ) 


sinh(2Kgv ) 
2KgV o _ 


(3-31) 


equation (3-30) can be modified and substituted into equations (3-28) , with the 
result 


V 


4nkT 


rs 


(i) 


3/2 exp (-Kv ) 


Q rs ( s) 


exp (-Kg 2 ) g^d 


[ sinh(2Kgv o )~j 

-^T J- 


3-13 


(3-32) 


Since tp represents the integral in equation (3-18) , that equation may 


now be written 


rs 


exp(-Kv‘ 


-n n 


rs 


r s 2 

JttK v 

f o 


} ” 
n J 


2 4 I sinh(2Kgv o )' 

o , Q rs (g)exp(-Kg )g d| ~ 2Kgv” 


|“ s inh ( 2KgV Q ) ~j 

re J' 


(3-33) 


An effective collision frequency for momentum transfer is now defined by 
the relation (Stubbe, 1968) 


' E -n v y v 
rs r rs rs o 


(3-34) 


Comparison with equation (3-33) identifies the collision frequency as 


v (v ) 
rs o 


= 2 


/r f . , , , „. 2 , 4,i slnh(2K *v 


2 

/ 

o 


j Q rs (g)exp(-Kg 2 )g 4 d l 2Kgv 


[ sinh(2Kgv o )“| 

~^ r o J 


(3-35) 


Equation (3-35) formally generalizes the collision frequency obtained by 
Banks (1966a) for two gases at different temperatures to include the effects 
of average relative motion between the gases. This equation is equivalent to 
the expression obtained by Morse (1963). It is also equivalent to the series 
form obtained by Stubbe (1968) , as may be verified by expanding the hyperbolic 
sine in equation (3-35) in a power series. In the Appendix it is shown that 
in the limit 0, equation (3-35) has the form 

oo 

v rs (0) - | n s /|K 2 | Q rs (g) exp (-Kg 2 ) g 5 dg , (3-36) 

o 

which is the expression obtained by Banks (1966a) . Equation (3-35) is evalu- 
ated analytically in the Appendix for a number of forms of the momentum trans- 
fer cross section, corresponding to collisional interactions of general in- 
terest. Specific collision frequencies of aeronomical interest are compiled 
in Section V. 
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3.3.3 Energy Transfer 

Total energy transferred to species r by species s through elastic colli- 
sions is calculated from equation (3-16). When displaced Maxwellians, equa- 
tions (3-17), are assumed for the distribution functions, equation (3-16) may 
be written 


m m 

„ f r s 

E = - y n n I r—r- 

rs rs r s l , 2, 2 


,3/2 


4tt k T T 
r s, 




g Q rs (g)exp[-(Ag + Da • g + F]g 


• | £ exp [— (BG 2 + 2 • (Cg + E(£)]d£, (3-37) 

where the constants are as defined in equations (3-20), Denote the integral 
over (5-space by Since the exponent in is the same as previously, 

■fci Cj 

equation (3-22) can be used to write 

2 

(3-38) 

If variables are changed to £ = 2 + ?/2B, equation (3-38) becomes 

(3-39) 


t £ = exp(S 2 /4B) (5 exp [-B ((5 + — ) ]d<5 


t = exp(S 2 /4B) 


(t " |g) exp(-B£ 2 )d£ . 


The integral of the first term vanishes since the integrand is an odd function 
of each cartesian component of J over the symmetric, infinite range. Integra- 
tion of the second term then gives 


t /,- 3/2 


t E - - exp(S 2 /4B) (f) 
by equation (3-23). 


(3-40) 


With this result equation (3-37) becomes 


V5/2 2ttk T T 


E = n n y 
rs r s rs \ it 


IT K. I 1 f 

r s > 

mm I ® 

r s J 


[Cg + Efl]gQ (g) 

la 


• exp [-K(g - v q ) ]dg . 


(3-41) 
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Designate the first and second terms on the right by 1^ and respectively 
and consider first. The integral in is basically the same as I p ; so 
from equations (3-32) through (3-35) can be evaluated immediately: 

I 9 =* n n rS 2 irk 2 T T ) E | •( g g Q (g)exp[-K(g - v ) 2 ]dg 
2 rsmm r s \tt/ & & x rs r ° o ° 

r s J 




rc (v ) » 

rs o 


where E has been evaluated from equation (3-20) . 


(3-42) 


For 1^ the integration over angles can be carried out directly: 

| sinh(2Kgv ) 

2tt I exp[2Kgv Q cos0]d cos0 = 2ir Kg~v — ~ 

-1 ° 


(3-43) 


1^ is then given by 


/ ‘"rs \ /K\ 3/2 exp(-Kv o ) r 4 

-n n I -r ) (~ ) 4iik(T - T ) g Q (g) 

r slm r + m / \tt / r s v q J 6 rs 6 


exp(-Kg )sinh(2Kgv Q )dg . 


(3-44) 


By differentiating with respect to v , the following equation can be 
verified from equations (3-31) and (3-35) : 


^ g^ Q(g) exp(-Kg 2 )sinh (2Kgv Q )dg 


= /¥{ 2 K tv o ex P (Kv o )v rs (v o )] + 2K exp(Kv 2 )v rs (v o ) 


(3-45) 
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When this equation is substituted into equation (3-44) and part of the dif- 
ferentiation carried out, the result can be written 




rs + (3kT + y v^)v 

rs o rs 


] 


(3-46) 


With equations (3-42) and (3-46) , equation (3-41) can be written 

3v 


rs 




rs 


-> 

? 3 • v 

+ [3kT + y v + 


rs o T - T 

r s j 


kT * 
o 8v 


-]■ 


(3-47) 


This is the desired result. 


Equation (3-47) expresses the energy gained by species s due to elastic 
collisions in terms of the effective collision frequency for momentum transfer. 
This formulation in terms of momentum transfer collision frequency is original, 
and it is useful in simplifying both calculations and interpretation of colli- 
sional energy transfer. Expressed in this way, terms which depend on tempera- 
ture differences or average velocity differences or both are clearly separated. 
The high velocity correction (Sv^/Sv^) term is an isolated term as opposed to 
a set of factors (e.g. , Horwitz and Banks, 1973). This permits a more direct 
assessment of the importance of the correction involved. Once the collision 
frequency is determined, the additional term is found by differentiation, 
which is simpler than performing the additional integral required by other 
formulations. This is illustrated in detail by the analytical evaluations in 
the Appendix. 

Consideration of some special cases aids interpretation of equation (3-47) . 

If the species flow velocities are the same (v = v ) , then v = 0. In the 
r r s o 

Appendix it is demonstrated that in this limit ^ v rs (^)/^ v 0 0- Equation (3-47) 

then has the limiting form 

z = - n (-— ■ I s —) 3k (T - T ) V (0), (3-48) 

rs r\m+m/ r srs 

\ r s/ 
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where v rg (0) is given by equation (3-36). This is the result derived by Banks 
(1966a) for these conditions. It can therefore be concluded that the term 
3kTv rg in brackets in equation (3-47) represents the transfer of random energy 
(heat) from one gas to the other due to temperature differences. 


For = T g the energy transfer given by equation (3-47) is 

E = n (m + m ) I — v^ - v • v I v (v ) , 

rs r r slm^ o r olrs o 


(3-49) 


where the definitions of T, and v q (equations (3-20a) , (3-26a) , (3-26b) 

have been used. Comparison with corresponding results of Morse (1963) shows 

that the second term in equation (3-49) constitutes the difference between 

taking the total energy moment and the random (thermal) energy moment of the 

collision integral. This term is seen to represent the rate of work done on 

species r by the frictional force - (m + m )v v due to the difference in 
r J r s rs o 

flow velocities. Thus the term (g • v Q )/( T r - T g ) in equation (3-47), which 
is the right side of equation (3-49), is seen to have two components, one 
representing a transfer of directed energy and the other a conversion of 
directed energy to random energy. The collision of two monoenergetic beams is 
a limiting example of this case. 


The remaining terms in brackets in equation (3-47), v^kT Sv^/Sv^ anc * 

p v^v , correspond to cross terms in the sense that both conditions T 4 T 
rs o rs • r s 

and v^_ 4 v g must be satisfied for a nonvanishing contribution. The first of 
these depends on the interaction potential for its sign as well as its magni- 
tude. Neither term has a straightforward physical interpretation, due to the 
dual requirement for a finite value. 


In the energy conservation equation (2-41), the elastic collision terms 
are (E^ - v^ • ? r ) . The first term is the rate at which total kinetic energy 
is transferred to species r through collisions. From this is subtracted 
(second term) the rate at which work is done on species r by collisions, that 
is, the rate of transfer of directed energy. What remains is the random energy 
part, which is consistent with the left side of equation (2-41). The collision 
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V 


terms can be evaluated in terms of the momentum transfer collision frequency 
from equations (3-34) and (3-47) (with a summation over all species). The 
result, after lengthy but straightforward algebra, is 


rs 


. t 

r rs 


n ij 
r M rs 

m + m 
r s 


m 


m v v - k(T - 
° o rs r 


T )(% + v 

s y rs o 3 v q J 


(3-50) 


Since, from discussion above, the right side of equation (3-50) must repre- 
sent thermal energy, the physical meaning of the terms can be obtained. In the 
case of thermal equilibrium (T^ = T g ) , but v q 4 0, only the first term remains. 
This term thus represents the conversion of directed energy to thermal energy, 
which is a generalization of the Joule heating usually associated with electric 
fields. In this form, Joule heating (or frictional heating) is seen to result 
from any forces which cause the gases r and s to move with different average 
velocities; and this difference in velocities may be due to the nature of the 
force (for example, electric field) or due to conditions of the medium (for 
example, a third species which couples preferentially with species r or s) . 

From the form of equation (3-50), it can be inferred that for a two-component 

2 

system, the total rate of Joule heating is n u v v with a fraction 

° r M rs o rs 

m s /(m r 4- m g ) going to species r (Banks and Kockarts, 1973, Chapter 22). 

If v q = 0 but T^ 4 T , only the second term in equation (3-50) remains. 

This gives the rate of transfer of thermal energy to species r due to the 
temperature difference. For the third term not to vanish three conditions 
must be met: (1) T^ 4 T g , (2) v^ 4 v g , (3) v rs must vary with velocity. This 

term represents a correction term for large relative average velocities, but 
only if the gases are in thermal nonequilibrium, and only if v rs is velocity 
dependent. 


3.4 FINAL FORM FOR FORMAL CONSERVATION EQUATIONS FOR MULTICOMPONENT CASES 

Results of Section II for the left sides and Section III for the right 
sides of the conservation equations can now be combined to obtain the final 
forms. For multicomponent gases the species conservation equations are: 
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(i) Continuity Equation (equation (2-16)) 


3n r - p r 

+ V • (n v ) - — 

at r r J m r 

(ii) Momentum equation (equations (2-40) and (3-34)) 


3v 


at 


— + (v_ • V)v_ + ~~ZT V • p - < ^r > 


r n m 
r r 


- I 


rs 


(v - V ) V 


m rs r s n m r 
s r r r 


(3-51) 


(3-52) 


(iii) Thermal energy equation (equations (2-41) and (3-50)) 

3 9T r 3 -► -* = -*■ 

— n k — h — n kv • V T + V • q + p:Vv 

2 r at 2 r r r n r *r r 




rs 


m -t- m 
s \ r s) 


m v v - k (T - T ) 
s rs rs r s 


(»- + ~ £)] 


5T p >Vr - 3/2 k I r J 


(3-53) 


where the notation v = v - v has replaced the v previously used. 

rs r s r o r 


These equations are formal in the sense that specific interactions must 
be known before can be computed from equation (3-35) and must be speci- 

fied. In addition, there is a closure problem for the system of equations. 
Although these five equations are sufficient to solve for n^, v , and T^, the 
additional unknowns associated with p^ and q^ require considerably more infor- 
mation than is given in these equations. And finally, as with such partial 
differential equations, initial and boundary conditions must be specified. 

Part of the additional information required can be determined from general 
ionospheric conditions, and some simplification of the equations can be obtained 
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by considering such. These matters are treated in the next section. Some in- 
sight can be obtained from investigations of transport properties under con- 
trolled conditions, obtaining results which can be applied more generally. 

This aspect is examined in Section V. Finally, certain information must await 
the particular problem to be solved for specification. This latter category 
is by its nature a catchall, since every variable which has not been tied down 
by previous efforts must be approximated or specified by assumption before the 
problem can be solved. Section VI illustrates this procedure in a numerical 
calculation. 
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Section IV 


if 


THREE-FLUID IONOSPHERIC MODEL 

4.1 CONSERVATION EQUATIONS IN THE THREE-FLUID APPROXIMATION 

In the three-fluid approximation, macroscopic properties of the electrons, 
ions, and neutrals are treated separately, while differences of these proper- 
ties among the various ion or neutral species are assumed negligible. Specifi- 
cally, all ion species are assumed to have a common temperature and flow 
velocity v_^ , and all neturals are similarly assumed to have common temperature 
and flow velocity v^. With these assumptions, discussed in Section I, the 
relevant conservation equations for ions and neutrals can be obtained by 
summing over the species comprising these gases. Electron equations are 
given directly by the species equations since a single species constitutes the 
gas. 


From equation (3-51) , the electron continuity 


UdLlUU 


3N 
e 

at 


+ v 


(N v ) = -S 
e e m 


(4-1) 


For number densities, the notation N will hereafter denote the total number 
density of the (ion or neutral) gas, while n will denote a species number 
density. These are then related by 


N - l n, • 
j 


(4-2) 


To avoid cumbersome notation, subscript i will denote ion and n, neutral; 
summations will imply summation over all relevant ion or neutral species, as 
specified by an i and/or an n. Thus, these subscripts serve the double 
function of summation index and gas indicator. In practice, this should cause 
less confusion than a proliferation of subscripts. Summing equation (3-51) 
over all ion species gives 


3n. 

1 X 


(n i V 
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From charge neutrality, * N e » 80 that interchanging the order of summation 
and differentiation allows this equation to be rewritten as 


8N 
e 

3t 


+ n Pi 

+ v . (n v , ) = y — 

e i h m. 


(4-3) 


where v.^ was factored out of the summation. In a similar manner the neutral 
continuity equation is found to be 


8N 


— - + v • (n v ) = y — 

3t n n L m 


(4-4) 

n n 

If a problem involves composition change as one of the significant physical 
processes, this formulation is inadequate; a separate continuity equation 
(3-51) must then be used for each ion and neutral species 


Equations of motion are obtained from equation (3-52). For electrons 
the equation may be written down directly: 

3v i <f > 

e , / -> r-rX"*' 1 1 r-» © 

- + (v • V)v + ~ — — V • p — 

dt e e N m e m 

e e e 


s=i ,n e \ e/ 


(4-5) 


where the approximation y 


es 


m has been used for the reduced mass of the 
e 


electron with any ion or neutral particle. To obtain the momentum equation 
for ions and neutrals, it is convenient to multiply equation (3-52) by the 
number density and mass before summing. The ion equation of motion may then 
be written 


V n. m. — — + y n. m. (v. • v)v. + J V • p. - J n , <t . > 

£ l l 3t £ l i l l £ r i £ i l 

= - Y n. y. v. (v. - v ) - v. V p 
L i is is l s l £ K i 

i,s=e,n i 


(4-6) 


Let the mean mass per particle m of a gas (ion or neutral) and the mass density 
p be defined by 


= N m E y n m 
^ s s 


(4-7) 
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Similarly, let the mean effective collision frequency for momentum transfer 
between gas r and gas s (if r represents ions, then s represents neutrals and 
vice versa) be defined by 


= — y n p V . (4-8) 

rs p L r rs rs 
K r r,s 

With these definitions and the fact that the three-fluid approximation allows 
the flow velocities to be factored out of the summations, equation (4-6) can 
be rewritten as 


3v 

il 


1 + (v ± • V)v. + 77 v • I P ± - 77 I v*i 
i i K i 


(4-9) 


It 

P i i 

Following the same procedure for neutrals gives the neutral equation of motion 


v . 

= - v\ ( v - v ) - v". ( v . v ) Y p . 

mi n le l e p . h l 


3v .. ^ 

t — + (v • v) v + — v • y p - — y n <F > 
at n n P L r n o L n n 
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= - v 


n n 


(v - v.) - v (v - v ) — Y p 

n l ne n e p L M n 


ni 


(4-10) 


n n 


It is helpful at this point to establish the relation between v rg and 

v - The total momentum transferred to gas r due to collisions with gas s is 
s r 

given by (cf equations (3-34) and (4-6)) 


y n y v (v - v ) = p V (v-v) 

^ r M rs rs r s F r rs r s 
r,s 

by equation (4-8). From equation (3-35) it can be seen that 


(4-11) 


n v = n v 
r rs s sr 


(4-12) 


When equation (4-12) is used on the left side of (4-11), the result obtained 
is 


) n u \> -pv = o v 
L s M rs sr H s sr rs 


(4-13) 
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The distinction between equations (4-12) and (4-13) is sometimes confused. It 
can be seen that the proper relation depends on whether the momentum transfer 
term in the momentum equation is written as in equation (4-6) (use equation 
(4-12)) or as in equation (4-10) (use equation (4-13)). Consistent use of the 
notation Introduced above should eliminate any confusion in the remainder of 
this study. Application of equations (4-13) and (4-8) gives the following 
specific relations: 
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(4-l4a) 


(4-l4b) 


(4-14c) 


This allows all collision terms to be expressed in terms of the three collision 
frequencies v e ^, and v in or, alternatively, their averaged counterparts v 

v en* '’in' 


The energy equations are obtained from equation (3-53) . For electrons 
this can be written 


3T 

4 N k + ^N kv •• V T + V • q +p:7v 
2 e 3t 2 e e e e e e 
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m 12 e e 2 el 


(4-15) 


where v = v - v . Equations for ions and neutrals follow from proce- 
rs 1 r s 1 

dures used above. For ions the result is 
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where equation (4-12) and the fact that v^ s * v gr have been used freely. The 
corresponding equation for the neutral gas is 

3 X 

~ N k + 4** k v • V T + V • J q + 7 p :V v 
2 n 3t 2n n n n n r n n 
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+ 7 v z y p - 4 k T y— . (4-17) 

2n^^n2n^m 
n n n 

This completes the set of conservation equations in the three-fluid approximation. 


4.2 IONOSPHERIC CONDITIONS AND CORRESPONDING APPROXIMATIONS 

Representative ionospheric conditions were indicated in Section I, Figures 
1-1 through 1-4. Here the Interest is not so much in typical conditions as in 
the extremes likely to be encountered in a wide variety of potential applica- 
tion's. The objective is to strip the conservation equations of unnecessary 
terms without limiting the scope of applications. This is accomplished by con- 
sidering the range of values of variables and parameters likely to be encoun- 
tered under usual and unusual conditions, then, from order- of-magnitude esti- 
mates, discarding those terms of negligible contribution thrpughout the range. 
At this stage a conservative approach is taken — borderline or questionable 
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terms are retained. There is ample opportunity for further reduction of the 
equations when a particular problem is defined and physical conditions are 
known more precisely. 

Extremes of ionospheric variables generally correspond to different 
altitudes, different geographic locations, and different times. Values given 
are based on observations or inferences from observations and should be taken 
as indications of order of magnitude rather than as absolute limits. With 
that caution, extreme values are presented in Table 4-1. These values will be 
used as the basis for order-of-magnitude estimates in the remainder of this 
section. 


Table 4-1. EXTREME VALUES FOR SELECTED IONOSPHERIC PARAMETERS 


Masses : 

Vv 10 ~ 4 ’ 

s = i , n 

Number Densities: 

N e /N n s 10- 3 


Flow Velocities*: 

10 2 .< v < lo 5 

cm s' 1 

Temperatures: 

10 2 < T < 10 4 

°K 

Spatial Scales*(L) 

10 6 < L < 10 8 

cm 

Time Scales*: 

T > 1 s 


Electric Fields: 

0 < |t| < 150 mV m' 1 


^ Lower limit is a practical limit of interest and observability . 
* Restriction - not a natural extreme . 


4.2.1 Continuity Equation 

As noted in Section III, the term p/m on the right side of the continuity 
equation is a phenomenological term in this study because only elastic colli- 
sions are treated systematically. It is appropriate at this point to discuss 
this term and its role in the ionosphere more specifically. Rather than a 
single term it actually represents the net sum of all processes which add or 
remove particles of the appropriate type to or from a given volume element. 
This can be represented mathematically by 
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p s ,m s " l (P - L ) , 

s v s s' 


s = i, n, e 


(4-18) 


where P denotes particle production processes and L particle loss processes. 

For ionization, the electron and ion processes are the same; and since 
only singly ionized ions are present, ions and electrons are produced and lost 
at the same rate. In the notation of equation (4-18), this is written 


p - y p. , l = y l, . 

e £ i * e j i 


(4-19) 


The "P" terms originate primarily from photoionization of neutral atoms and 

o 

molecules by extreme ultraviolet radiation (A < 1029 A for primary constitu- 
ents). An extensive set of chemical kinetic processes modifies the ion compo- 
sition from that due to photoionization. Photoelectrons produced in this way 
are not the thermal electrons of interest here, but since they quickly become 
thermalized, the photoionization rate is generally taken to be the thermal 
electron production rate. Collisional ionization is another source of thermal 
plasma, but its importance is restricted primarily to high latitude regions 
where large numbers of relatively high energy (>1 keV) charged particles fre- 
quently penetrate to the altitudes of interest here (90 km to 800 km) . 


The important loss process for charged particles is dissociative recombi- 
nation, represented schematically by 

X Y + + e -► X* + Y* 

where the asterisk denotes a possible excited state. Of secondary significance 
is radiative recombination 

X + + e ->• X* ->■ X + hv 

which may play a significant role only in the upper F region where atomic 
neutral species predominate. Three-body recombination is important only below 
the E region where three-body collisions can occur frequently enough to be 
significant. 
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Despite the apparent coupling between the source-loss terms of the 
charged particle and neutral particle continuity equations, they are indepen- 
dent, except for possible coupling in the transport terms. This results from 
the low degree of ionization, particularly at altitudes where production and 
loss rates are highest (110 km to 250 km); that is, the number of particles 
involved in ionization and recombination is negligible compared with the total 
neutral number densities. More important is the photodissociation of molec- 
ular species, especially molecular oxygen. For many problems it is adequate 
to set this term to zero. Such considerations, however, are best left to 
particular applications. For a detailed discussion of these terms the work by 
Banks and Kockarts (1973) (hereafter denoted B&K) may be consulted. 


Finally, one of the consequences of the assumptions of (1) charge neutral- 
ity, (2) single ionization and (3) no negative ions can be obtained from the 
continuity equations. Use of equations (4-18*) and (4-19) in (4-1) and (4-3) 
gives 


9N 

— ^ + V • (N v ) = P - L 
9t e e e e 


and 


9N 

— ^ + V • (N v.) = P - L 
3t e i e e 

Current density is defined as 


(4-20) 


(4-21) 


-tv “*■ 

1 = ) q n v 

J L . n s s s 
s=e,i 

= e N (v - v ) (4-22) 

e v i e 

where -e is the electron charge and the latter result follows from the assum- 
tions. If equation (4-20) is subtracted from (4-21), the result can therefore 
be written 


V • j = 0 (4-23) 

where the vanishing of the time derivative of the space charge density, normally 
found in this equation, is a direct result of the charge-neutrality assumption. 
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For completeness the neutral continuity equation, 
(4-20) and (4-21), is written as 
8N 


n 


at 


+ v 


(N v ) — P - L 
n n n n 


corresponding to equations 


(4-24) 


4.2.2 Equations of Motion 

A. distinctive difference between the equations of motion of charged and 
neutral particles is the external force term. Electromagnetic forces, as 
noted previously, play an extremely important role in differentiating the 
motions of the charged and the neutral particles. Other forces are also 
operative in the upper atmosphere, however, and their importance must be 
assessed. These forces are the gravitational force and the coriolis force, 
the latter a consequence of working in an Earth-fixed, rotating coordinate 
system. For the neutral gas, the external force term may be written 

~ l n <F > = g + 2v x (4-25) 

p u n n n h 

n n 

where g is the gravitational acceleration, ft is the angular acceleration of 

£j 

the Earth and the brackets around F^ indicate that the average velocity, 
rather than particle velocity, is to be used. The gravitational term has a 
unique quality among all the terms in the momentum equations: for a given 

location it is constant in magnitude and direction. This makes it difficult 
to neglect in general because a situation can be conceived in which every term 
is, at some time or other, smaller than the gravitational term. So gravity 
should always be included unless it can be shown to be insignificant in a 
particular problem. No such general statement can be made of the coriolis 
term; it must be examined in the context of the other terms. 


External force terms for the charged particles are similar to equation 
(4-25), with the addition of the appropriate electromagnetic terms: 


<t > 

e 

m 


= - — (2 + — x + g + 2v x ft 

m \ c } e E 


^ y n.<l> = 

p i i 1 1 m. ' C / 

x 


+ g + 2v ± x 
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The angular velocity of a charged particle gyrating about a field line is 
given by 



(4-26) 


the magnitude of which is the gyrofrequency . Using this notation in the 
above equations gives 
-+ 

<F > ^ ^ 

— — = - ai — c - v xw + g + 2v x Q , (4-27) 

m e B e e ° e E ’ 

e 


— y n.<F.> = oj. c + v. x oj + g + 2v. x ft , (4-28) 

p . 7 X X lr> X 1 X Hi 

X X 

where the bar on the ion gyrofrequency denotes an average value obtained by 
using the mean ion mass dk. Note that in the ion equation, for example, 
the magnetic and coriolis terms can be combined to obtain 


F. + F. = v. x (oj. + 2ft ) 
xmag xcor x x E 


(4-29) 


— „ 2 -1 -5-1 

Since oj. 10 s while ft^ = 7. 29 x 10 s , the coriolis force is negli- 
x E 

gible compared to the Lorentz force. The same holds true for the electron 


equation since o) € 
neutral equation. 


10 7 s- 1 . 


So the coriolis force is retained only in the 


The momentum equations may now be written with external force terms 
written out explicitly and the notation of equations (4-14a through 4-l4c) 
incorporated in the collision terms. From equation (4-5) the electron 
momentum equation 


^ V e ->■ 1 = ^ 

—— + (v • V) v + — V • p + 03 ~ C + V X 03 - g 

3t e e N m *e e B ee° 

e e 

_ v 

= - v . ( V - v.) - v (v - v ) - rp (P - L ) . (4-30) 

ex e x en e n N e e 
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Normally, the last terra on the right side is negligible because 


(P - L )/N « 1 s' 1 , 

e e e 

while 

“ - > , n 2 -1 

v or v . ~ 10 s 

en ei 

At certain times the ratio may approach the order of one, but collision or 
other terms would still dominate. So this term will be discarded. The other 
terms require closer examination. 


To this point little has been said of the pressure tensor, p g (also called 
the momentum flux tensor), since its definition in equation (2-19). It was 
noted, following its definition in equation (2-20), that the scalar pressure 
p g is one third of the trace of the pressure tensor. For convenience, these 
definitions are rewritten in Cartesian tensor notation: 

P ag = n m<c a c g > ct.S = 1, 2, 3 (4-31a) 

p = 4 n m<c c > a = 1, 2, 3 (4-31b) 

3 a a 

where the summation convention is used for repeated indices. From these, the 
viscous-stress tensor (also called the zero-trace pressure tensor and shear- 
stress tensor) is defined by 


T ae 5 p ae ■ p 6 «e (4 ' 32: 

where 6 is the Kronecker delta. Chapman and Cowling (1970, Section 7.41), 

otp 

following rigorous procedures, obtain an expression which can be written in 
the form 


T a8 = - * 


+ 3 n 6 «6 V * V 


(4-33) 


where n is the coefficient of absolute viscosity. That derivation also pro- 
vides a rigorous means for computing n; however, for present purposes that is 
less useful than the simple mean-f ree-path expression derived from physical 
arguments . 
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These expressions can now be considered in the context of the electron 
momentum equation (4-30). The question of immediate interest is whether the 
anisotropic part of pressure tensor is negligible compared with the scalar 
(isotropic) pressure. To determine this, an order- of-magnitude estimate is 
made based on a mean-f ree-path expression for n. Expressions of Mitchner and 
Kruger (1973, Section 11-12), obtained for partially ionized gases, are 
employed. With factors of order unity ignored, the desired ratio can be 
written, for purposes of estimating order-of-magnitude only, as 


(T e ) 

ag 

„ n e 

3v 

e 

a 

P e 

‘ P e 

9x g 


(4-34a) 


The coefficient of viscosity is given approximately by 


_ _ — . m 

n - N m c l 
e e e e e 


(4-34b) 


where c e is the mean electron thermal velocity and SL^ is the electron momentum- 
randomization mean free path, which can be written 

c 


m 


e — — 

v . + v 
ei en 


(4-34c) 


When these expressions, together with 


p=NkT=Nmc 2 , 
r e e e e e e 

are introduced into equation (4-34a) , the result is 


(4-34d) 


<T e ) fi 
a6 
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3(v e ) 
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(v . + V ) 
ei en 

3x g 


(4-34e) 


A conservative estimate is obtained by taking maximum values in the numer- 
ator and minimum values in the denominator. Possible extremes would be v . + 

-1 -1 61 
v = 10 s , A(v ) - 1 km s and Ax - 10 km, very unlikely to occur in this 

en ® ot p 

combination. From this it appears that 


(T e ) 


Si 


s 10- 2 


(4-34f) 
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so that only the scalar pressure part of the pressure tensor need be retained 
in the electron momentum equation. 


Rewriting equation (4-30) to include only the retained terms gives 


3v 
e 

at 


+ (v 


V)v e + 


TZ Vp + 0 J “ c + v x m - g 

Nm r e eB e e 


= - v . (v — v . ) — v (v - v ) 
ei e i en e n 


(4-35) 


It is now convenient t d estimate the range of values each of these terms may 

be expected to have based on values listed in Table 4-1. These are given in 

1 4-15-1 

In the lastl two cases, values of 10 s and 10 s were used for 


Table 4-2 
maximum values of v 


ei 


^nd v respectively. 


Several considerations apply to the judicious use of these numbers. 

First, some terms can vanish or almost vanish, and in a number of problems 
this can be reasonably assumed. However, maximum values for these terms are 
much less likely. Second, values for some of these terms scale with the 
electron velocity, making a better comparison possible. Third, directions must 
be taken into account explicitly since the geomagnetic field makes the ionos- 
phere quite anisotropic. 


Individual terms can now be examined in light of these considerations. 
Because of the anisotropy due to the geomagnetic field, it is convenient to 
think in terms of components parallel and perpendicular to the geomagnetic 
field. This differentiation applies particularly to terms (4) and (5) , which 
are potentially the largest. Since u> e is parallel to the magnetic field 
(equation (4-26)), term (5) affects only transverse components of electron 
velocity. Because the value of this term varies directly with electron 
velocity, it may be compared with terms (7) and (8) which may vary somewhat 
similarly. Since remains virtually constant while factors in terms (7) and 
(8) can become much smaller, these latter terms can be neglected in equations 
for the transverse components. Similarly, terms (1), (2), and (6) are negli- 
gible for transverse components (radio frequency driving forces would be 
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Table 4-2. ORDER-OF-MAGNITUDE ESTIMATES OF EXTREME 
VALUES FOR TERMS IN THE ELECTRON 
MOMENTUM EQUATION 


TERM 

NUMBER 

RANGE OF VALUES 
(cm s~ 2 ) 

(1) 

-2 3 *e 5 

10 1 io 5 

(2) 

10- 5 < |(v e • V)v e l 5 10 4 

(3) 

io 4 < N \ V p <: 10 8 
N e m e K e 

(4) 

i 12 

0 < a) ■=■ c <10^ 
— e B 

(5) 

IO 9 < |v c x J e | s 10 12 

(6) 

III 1 io 3 

(7) 

0 i F e j(“ e - *j)l s 10 9 

(8) 

°i|v e „(7 e - VI * "> 10 
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required to make the first term significant). Term (3) can also be neglected, 
although it appears to be a borderline case, because there is little likeli- 
hood of the minimum value of term (5) occurring at the same location* as the 
maximum of term (3). The former corresponds to a low or middle latitude 
location and the latter to high latitudes. This reduction leaves only terms 
(4) and ^5) , which can be solved explicitly for the transverse components of 
v e , giving 



l x i 


c 


(4-36) 


This is frequently termed the E x B drift velocity or the Hall drift. Although 
equation (4-36) is expressed in Gaussian units, in the literature measured 
values of E are generally given in units of millivolts per meter. A conven- 
ient conversion for estimating drift velocities is that an electric field 
(transverse to 2) of 1 mV m ^ results in a drift velocity of about 20 m s ^ 
for B = .5 gauss. 


For the equation of the component parallel to the geomagnetic field, the 
situation is entirely different because term (5) is not part of the equation. 
It may be recalled from discussions in Section I that the existence of an 
electric field component parallel to the geomagnetic field is controversial. 
For this reason term (4) is retained in the equation for possible use, but it 
is not used to eliminate other terms. Since terms (1), (2), (7), and (8) 
scale with electron velocity, comparisons among these terms allow (1) and (2) 
to be dropped. Gravity is parallel to the magnetic field only at the poles 
where the (vertical) pressure gradient will be closer to the maximum value, so 
term (6) can be dropped. The resulting equation can again be solved explic- 
itly giving 



(4-37) 
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Equations (4-36) and (4-37) together comprise the electron momentum equations 
for the E and F regions of the ionosphere, within the few restrictions 
adopted. 


Since rather similar considerations apply to the ion equation, it can be 
treated more briefly. From equations (4-9), (4-14), (4-18), (4-19), (4-28), 
and (4-32) , the ion momentum equation can be written 


3v, 

iir + 


- v i X to i 


v ->■ 1 1 r = 

V)v. + — Vp + — V • 7 x . 

1 p i p i i 1 

— -* v m e 

~ v in (v i ' V ~ ~ v ei (v i 



(4-38) 


Corresponding to equation (4-34e) , the ratio of the anisotropic part of the 
ion pressure tensor to the isotropic part is given by 


(T±) - 
cxB 

1 

3(v ± ) 

a 

p i 

V . . + V . 

xi xn 

3x e 


(4-39a) 


As before, a maximum estimated value is obtained for (v. . + v. ) 
_1 ii in 

A(v, ) - 1 km s and Ax D - 10 km, with the result 

1 a e 


.1 s 1 , 


(T i } 


o6 


< 1 


(4-39b) 


So in the following order-of-magnitude estimates, it is sufficient to consider 
only the scalar pressure gradient - the viscous-stress tensor divergence will 
be retained or discarded along with the pressure gradient. 


Terms in equation (4-38) are evaluated from Table 4-1 as before, to 
obtain the estimates given in Table 4-3. Because the ion gyrof requency is 
much smaller than the electron gyrof requency , the electromagnetic terms do not 
dominate the ion momentum equation as completely as in the electron case. 
Comparison of terms (1), (2), (7), (8), and (9), all of which scale with v^ 
for the most part, allows all but term (7) to be dropped. Thd^oth^r terms do 
not scale together, and the ranges overlap sufficiently that \ none can be 
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Table 4-3. ORDER-OF-MAGNITUDE ESTIMATES OF EXTREME VALUES 
FOR TERMS IN THE ION MOMENTUM EQUATION 

TERM RANGE OF VALUES 
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totally ignored. Although the gravitational term is borderline, compared 
with term (6) (transverse components), it is convenient to retain it and write 
a single vector equation. The remaining tferms can be placed in the form 


->- ^ -v — ->- 1 1 p — ""j 

V. = V +— 0). — C + V, x w. + g Vp V • ) T, 

1 n V in L 11 P i 1 p i i *J 


(4-40) 


which is an implicit partial differential equation for v^ This is the final 
form for the ion equation of motion. 


Similar procedures are followed in treating the neutral momentum equa- 
tion. From equations (4-10) , (4-14) , (4-18) , (4-25) , and (4-32) this equation 
may be written 


at 


(v • V) v + — Vp + — V • 7 x - g - 2(v x ft ) 
n n p *n p L n & v 
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n E 

(P - L ) 


N m. N m 
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N m 
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N m 
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The ratio of the magnitude of the viscous stress tensor to the scalar press- 
ure, analogous to equation (4- 39a) for ions, has an order of magnitude given 
by 


(t ) 

c»e 

1 

3( V 

Ot 

P n 

’ N 

8x e 

zr- , v 


N V in + V nn 


n 


(4-41a) 


-3 -1 — 

A value of approximately 10 s is used for v at 800 km, based on a hard- 

nn 

sphere collision assumption. For neutrals, the scale size for significant 
variation is larger than for electrons and ions because the smaller scale 
sizes are associated with the geomagnetic field at high latitudes. These 
affect the neutrals only through collisions with charged particles; and because 
of the number density disparity, relatively long times are required to make 
significant changes in the neutral times. During this time, motions of the 
neutral gas cause the effects to diffuse to surrounding regions, that is, 
enlarging the spatial scale for significant variation. A conservative value 
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of 100 km is used for Ax D . Together with an extreme value of 1 km s 

p 

ACv^)^, equation (4-41a) yields 


for 


(T n> 


oB 


10 


(4- 41b) 


As before, the scalar pressure term will be used to represent all pressure 
tensor terms in the following order-of-magnitude estimates, with due recogni- 
tion of the above result in the event that the pressure gradient term appears 
negligible. 


Order-of-magnitude estimates, based on values in Table 4-1 and the above 
considerations of spatial and time scales, are presented in Table 4-4. The 
estimate for term (8) is based on a phbtodissociation rate constant of 10 s' 
for molecular oxygen (B&K, 1973); this term is clearly negligible. Also, term 
(7) may be neglected compared with term (6). 


For the neutral gas, gravity plays the prominent role in developing aniso- 
tropies in the motions, much as the magnetic field did for the charged parti- 
cles, but in a qualitatively different way. Gravity clearly separates the 
vertical and horizontal directions for different treatment. Maximum observed 
vertical neutral velocities are on the order of tens of meters per second 
(Rees, 1969; Rieger, 1974), while horizontal velocities in excess of 1 km s ^ 
have been observed (Wu , et al. , 1974). In addition, vertical scale lengths are 
of the order of tens of kilometers, while horizontal scales are at least an 
order of magnitude larger. When these considerations are taken into account, 
the vertical component equation reduces to terms (3), (4), and (6), with all 
terms except (4) contributing to the horizontal component equations. It may 
appear that term (5), the coriolis term, should be dropped. However, this is 
the only term which explicitly couples the different directions, resulting in 
a rotation of the wind direction. For sufficiently long-lived winds this 
effect is significant; hence it is retained in the neutral equation of motion. 
The final vector equation of motion for the neutral gas is then 
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Table 4-4. ORDER-OF-MAGNITUDE ESTIMATES OF EXTREME 
VALUES FOR TERMS IN THE NEUTRAL MOMENTUM 
EQUATION 
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(4-42) 
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This completes reduction of the momentum equations to forms appropriate to 
ionospheric conditions. 


4.2.3 Energy Equations 

The procedure for examining the energy equations is quite similar to that 
for the momentum equations, which is by now familiar. Only the heat flux 
vector has not been encountered in the previous equations. This was defined 
in equation (2-30) as 

+ _ 1 2 + 
q = — n m<c c> 

where c is the random (thermal) velocity. Chapman and Cowling (1970, Section 
7.4), in a rigorous derivation, show that the heat flux vector can be written 
in the form (second approximation) 

q = - K T VT , (4-43) 

T 

where K is the thermal conductivity coefficient (sometimes designated by X). 


With this result the electron energy equation may be written from equa- 
tions (4-15), (4-18), (4-32), and (4-34f) as 
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Some of these terms can be ignored forthwith by considering the relative 

magnitudes of flow and thermal energies. Maximum flow velocities are of the 

order of 10^ cm s ^ while electron thermal velocities are of the order of 10^ 

cm s \ Since energy varies as the square of the velocity, the ratio of flow 

— 4 

to thermal energy is of order 10 ; therefore, for electrons, flow energy can 

be neglected in comparison with thermal energy. Thus, all terms containing 
2 2 

V e ° r V es Can be ^ ro PP ec * f rom the right side of equation (4-44). This is why 
Joule heating (conversion of bulk motion to random (thermal)) motion is negli- 
gible for electrons. 


In the notation of Section III and the Appendix, the small ratio of flow 

2 

to thermal energy corresponds to the condition Kv q << 1. From equations in 
subsection A. 1 of the Appendix, in this limit the ratio of terms in brackets 
on the right side of equation (4-44) is given approximately by 


/ 3v \„ 2 2 r_ K 

- 3 * L 1 - 


Q(g) g 7 exp (-Kg 2 ) dg 


Q(g) g 5 exp (-Kg 2 ) dg 


] 


(4-45) 


If the momentum transfer cross section, Q(g), is assumed to vary as an integer 
power law n, within the range covered in the Appendix, then the ratio of 
integrals is given by 


7 7 n = even 

Q(g) % exp (-Kg ) dg 

Q(g) g 5 exp (-Kg 2 ) dg 

n = odd 

(Gradshteyn, I. S. and Ryzhik, I. M. , 1965, Integrals 3.461-2,3). The 
bracketed factor in equation (4-45) is then of order unity or less so that 



Terms containing the partial derivatives of collision frequencies may there- 
fore be dropped from the right side of equation (4-44). 


r + 1 


r = 0, 1, 


2r + 1 
2K 


(4-46) 


r = 1, 


With these reductions, the remaining terras of the electron energy equa- 
tion are 
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(4-48) 
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Order-of-magnitude estimates of terms in this equation are presented in Table 
4-5* The estimate for term (3) is based on a mean-free-path expression ob- 
tained by Mitchner and Kruger (1973, Section 11-12): 
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N k c l 
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(4-49) 


Term (7) is a consequence of the formalism, predicated on electrons being 
created with zero energy, thus requiring an average energy of 3/2 kT g from 
the electron gas. For now, this term will be retained; the actual physical 
situation will be discussed later in this section. For the other terms, it 
appears that terms (2) and (4) are negligible compared with the others, 
leaving 
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The ion energy equation does not have the same simplifications because 
flow energies are comparable to thermal energies, when flow velocities approach 
1 km s \ Consequently, neither the flow energy terms nor the derivatives of 
the collision frequencies can be neglected outright. From equations (4-16) , 
(4-18) , (4-19) , (4-32) , and (4-43) the ion energy equation can be written 
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Table 4-5. ORDER-OF-MAGNITUDE ESTIMATES OF EXTREME 
VALUES FOR TERMS IN THE ELECTRON ENERGY 
EQUATION 
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(4-51) 
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As with the ion momentum equation the scalar pressure is taken to represent 
the entire pressure tensor in order-of-magnitude estimates. These estimates 
are presented in Table 4-6. 


Only one term appears to be negligible under all likely conditions, term 
(5). Term (7), the generalized Joule heating term, has potentially the largest 
magnitude, depending on the magnitude of the ion-neutral relative flow veloc- 
ity. Because the relative flow velocity is negligible in many problems, it 
cannot be used alone as a basis for discarding other terms. Term (9) will be 
discussed later with the corresponding electron term. There appears to be no 
simple means of estimating the magnitude of the collision frequency deriva- 
tives short of a detailed numerical calculation. These terms will therefore 
be retained until there is sufficient reason to drop them. The reduced ion 
energy equation is then 
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The primary differences between the electron and ion energy equations were 
due to the disparity in masses. From that standpoint the neutral energy equa- 
tion is similar to the ion equation. Here the main differences are due to the 
disparity in number densities. From equations (4-17) , (4-18) , (4-32) , and 
(4-43) , the neutral energy equation can be written as 
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Table 4-6. ORDER-OF-MAGNITUDE ESTIMATES OF EXTREME 
VALUES FOR TERMS IN THE ION ENERGY 
EQUATION 
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Order-of-raagnitude estimates for terms in this equation are given in 
Table 4-7. 


From these values it appears that collisions with charged particles play 
a very small role in neutral particle energetics, with the exception of the 
Joule heating term (7). Terms (5), (6), (8), and (9) are negligible, leaving 
the following equation: 
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(4-54) 


Actually this equation is far from complete due to limitations of the formal- 
ism, specifically, omissions of inelastic collisions and interactions with the 
radiation field. These subjects are discussed in a limited way in the next 
subsection. 


4.3 INELASTIC PROCESSES IN THE IONOSPHERE 

At various points throughout preceding discussions allusions have been 
made to factors which, by assumption, have been excluded from the system under 
consideration. These include inelastic collisions, interactions with the 
(solar) radiation field and elastic collisions with high energy, nonthermal 
(non-Maxwellian velocity distribution) particles. Formal inclusion of these 
effects has been accomplished by S. T. Wu (1969, 1970) using the Chapman- 
Enskog method of successive approximation. The resulting increase in formal 
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Table 4-7. ORDER-OF-MAGNITUDE ESTIMATES OF EXTREME 
VALUES FOR TERMS IN THE NEUTRAL ENERGY 
EQUATION 
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complexity makes use of this approach somewhat formidable for numerical cal- 
culations. For applications which exercise the full extent of the formal 
framework, this added complexity may be no more than is required for a com- 
plete treatment in a less rigorous and consistent manner. But for the ionos- 
phere, a simpler approach appears more efficient. Furthermore, the method 
formulated by S. T. Wu (1969, 1970) is based on expansions of species tempera- 
tures and flow velocities in a center-of-mass system for the gas as a whole. 

It has been noted previously that such an approach is less useful for the 
ionosphere than explicit use of distinct equations for various gas*3, due to 
the relatively loose coupling observed among them. Because large differences in 
temperature and flow velocities can occur among ionospheric electrons, ions 
and neutrals, it is not clear that the expansions used by S. T. Wu (1969, 

1970) (to first order) would converge, so that truncation with the linear 
terms cpuld yield inaccurate results. For these reasons the approach adopted 
for this study has been to formulate a complete, consistent treatment based on 
elastic interactions, with nonelastic effects included phenomenologically when 
and where they are required. 

A thorough examination of these topics exceeds the scope of this investi- 
gation; such a treatment is given in the recent two-volume work by Banks and 
Kockarts (1973) . The purpose here is to discuss briefly those physical pro- 
cesses which may affect the conservation equations. In a few limited instances, 
explicit expressions are provided to allow these processes to be taken into 
account directly in the conservation equations. Generally, however, the 
treatment must be limited to qualitative discussion, because no convenient 
closed form expressions of direct applicability are available. In these cases 
recent reviews or model calculations are cited, where information required for 
numerical calculations can be obtained. As in previous discussions of the 
conservation equations, the continuity, momentum and energy equations are 
discussed in turn. 

Because the source/loss terms of the continuity equation were introduced 
into the momentum and energy equations in the process of eliminating time 
derivatives of number densities in those equations, it was expedient to 
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examine those terms preceding the order- of-magnitude estimates of the previous 
section. For present purposes, nothing need be added to that discussion of 
the processes whereby particles are added or lost to the various gases. How- 
ever, the concomitant effects on the momentum and energy equations require 
elaboration. First, momentum transferred from solar radiation to upper atmos- 
pheric gases is negligible. Since photoelectrons tend to be emitted in the 
direction of polarization of the photon (Heitler, 1954), and since solar 
radiation incident on the upper atmosphere is predominantly unpolarized, there 
is azimuthal symmetry in the emission of photoelectrons about the direction of 
the solar radiation. The result is that on the average essentially no net 
momentum is added to any of the gases as a result of photoionization. So 
newly formed ions on the average should have the average neutral velocity, 
resulting in an effective transfer of momentum from the neutral to the ion 
gas. Isotropy in the recombination process results in the reverse effect. In 
general, net production or loss involves only a small portion of the ion 
population, particularly at altitudes where ion and neutral flow velocity 
differences are largest. Even though the reasoning here differs from that 
used in the previous analysis of this term, order-of-magnitude estimates are 
the same in both cases. So the previous results, that the source/loss terms 
of the continuity equation have negligible effects in all the momentum equa- 
tions, remain unchanged. 

Different considerations apply to the energy equations. In photoioni- 
zation, since momentum must be conserved, the mass disparity between the newly 
formed electron and ion causes the photoelectron to carry away all but about 

one part in 10^ of the photon energy exceeding the ionization potential. The 

-3 

energy left to the ion (-10 eV) is small relative to thermal energies and, 
from the discussion of momentum, randomly directed. Newly formed ions there- 
fore have a temperature which may be expressed as T^ + AT, where, AT << T^. 

In the ion energy equation (4-52) , the last terms reflect the effects of 
deviations of the average properties of added and lost ions from the average 
ion gas properties. For the lost ions, data on the energy dependence of rate 
coefficients for loss processes presented by Ferguson (1974) in a recent 
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review indicate that such deviations are small. The last term in equation 
(4-52) should therefore be replaced with 


% - «r ! k < T „ + 4T - V 

r e 


(4-55) 


For altitudes below about 300 km, - T^, so that the temperature factor is 

small where the source factor P is maximal; above 300 km where temperature 

e * 

differences become large, P e / N e is quite small. In terms of previous order- 
of-magnitude estimates, however, this term cannot be dropped as negligible 
under all conditions. It should be noted that equation (4-55) applies only to 
photoionization. For collisional ionization, as in auroral activity, a 
thorough reformulation would be required. 


For the electron energy equation a different physical situation prevails. 
Photoelectrons are generally created with several tens of electron volts of 
energy, which places them outside the energy range of particles considered 
part of the system, as discussed previously. The energy they provide to the 
thermal gases is treated as if it were due to an external heat source. Thus, 
photoelectrons are placed in the incongruous roles of being counted as new 
particles in the system in the continuity equation, but being considered as an 
external heat source in the energy equation. These roles are somewhat recon- 
ciled by the fact that at sufficiently low altitudes thermalization of the 
photoelectron takes place rapidly; that is, through inelastic collisions with 
neutral particles and finally through elastic collisions with thermal electrons 
the photoelectrons energy is quickly reduced to less than a few electron 
volts. When the delay between creation and thermalization is short, few 
problems of consistency arise in the time evolutions of the conservation 
equations over the time scales of interest. 

More serious is the distance the photoelectron travels before it is 
thermalized. Photoelectrons created above ~300 km and moving upward approxi- 
mately parallel to the geomagnetic field may suffer no collisions to inhibit 
their motion until they traverse the entire field line to the conjugate iono- 
sphere. Charge neutrality assures that the local rate of photoionization 
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remains the local rate of electrpn production, because newly created ions will 
not travel far. However, photoelectron energy may be deposited nonlocally. 

In a recent review, Cicrone (1974) has concluded that photoelectron transport 
and energy loss within the ionosphere (local effects) are well understood; 
however, understanding of nonlocal effects due to photoelectron transport 
between conjugate ionospheres along geomagnetic field lines remains rudimen- 
tary. For detailed discussions of photoelectron properties with effects on 
ionospheric gases, and for problems requiring soured functions due to photo- 
ionization, the review by Cicerone (1974) and references therein can be con- 
sulted. Such source functions should replace the final term in equation (4- . 
50) which does not adequately represent the physical situation. 

In previous discussion of the neutral continuity equation it was noted 
that photodissociation is the primary source term. When order-of-magnitude 
estimates were made for this process, effect in the energy equation were 
found to be negligible. However, once again this estimate failed to consider 
nonelastic effects, in particular excess energy of the photon beyond that 
required for dissociation. This contribution will be compared with those of 
other nonelastic processes below; here it is sufficient to note that although 
the final term in equation (4-53) is dropped, the physical process behind it, 
photodissociation, does make a contribution in another form. 

The conclusion to be drawn from these discussions is that in the ionos- 
phere, the source/loss terms of the continuity equation do not contribute 
significantly to the momentum equation. They do have significant effects in 
the energy equation, although not in the form suggested by the formalism. 

This is due to the neglect of inelastic processes in the formalism, and these 
processes are the basis for effects in the energy equation. 

To this point, discussion has been centered on particle sources and 
losses and the associated effects on momentum and energy equations. Sources 
and sinks of momentum and energy themselves must now be considered. Source/loss 
terms for the momentum equations are already included in the form of the 
external force terms; these require no further discussion. However, in the 
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energy equations, as in the continuity equations, major source/loss mechanisms 
are not included within the elastic-collision formalism developed above. These 
involve energy transfer through inelastic collisions and radiative sources and 
sinks of energy. The former is important primarily for electrons and neutrals, 
the latter for neutrals alone. This separation of effects can be understood 
from the properties of the gases in the upper atmosphere. 

The importance of inelastic collisions for electron energy transfer is 
another consequence of the disparity between electron and heavy-particle 
masses. Equation (3-48) shows that energy transferred per elastic collision 
between ions and neutral particles is equivalent to the difference in mean 
particle energies, while for electrons and neutrals it is this energy differ- 
ence reduced by a factor of (m^/m^) . Consequently the elastic collision is an 
efficient means for the transfer of energy from ions to neutrals, but rather 
inefficient for electrons to neutrals. In this context the relative magnitudes 
of elastic and inelastic cross sections must be considered. Elastic cross 

sections are of the order 10 ^-10 ^ cm^ while inelastic cross sections are 

—17 —18 2 

generally of the order 10 -10 cm or smaller at energies of an electron 

volt or less. Since the most efficient inelastic collisions for energy trans- 
fer will be those with transition energies near the kinetic energies of the 
colliding particles, the relative magnitudes of these cross-sections indicate 
that energy transfer for ions through inelastic collisions is negligible 
compared with elastic collisions. However, for electrons the smaller inelastic 
cross section is more than compensated by increased energy transfer per inelastic 
collision, making this the predominant means of electron energy transfer, up 
to heights. where electron-neutral collisions become less important than electron- 
ion collisions. 

The importance of a given inelastic process for electrons depends on the 
transition energy for the process and the electron energy distribution. Con- 
servation of energy requires that the kinetic energy of the collision partners 
equal or exceed the transition energy in order for the cross section to be 
nonzero. Thus, the significant transitions are those with energies of the 
order of or less than the mean thermal electron energies — a few tenths of an 

4-33 


I 


electron volt. For the primary neutral species of the upper atmosphere, 
rotational and vibrational excitation of .and 0^ and excitation of the fine 
structure levels of the ground state of atomic oxygen meet this criterion. 

For unusually high electron temperatures (T^ > 4500°K) , excitation of elec- 
tronic states of atomic and molecular oxygen becomes significant. These and 
other inelastic collisions have been reviewed by Dalgarno (1969). The appli- 
cation of these processes to electron energetics has been examined by a number 
of investigators, notably Dalgarno, et al. (1968). To facilitate computer 
calculations, these investigators have developed analytic expressions for 
electron energy transfer rates due to inelastic collisions; these will be 
extracted from the literature and presented here' without detailed discussion 
of their derivation. 


For energy transfer due to electron collisional excitation of the rota- 
tional levels of , the expression given by Dalgarno (1969) is 


N 2 (R0T) = q L - - N e n„ 2 

x (2.9 x 10 -14 )(T - T )/T 1/2 eV cm -3 s _1 (4-56) 

e n e 

where n is the number density of N . The corresponding energy transfer rate 
2 

for (> 2 , also given by Dalgarno (1969), is 
0 2 (R0T) = Q l - - » e 

(6.9 x 1(T 14 )(T - T )/T 1/2 eV cm -3 s _1 . (4-57) 

e n e 

me energy transfer rate due to electron excitation of the vibrational states 
of N^, as presented by Stubbe and Varnum (1972), is 



(4-58) 
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where 


f - 1.06 x 10 4 + 7.51 x 10 3 tanh[l. 10 x 10~ 3 (T e - 1800)] 

a n d , 

g = 3300 + 1.233(T e - 1000) - 2.056 x lCf 4 (T e - 1000) (T & - 4000) . 


The corresponding, expression for O^, given by Prasad and Furman (1973), 


is 


°2<VIB>= \--» e V 


x(5.196 x 10 1 ^) exp 


V 700 T e ) J 


x ( 1 - exp 


[- 2770 (■ VvO) 


XT -1 

eV cm s 


(4-59) 


where 


f = 3300 - 839 sin[l. 91 x 10 4 (T - 2700)] . 

e 


3 1 

Electronic excitation of atomic oxygen (0( P) 0( D) : transition 

energy = 1.97 eV) becomes significant for high electron temperatures such as 
are found in aurora and stable mid-latitude red arcs (SAR-arc) . An expression 
for energy transferred by electrons to 0 by this method, also provided by 
Stubbe and Varnum (1972) , is 


°(D): Q L =- N e " 0 

-12 I / T e ~ 3000 ^ 

x (1.57 x 10 ) exp 1 " 


3000 T 


[■( 

1 - e«p [- 22713 (Vvll) 




,7 "3 "I 

ev cm s 


(4-60 


where 


f = 2.4 x 10 4 + .3(T - 1500) - 1.947 x 10 5 (T - 1500) (T - 4000) 

e e e 


Excitation of electronic states of molecular oxygen is less likely to be 
significant because at altitudes where molecular oxygen is more important, 
electron temperatures are not likely to be high. However, conditions for 
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these processes to be significant might be found in aurora. Therefore, the 
following expressions obtained by Prasad and Furman (1973) are presented: 


where 


(0 0 (^A ): transitional energy = .98 eV) 

^ g ^6 


W ! Q L-- N e ' , 0 2 

x 1 1 - exp [-1U00 (V/-")] ) 


f = { 13200 + 1410 sin[2. 41 x 10 4 (T - 500)] 

e 

) 

1 + exp[(T e - 14011) /1048] 


3 — i + 

(0 n ( I ) 0 ( Z ): transitional energy = 1.64 eV) 

2 g 2 g 


where 


0 (X ): Q = - N n 

z g L e CL 


x (1.616 x 10 16 ) exp 


1 - exp 


r n - 1500 \-| 
|_ 8 1500 T, /J 

[ n4oo (W)]| 


= ( 19225 + 


560 sin[3. 83 x 10 (T - 1000] 

e 


x |l + exp[(T e - 16382) /1760] 


(4-61) 


(4-62) 


Fine structure excitation of ground state atomic oxygen requires a more 
extended discussion (1) because its cross section is of the order of the 
elastic cross section, (2) because it is thought to be the dominant mechanism 
whereby electrons lose energy to neutrals, and (3) because the cross sections 
have recently been recalculated, casting some doubt on the previous assertion. 
Atomic oxygen has a T ground state with nondegenerate total angular momentum 
states, the J = 0 and J = 1 states lying .028 eV and .02 eV above the J ^ 2 


4-36 


ground state respectively. Dalgamo and Degges (1968) first suggested that 
this mechanism might be an efficient means of cooling ionospheric electrons. 
They presented electron cooling rates based on collision strengths computed by 
Breig and Lin (1966). The importance of this mechanism for electron energetics 
was confirmed by Dalgarno, et al. (1968). New calculations of cross sections 

9 

for this process were reported at a recent American Physical Society meeting 
by Tambe and Henry (1973). Cooling rates calculated from these new cross 
sections are about a factor of three smaller than for the previous ones (G. A. 
Victor, 1974). Since these recent developments have not yet appeared in the 
open literature, their implications for electron energetics are not yet clear. 
However, it appears that some caution is in order in making and interpreting 
calculations involving the electron energy equation. 


In the absence of an available electron cooling rate expression based on 
the new atomic oxygen fine; structure excitation cross sections, that based on 
the previous cross sections is provided, recognizing that it is likely to be 
superseded in the near future. A detailed derivation of the cooling rate due 
to fine structure (FS) excitation of atomic oxygen, demonstrating the book- 
keeping required for a multi-state system, is presented by Comfort (1970). 

That expression, determined from the collision strengths of Breig and Lin 
(1966), is 

N e n Q (6.217 x 10 -13 )T e l/2 

0(FS) : Q l ° " 5 + 3 exp (- .02/k T ) + exp(- .028/k T *) 

n n 

X | 6 [exp (- .02/k T e > - exp(- .02/k T n > ] 

x [1.323 - 1.028/f 2 ] + 2.8[exp(- .028/k T £ ) 

- exp(- .028/k T )] [.9911 - .9555/f 2 ] + .8 exp(- .02/k T ) 
n n 


x [exp(- .008/k T fi ) - exp(-.008/k T^) ] 

X [1.505 - 1.510/f 2 ] eV cm -3 s -1 | , (4-63) 


where f = 8k T +1 
e 
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As noted above, the cross section for excitation of the fine structure 

levels of ground state atomic oxygen is of the order of the elastic momentum 

“16 2 

transfer cross section (~10 cm ). The basis for ignoring momentum transfer 
in inelastic collision was previously stated to be that inelastic cross sections 
were much smaller than momentum transfer cross sections; the frequency of 
inelastic collisions would therefore be so much smaller than that of elastic 
collisions that the total momentum transferred through inelastic collisions 
would be negligible by comparison. Since for this inelastic process, such 
reasoning is invalid, Schunk and Walker (1970) have suggested that an effec- 
tive increase in the momentum transfer collision frequency, due to fine struc- 
ture excitation of atomic oxygen, could have a significant effect on iono- 
spheric electron transport properties. It would therefore be desirable to 
take this into account quantitatively. However, Itikawa (1972) has discussed 
the effect of inelastic processes on the momentum transfer cross section and 
concluded that the presently available formalism is inadequate to treat this 
problem. Moreover, given the uncertainties in the fine structure excitation 
cross sections noted above, it appears unprofitable to attempt a quantitative 
estimate of such effects at this time. This effect must therefore remain an 
uncertainty in any calculated results, to be taken into account in their 
interpretation. 

For ions, neither inelastic collisions nor radiative interactions play a 
significant role in the energy balance. Excitation energies for ions are too 
high for thermal particles of the upper atmosphere to interact yith them 
inelastically to any significant degree. Similarly, too few ions absorb or 
emit radiation to affect the energy balance appreciably. Interest in ion 
excited states is restricted primarily to auroral emissions, as indicators of 
interactions with nonthermal particles. 

Neutral energetics, on the other hand, depend to a high degree on radia- 
te 

tive and inelastic processes. Chandra and Sinha (1973) have recently presented 
model calculations of the neutral thermal energy budget. Results of their 
calculations for local noon show that below 300 km primary heat sources for 
neutrals are inelastic collisions with photoelectrons and photo-dissociation 
by solar ultraviolet radiation. Above 300 km inelastic interactions with the 
thermal plasma deliver energy comparable to the photoelectron contribution. 
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while photo-dissociation is a negligible heat source. Energy losses are due 
primarily to thermal conduction and, to a lesser extent, infrared emissions. 
Smaller heat sources in the calculation include chemical energy (energy 
released in dissociative recombination processes) and neutral wind energy 
(Joule or frictional heating in collisions with ions due to flow velocity 
differences). At night, of course, heat sources based on solar radiation 
vanish, resulting in a very different heat budget. Also, since the calcu- 
lations of Chandra and Sinha (1973) are based on middle latitude conditions, a 
rather different set of relationships can be expected at auroral latitudes. . 
This simply emphasizes that the importance of a given energy source depends on 
the conditions associated with each problem. Each problem, therefore, requires 
a detailed review of energy sources and sinks in order to include all that are 
of possible importance without including unnecessary computational burdens. 

The paper by Chandra and Sinha (1973) and B&K (1973) , as well as references 
therein, may be consulted for quantitative expressions and detailed discussions. 

This discussion of inelastic and radiative processes has been necessarily 
incomplete and, for the most part, qualitative. Its primary purpose has been 
to note physical considerations of importance which are not contained in the 
formalism previously developed. In many cases, these involve complicated 
problems requiring extensive computations in their own right. Since they are 
generally treated in this manner, it is more practical to incorporate results 
of such treatment into a calculation phenomenologically than to include every- 
thing into the formalism and solve it in one computational effort. This 
approach is illustrated in the problem treated in Section VI. References 
cited in the discussion above should provide reasonable resources for obtain- 
ing the detailed information required for specific calculations. 

i 

4.4 CONSERVATION EQUATIONS APPROPRIATE TO THE IONOSPHERE 

The results of this section are summarized in the final conservation 
equations which have evolved from consideration of the physical conditions in 
the ionosphere. These equations represent the starting points for the investi- 
gation of particular ionospheric problems. For convenient reference they are 
collected here. 
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4.4.1 Continuity Equations 


3N 

Plasma: -r- 5 - + 7 * (N v ) - P - L 
o t e e e e 

3N 

Neutral: -r- 11 + V • (N v ) - P - L 
ot n n n n 


4.4.2 Momentum Equations 
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4.4.3 Energy Equations 
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(4-64a) 

(4-64b) 

(4-65a) 

(4-65b) 

(4-65c) 

(4-65d) 

(4-66a) 
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(4-66c) 


n x ,n 


In the electron and neutral energy equations, the final terms represent sums 
over inelastic and/or external energy production (P) and loss (L) processes, 
as discussed in subsection 4. 3 of this section (for example equations (4-55) 
through (4-63)). 
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Section V 


TRANSPORT PROPERTIES 

5.1 MACROSCOPIC APPROACH TO TRANSPORT PROPERTIES 

In the derivation of the conservation equations (Section II) , a set of 

velocity moment equations is generated by multiplying the Boltzmann equation 

2 

by successively higher powers of velocity (terminated at v here) and inte- 
grating over velocity space. An infinite number of moment equations could be 
generated in this way. However, the second term on the left side of the 
Boltzmann equation (2-2) contains velocity as a factor, resulting in a velocity 
moment one order higher than the moment of interest. Effectively, this intro- 
duces an additional variable into the system of equations. As a consequence, 
any finite sequence of velocity moment equations generated from the Boltzmann 
equation is an indeterminate system of equations. However, if the highest order 
moment can be related to lower order moments, essentially independently of the 
system of moment equations, the system can be closed and solved, in principle. 
This is one of the mathematical motivations for studying transport properties. 

Physical motivation stems from the possibility of better understanding 
the relative importance of different mechanisms which convey particle proper- 
ties from one location to another. At the elemental level, a particle can be 
thought of as having certain average properties, characteristic of its loca- 
tion. In the context of a Boltzmann binary collision, a particle moves undis- 
turbed (except for external force fields) from one collision to another; these 
collisions determine the length of a free path. Thus, the particle carries 
the average properties characteristic of the beginning of the free path to the 
other end of the free path. For this reason, Chapman and Cowling (1970) use 
the term M free path M phenomena as an alternate to "transport" phenomena. 

Clearly, if average particle properties at both ends of a free path are the 
same, nothing is changed. So the effects of transport depend on a spatial 
gradient in the property of interest. With a spatial gradient, the total 
effect of transport then depends on the length of the free path, which in turn 
depends on velocities, number densities, and collisional interaction potentials. 
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Using mean free path arguments of this nature. Chapman and Cowling (1970) 
derive expressions for diffusion, viscosity, and thermal conductivity coeffi- 
cients, which relate to gradients in number density, flow velocity, and tempera- 
ture, respectively. These expressions are useful for physical understanding 
and order-of-magnitude estimates. However, for rigorous derivations with 
accurate results. Chapman and Cowling (1970) (and others) use the Boltzmann 
equation. 

For charged particles, the effects of electric and magnetic fields must 

be considered in addition to spatial gradients. An immediate consequence of an 

external magnetic field is that it renders the medium anisotropic. This is 

easily seen in terms of the simple free path arguments^ For a strong magnetic 

field (a), >> v) , the effective free path in the transverse direction is no 
b 

longer determined by collisions, but by the gyroradius . In the parallel direc- 
tion, of course, collisions are still the limiting factor. Furthermore, the 
Lorentz force causes the transverse component of external forces to drive 
particles in a direction transverse to both the driving force and the magnetic 
field, for example, equation (4-36). The situation is more complicated when 
the collision frequency and the gyrofrequency are comparable, for example, 
ions in the E region; both collisions and the magnetic field then affect the 
transverse free path. The mathematical consequence of these physical effects 
is that each transport coefficient must be replaced by a three dimensional 
tensor of rank two. 

Shkarofsky (1961) has developed a method for calculating electron trans- 
port tensor components, explicitly including the effects of a magnetic field 
of arbitrary strength and an arbitary degree of ionization. Use of this method 
in calculating transport tensor components is rather difficult for the general 
case. However, if electron-neutral collision frequencies vary with velocity 
according to a power law r, Shkarofsky (1961) presents results with all the 
integrals evaluated. Although these results still require lengthy calcula- 
tions, they represent a considerable reduction of effort over the general case. 

With this motivation, methods have been developed for determining an 
effective power law for electron-neutral collision frequencies of arbitary 
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velocity dependence. This makes complete results of the Shkarofsky formula- 
tion immediately accessible. Derivation of this technique is presented first. 
This is followed by a discussion of specific elastic collision frequencies for 
upper atmospheric gases . The power law approximation method is applied to the 
electron-neutral collision frequencies as an example of its use. With these 
parameters available, expressions for electron transport tensor calculations 
are extracted from the derivation by Shkarofsky (1961) • These are examined 
briefly for consistency with order-of -magnitude estimates in the previous 
section. Expressions of more limited applicability are presented for ions 
with a discussion of the consequences. Finally, available expressions for 
neutral transport coefficients are presented. 

5.2 POWER LAW REPRESENTATION FOR COLLISION FREQUENCIES 

As noted previously, calculations of transport coefficients by the method 
of Shkarofsky (1961) are considerably simplified if the collision frequency 
varies as a power law in velocity. Recognizing this, Schunk and Walker (1970) 
presented a technique for determining an approximate power law representation 
for electron-neutral collision frequencies when the velocity dependencies of 
the momentum transfer cross sections are expressed as polynomials. An alter- 
nate technique has been developed which offers the following advantages over 
the Schunk and Walker (1970) method. It gives a power law representation for 
arbitrary velocity dependence; It has a simple extension to multi-species 
gases; and it is derived straightforwardly without approximation from the 
defining equation. Most importantly, comparisons of this method and the Schunk 
and Walker (1970) method with results of exact calculations in a model study 
for the ionosphere demonstrate this method to be the more accurate (Comfort, 
1975). 

In the following derivation, it is convenient to adopt the notation of 
Shkarofsky (1961). This permits an orderly introduction of that notation with 
relationships to notation used here; in the later presentation of transport 
tensor components, use of this notation will facilitate reference to the source 
work . 
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From equation (A-3) in the Appendix, it is seen that beyond the first term 

in the expansion of the momentum transfer collision frequency, all terms are 

of order (Kv q ) or higher. For electron-heavy particle collisions K - m e /2kT e , 

which is approximately the reciprocal of the electron mean square thermal speed 

Since v q is the magnitude of the relative flow velocity [ | (j = i, n), 

7k V is the ratio of relative flow velocity to mean thermal velocity, which is 
V ° -2 

of order 10 or less in the ionosphere. Thus, electron collision frequencies 
may always be computed in the limit v q -* 0, given by equation (3-36) . This is 
consistent with the expressions of Banks (1966a), Shkarofsky (1961), and 
Itikawa (1971, 1973). 


When the isotropic part of the electron velocity distribution is assumed 
to be Maxwellian, the total effective electron-neutral collision frequency for 
momentum transfer is given by Shkarofsky (1961) as 

00 

<v > = — — v (v) exp(-Kv^)dv, (5-1) 

m 3f J m 

* o 

where K is as above, v is the total electron velocity (approximately the total 
relative velocity g of Section III) , and (v) is the velocity-dependent colli- 
sion frequency. For a single neutral species, v^(v) is related to the momentum 
transfer cross section by 


v m (v) = N n v ^(v) 


(5-2) 


With this definition, equation (5-1) is seen to be equivalent to equation (3-36) 

hence <v > = v . If v (v) varies with v as a power law of the form 
m en m r 


v (v) - c v 
m 


(5-3) 


where c is proportionality constant and r need not be an integer, the integra- 
tion in equation (5-1) can be carried out immediately, yielding 


<v > = 
m 


K-r /2 r(i-P)/r(f> 


(5-4) 
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For given by equation (5-3), Shkarofsky (1961) gives the following 
relation (his equation (99)): 


(W " V 2 < r + 5) , 

m a m a 


(5-5) 


where W = Kv and, for arbitrary function (j> , the average ( ) is defined by 


(<t») a = Air 


(f f I 


exp(-Kv 2 )v 2 dv . 


(5-6) 


Equation (5-5) can be verified by straightforward integration; however, an 
alternate method indicates a convenient method for determining the effective 
power index r. From equations (5-1) and (5-6), the denominator and numerator 
of the left side of (5-5) are given by 


(Wv ) = 3/2 <v > 

m a m 


(5-7) 


and 


(W 2 v ) - 4 =. [ v v 6 exp(-Kv 2 )dv . 

m a yjw J m 


(5-8) 


Let I(K) denote the integral in equation (5-1); thus from that equation. 


I(K) 


3yfT -5/2 

-S- K <v_> • 


(5-9) 


Then the integral in equation (5-8) may be written as 

oo 

j V m v6 ex P<- Kv2 > dv - - Ik i(k > " - Ik (k " 5/2 <V m >) * (5_10) 


in which the order of integration and differentiation has been reversed. 
Equation (5-8) may therefore be rewritten 


* - I 1(7/2 « <K ' 5/2 'V> • 


(5-11) 
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If v varies as a power law in v, as in equation (5-3), equation (5-4) 
m 

can bemused in (5-11) and the differentiation carried out. The, result, to- 
gether with equation (5-7), verifies equation (5-5). It should be emphasized 
that the steps leading from equation (5-8) to equation (5-11) depend on a 
Maxwellian velocity distribution for electrons, but place no particular 
restriction on the velocity dependence of v^. 


In the following, it is assumed that <v > is a known function of K (that 

m 

is, T e )> based on a knowledge of v m and equation (5-1). Equation (5-5) is now 

taken to define an effective value of r for those cases in which v does not 

m 

have the simple form of equation (5-3) . Solving equation (5-5) for r yields 



Substitution of equations (5-7) and (5-11) into (5-12) and performing part of 
the differentiation results in 


r 


2K 3 < > 

<v > 3K V m 
m 


(5-13) 


which provides a convenient means for evaluating r, once <v > is known as a 

m 

function of T^. It is noted that result (5-13) follows from the defining 

equation (5-5) without approximation; and the functional form of <v > is 

m 

unrestricted . 


Although equation (5-13) is sufficiently flexible to treat multi -component 
gases, a more convenient form for such gas mixtures is easily obtained. Let j 
designate the j ^ neutral species. For each neutral species an effective power 
law index r_. is computed from equation (5-13) . The total electron-neutral 
collision frequency is just the sum of the individual species collision fre- 
quencies 

<V = l <Vj> = l Cj K j/ r (-^ "2 -)/r(5/2) , (5-14) 
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where each <v j > is represented by equation (5-4), using the computed value of 
r . . When equation (5-14) Is used to evaluate the derivative in equation (5-13) , 

J 

the result may be written as 


y r 

< v i j 




(5-15) 


This specifies r simply as the mean of the power law indexes weighted by the 
corresponding species effective collision frequencies. 


A brief application of these results to ionospheric collision frequencies 
is given in the following subsection. Accuracy of this representation is dis- 
cussed in terms of ionospheric transport tensor component calculations by 
Comfort (1975) . It is found that for the collision frequencies presented in 
the next section, use of the above approximation results in deviations from 
exact calculations of less than .2 percent, in a model calculation. 


5.3 EFFECTIVE MOMENTUM TRANSFER COLLISION FREQUENCIES FOR THE IONOSPHERE 

Formal expressions for collision frequencies have been derived in Section 
III and the Integrals for particular velocity dependencies evaluated in the 
Appendix. Here explicit expressions for the different species of the ionosphere 
are presented in a form suitable for computer evaluation. Some of these are 
obtained directly from the literature; others are introduced here based on 
momentum transfer cross sections or other information in the literature. 


5.3.1 Elsctron-Nsutral Collisions 

The most recent reviews of momentum transfer cross-section information 
(experimental data and theoretical calculations) are those of Itlkawa (1971, 
1973), based on results available at the end of 1971. Itikawa (1971, 1973) 
defines two effective collision frequencies, one for computing a. c. conduct- 
ivities, the other for d. c. conductivities. That for d. c. conductivities is 
designated <v eff > and is defined as in equation (5-1) . Itikawa (1971) shows 
that <v eff > is appropriate for use in the transport calculations of interest 
here; accordingly those collision frequencies are presented here. 
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Collision frequencies are given by Itikawa (1973) as functions of in 
tabular form. If, for each neutral species, is written in the form 


m 


'J ■ "3 Jo 


(5-16) 


then equation (5-4) can be used to obtain 

vk/2 


<v > = 


-i m 

1 y 

/ 9 ^ L 


j T(5/2) k ^ 0 °jk 


m ■ m 


(5-17) 


where K = has been used. A fourth degree polynomial in 


1/2 


has been 

fit to the collision frequency data of Itikawa (1973), Table II). Equation 

(5-17) has then been used to determine the a., coefficients from the curve fit 

J & 

coefficients; these a., coefficients are listed in Table 5-1. The resulting 

J k 

polynomial approximations represent Itikawa 1 s (1973) tabulated collision fre- 
quencies within 2 percent over the range 200°K <_ T^ _< 5000°K. 


Table 5-1. COEFFICIENTS FOR ITIKAWA'S (1973) ELECTRON-NEUTRAL 
COLLISION FREQUENCIES (CGS UNITS) 


j 

a jl 

a j2 

a j3 

a j4 

N 2 

-2.32 x 10 -16 

6.60 x 10 -23 

-1.28 x 10" 30 

9.12 x 10 

°2 

1.50 x 10" 16 

3.87 x 10 -25 

3.24 x 10" 31 

-2.86 x 10" 39 

0 

5.38 x 10 -17 

4.32 x 10 -24 

5.10 x 10‘ 32 

-4.85 x 10" 40 

H e 

4.90 x 10" 16 

5.84 x 10 -24 

-4.69 x 10 -32 

8.02 x 10" 41 

H 

4.23 x 10' 15 

7.07 x 10' 24 

-9.76 x 10" 31 

6.50 x 10' 39 


As an application of the power law representation for collision frequencies, 
derived in the previous section, equation (5-17) can be substituted into equa- 
tion (5-13) and the differentiation carried out to obtain 


r 


j 


^ 

<Vj> r( 5 / 2 ) 



(5-18) 
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which is readily evaluated with the coefficients in Table 5-1. 


Numerical 

results are given by Comfort (1975). Equations (5-17), (5-18), and (5-15) are 
particularly convenient forms for computer evaluation. 


5.3.2 Electron— Ion Collisions 

Collisions among charged particles in plasmas have received extensive 
attention in the literature and in many texts, so that a detailed discussion 
need not be presented here. It is nevertheless worth noting that a variety 
of methods have been used, each focusing on a particular physical aspect of 
the process. To first order, all rely on either truncating the Coulomb poten- 
tial at the Debye length (R^) or using the screened Coulomb (or Debye) poten- 
tial (* exp (-r/R^)/r) with or without truncation at R^. The physical idea 
behind these approximations is that charges of a given sign are effectively 
screened by charges of the opposite sign, so that beyond a certain distance 
(Rp) there is negligible interaction. Among the approaches used to compute 
the effects of charged particle collisions are the binary collision concept, 
the Fokker-Planck equation, and two-particle distribution functions. These 
are examined in detail by Shkarofsky, et al. (1966, Chapter 7) and Burgers 
(1969, Chapter 4). 


To first order, the various approaches arrive at results which can be 
placed in the same form, with differences incorporated into the Coulomb 
logarithm (Jin A). Numerically, ionospheric collision -frequencies resulting 
from these methods agree within 10-15 percent. Since there appears to be 
little consensus favoring one result over others, and given the reasonable 
agreement among them relative to other uncertainties in ionospheric calcula- 
tions, the basis for selection here is simplicity and consistency with other 
factors in the study . 

With that rationale, the derivation of the electron- ion collision fre- 
quency is outlined in the context of the binary collision formalism developed 
in Section III. It is well-known that both classical and quantum mechanical 
treatments yield the Rutherford differential cross section for scattering in 
a Coulomb potential: 
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9 


(5-19a) 


a(g,X) 



sin (x/2) 


in the notation of Section III, with Z and Z the number of charge units for 

r s 

species r and s. The momentum transfer cross section for this process is then 
given by equation (3-12) as 


Q rs (g) " 2ir 


fa?) ! 


(l-cos x) 


sin 4 (x/2) 


(5-19b) 


Here the relation between the impact parameter and the differential cross 
section. 


bdb = -a (g,x) sinxdx. (5-19c) 

has been used. 


Since the integral in equation (5-19b) diverges at the lower limit, which 
corresponds to forward scattering and infinite impact parameter, the standard 
technique for avoiding this problem is to set a new lower limit x m » This is 
determined by assuming that because of screening by charges of the opposite 
sign, the potential of a given particle does not effectively extend beyond the 
Debye length 


*D 



(5-19d) 


With this assumption and the condition that densities be sufficiently high that 
the Debye sphere contain many charged particles, the integration in equation 
(5-19b) can be carried out and the result placed in the form 


Q rs (g) “ 16 it 


(Z Z e 2 \ 4 [y g 2 “| 

Ur.* 2 / 


Jin 


(5-19e) 



At this point, the usual approximation is to argue that the logarithmic 

factor varies sufficiently slowly, with velocity that the mean velocity can be 

used in its argument and the logarithm can be factored out of the integral in 

the collision frequency. For the electron-ion collision frequency, 

p *m and equation (3-36) is used for the integration over relative velocities, 
rs e 

giving the result 


ei 


= < v .> 
ei 


4 

3 


N e 
e 


* .\TT ^ 3/2 *nA, 


(5-20a) 


where the argument of the Coulomb logarithm (£n A) is given by 


A 




(5- 20b) 


This form for the collision frequency is consistent with those of Banks 
(1966a), Shkarofsky, et al. (1966) and Itikawa (1971) among others. The form 
for the Coulomb logarithm is one of the most common in the literature. How- 
ever, it neglects the effects of ion motion both in charge shielding and in 
contributing to the relative collision speed. Since the ion thermal speed is 
about one percent of the electron thermal speed, these effects are relatively 
small; but they contribute to the numerical differences discussed above. For 
present purposes, equations (5-20a) and (5-20b) are considered adequate. 

5.3.3 Ion— Neutral Collisions 

Ion-neutral collisions have special significance for several reasons. First, 
elastic collisions are an efficient means for transferring momentum and energy 
between ions and neutrals, so that the formalism developed in Section III 
applies directly. Second, it was shown above that all collisions involving 
electrons could be computed in the limit of vanishing relative flow velocity. 
However, this is not true for ions in general. Hence, the full mathematical 
framework developed in Section III, taking into account the effects of the 
relative flow velocity on momentum and energy transfer, must be used for ion- 
neutral collisions and only for those collisions. Third, because the masses 
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of ions and neutrals are comparable, electrostatic fields of external (mag- 
netospheric) origin give momentum to the charged particles; and this is 
transferred to neutrals primarily through collisions with ions. In other 
words, the ions are the agent through which maghetospheric electric fields 
interact with the neutral atmosphere. Some effects of this interaction are 
examined in the following section. 

Since there are five major neutral species and five major ion species, 
a total of 25 different collision frequencies are possible. However, all of 
these species are not important in the same altitude region. Thus, some 
possible combinations are not required. Those that are required may be 
conveniently grouped into two categories: collisions between ions and unlike 

neutrals, and resonance charge exchange collisions between ions and parent 
neutrals. Interactions in each category are the same, differing only in 
strength, so that only two different expressions are required to express all 
the collision frequencies. These are discussed separately below. 

5. 3. 3.1 Collisions Between Ions and Unlike Neutrals . The longest range inter- 
action between ions and neutrals results from an induced dipole potential of 
the form 


(p = -ae 2 /2r\ (5-21) 

where a is the dipole polarizability of the neutral atom. This interaction 
predominates at low temperatures, below about 300°K. Due to lack of experi- 
mental data for higher temperatures, little is known reliably of the actual 
shorter range interactions (Mason, 1970). Depending on the nature of these 
shorter range interactions, the momentum transfer cross section can either 
increase or decrease with increasing temperature. So in the absence of 
reliable information, the simple approach is taken: 'the polarization inter- 
action is assumed to predominate for all ionospheric temperatures. Corre- 
sponding to the induced dipole interaction, equation (5-21), Is the momentum 
transfer cross section (see Dalgarno, et al. (1958)) 
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, . 2.21TT /oe 2 \ 1/2 2 

Q rs <8) -— fc 


(5-22) 


where r designates an ion species and s a neutral species. This is recognized 
as the momentum transfer cross section for Maxwellian molecules, for which the 
collision frequency is independent of temperature. From equation (A-43d) the 
collision frequency for this cross section is 

/_ 2 \ 1/2 


2.217T 


(5-23) 


When values for the dipole polarizability given by Banks and Kockarts 
(1973, Chapter 9) are used in equation (5-23), collision frequencies may be 
written in the form 


v *An, 
rs rs s 


(5-24) 


where the A coefficients are given in Table 5-2, 
rs ° 


Table 5-2. COEFFICIENTS FOR ION COLLISIONS WITH UNLIKE NEUTRALS FROM 



DIPOLE POLARIZABILITIES 

OF BANKS AND KOCKARTS 

(1973) 

ION-NEUTRAL 

A r$ (lO-Ws- 1 ) 

ION-NEUTRAL A fS 

(10' 9 cm 3 s _1 ) 

o + - n 2 

.89 

0 + - H 

2.18 

o 

1 

+ o~ 

.70 

H e- N 2 

1 .84 

NO - N 2 

.90 

h:-o 

1 .29 

N0 + - 0 2 

.83 

h! - H 
e 

2.36 

N0 + - 0 

.71 

H + - N 2 

3.50 

0 + - N 2 

1.08 

H + - 0 2 

3.32 

<\j 

O 

+ 

O 

1.00 

H + - 0 

2.37 

a; 

□c 

i 

+ 

o 

.66 

H + - H e 

1 .33 
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The order corresponds to the altitude region of primary importance , with those 
listed first (left side) most important in the E region. Since H* and H + ions 
are negligible at altitudes of interest here, those collision frequencies are 
included only for completeness. 

A significant feature of these collision frequencies is that they are 
constant. This is associated with the type of interaction assumed (that of 
Maxwellian molecules) rather than any approximation involving dropping small 
terms. Since the derivative with respect to relative flow velocity vanishes 
for this type of collision frequency, that term will be absent from the ion 
energy transfer equation (4-66b) for these collisions. 

5. 3. 3. 2 Collisions Between Ions and Parent Neutrals . For collisions between 
ions and parent neutrals, resonant charge exchange interactions become more 
important than induced dipole polarization interactions as temperature 
increases (Mason, 1970). Such collisions are quite important for the exchange 
of momentum and energy between ions and neutrals. In a single collision, a 
fast ion becomes a fast neutral and a slow neutral becomes a slow ion. Because 
both momentum And kinetic energy of the particles are conserved in the collision, 
the collision may be treated as if it were elastic, by the formalism developed 
in Section III. Indeed, the small-angle charge exchange collision is indistin- 
guishable from a large angle elastic collision. Therefore, these collisions 
are discussed here in the context of elastic collisions and in that mathematical 
framework. 

Banks (1966c) has examined the transition from induced dipole polarization 
interactions at low temperatures to charge exchange interactions at higher tem- 
peratures. He finds that the transition takes place over a temperature range 
of a few hundred degrees. Outside this range, he determines that the individ- 
ual interactions dominate almost totally in their respective temperature 
regions. Banks (1966c) therefore takes the approach of extrapolating experi- 
mental or theoretical curves of cross section versus temperature for each 
interaction into the temperature region of the other. Where they cross defines 
the transition temperature and expressions for each interaction are used alone 
in their respective regions of dominance. 
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Stubbe (1968) takes an alternate approach. In a semi-classical calcu- 
lation, he corrects the contribution of charge exchange to the momentum transfer 
cross section by including contributions of induced Dipole polarization for 
impact parameters beyond a critical value. Solving the resulting equations 
numerically, he presents graphic results for the momentum transfer cross sec- 
tion of (0 + , 0) and (H + , H) for g £2.5 x lO^cm s \ as well as collision fre- 
quencies as functions of relative flow velocity, with T as a parameter, for 
|v. “ v n l £ 3 x 10^ cm s **. He also approximates the collision frequencies 

at zero relative velocities as power laws in temperature for other charge 
exchange collisions of aeronomical interest. 

In comparing his results for (0 + , 0) with those of Banks (1966c) , Stubbe 
(1968) finds his results greater than Banks' (1966c) ranging from about 35 
percent at 500°K (transition region) to 22 percent at 2000°K. Stubbe (1968) 
appears to interpret these differences at high temperatures as due to dipole 
polarization effects, ignored by Banks (1966c) at high temperatures. However, 
since they used slightly different charge exchange cross sections, exact 
agreement should not be expected. Stubbe* s (1968) own data indicate that at 
2000°K the polarization effect is less than 9 percent. In view of uncertain- 
ties in the source data, such differences are not necessarily the limiting 
factors in the accuracy of the results. 


For present purposes, it is not so much the assessment of accuracy as 
the form of the results that leads to the use of Stubbe's (1968) results for 
the present study. The desirable feature is a single representation of the 
velocity dependence of the momentum transfer cross section over the entire 
velocity range of ionospheric interest. Numerical results of this type are 
presented graphically by Stubbe (1968) in his Figures 6 and 7. Since these 
momentum transfer cross sections clearly vary inversely with relative velocity 
(g) , second order polynomials in (1/g) have been fit to these curves over the 
velocity ranges presented by means of least squares curve fits. Momentum 
transfer cross sections are then represented analytically in the form 


V 8) = V 




(5-25) 
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where the A ^ are curve-fit coefficients of appropriate units. Values of these 
coefficients are given in Table 5-3. These coefficients used in equation (5-25) 
reproduce the graphic values within 2 percent over the entire velocity range 
of the graphs. It is assumed that this expression can be used in collision 
frequency integrations for all velocities . 

Table 5-3. COEFFICIENTS FOR STUBBE 1 S (1968) MOMENTUM 
TRANSFER CROSS SECTIONS (CGS UNITS) 


COEFFICIENT 

0 + - 0 

H + - H 

A J0 

.67067 x 10' 14 

1.4056 x 10' 14 

A J1 

.33298 x 10 -9 

1.1572 x 10' 9 

A J2 

.11465 x 10' 4 

.64445 x 10' 4 


With the momentum transfer cross section represented by equation (5-25) , 
collision frequencies can be obtained directly from equations (A-43 c-e) . 
Thus, the collision frequencies have the form 


v. (v ) = n. 
in o j 


exp (-KV ) 
A + — 


k±i *■)•.•* a 


[(■■ 


+erf (TIT v o ) | K + 1 - TTT) VT + (2Kv o _1) ^3 

4Kv / o 2Kv 

0 / o 


i 


(5-26) 


where v = |v. - v L K = 4kT/m., T = 1/2 (T, + T ) and the A. f s are given in 
o 1 i n 1 j in j r • 

Table 5-2. 


When equation (5-26) is evaluated for (0 , 0) collisions, using various 
values for v q and T, comparison with corresponding graphic data presented by 
Stubbe (1968, Figure 9) shows exact agreement, within the precision of the 
graph. However, similar comparison for (H + , H) collisions (Stubbe, 1968, 
Figure 10) shows a temperature-dependent discrepancy. Since both calculations 
are made by the same computer program, changing only the A coefficients, 
the agreement for (0 , 0) validates the equation and the program. Similarly, 
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the same computer program computed and checked both sets of curve fit coeffi- 
cients; so the (0 + , 0) results validate that procedure. Review of the input 
data reveals no error there. It appears that for (H , H) collisions, the 
momentum transfer cross sections presented by Stubbe (1968, Figure 7) are in- 
consistent with his collision frequencies (his Figure 10). Collision frequen- 
cies computed from the cross sections by equation (5-26) are consistently 
larger than Stubbe ? s (1968) graphic values by amounts ranging from about 6 per- 
cent at 500°K to 18 percent at 3000°K (v =0). For present purposes, this 

+ ° 

is no problem since H is generally a minor ion below 800 km. In addition, 
the size of the discrepancy is within the uncertainties associated with ion 
neutral collision frequencies. 

In addition to the collision frequencies above, Stubbe (1968) provides 
the following approximate expressions, for 500°K <_ T _< 3000°K and v q = 0 : 

v(C> 2 > 0 2 ) = 1.17 x 10~ 9 (T/1000) * 28 n Q s -1 (5-27a) 

v(N 2 , N 2 ) = 2.11 x 10 -9 (T/1000) ’ 38 s" 1 (5-2 7b) 

v(N + , N) = 1.75 x 10“ 9 (T/1000)* 34 n N s' 1 (5-27c) 

v (He , He) = 2.92 x 10 -9 (T/1000)* 37 n Re s -1 . (5-27d) 

Of these only v(0*, 0 9 ) is likely to be significant below 800 km. For most of 
this region, v(0 , 0) is the single most important ion-neutral collision fre- 
quency because 0 is the major ion species throughout the F region (above abopt 
180 km, see Figure 1-lc) and resonant charge exchange is the most important 
collision process. 
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5.4 METHODS FOR COMPUTING TRANSPORT TENSOR COMPONENTS 

Rigorous calculations of transport properties are inherently difficult, 
requiring numerical solution of coupled sets of integro-dif ferential equations. 
As the number of independent parameters with which the transport properties 
vary increases, it even becomes difficult to place the numerical results in a 
form usable in practical calculations. Thus, when arbitrary degree of ioniza- 
tion and magnetic field strength must be considered in addition to number 
density and temperature as independent parameters, the problem becomes 
formidable indeed. This is why the method of Shkarofsky (1961) is so useful: 
it includes all these effects in a basically simple form. For the power law 
representation derived in Section 5.2., the more complicated mathematics has 
been worked out by Shkarofsky (1961) , who presents expressions requiring only 
the evaluation of determinants. Recent experimental data supports Shkarofsky ! s 
(1961) treatment of partial ionization (Albares, 1973). The primary disadvan- 
tage is that only electron motions are treated. 

Shkarofsky, et al., (1966) treat the ion pressure tensor in a manner 
similar to that for electrons. However, only a fully ionized plasma is con- 
sidered (ion-neutral collisions are neglected) and only transport associated 
with gradients in flow velocity (i.e. viscosity/pressure tensor) are treated. 
For the ionosphere, this set of limitations makes the results of little practi- 
cal value. At high altitudes where the fully ionized gas approximation is most 
nearly approached, the strong magnetic field limit (<v^> <<(J °^) applies to the 
ions, so that ion transport is restricted to the parallel direction. In this 
case, less complicated results for no magnetic field can be employed. Also, 
gradients in flow velocity are likely to be much smaller at high altitudes 
than in the E region (e.g. see Section VI) where the magnetic field competes 
with neutral collisions for dominance of ion motions. The fact of this 
competition, however, means that the fully ionized approximation is inappli- 
cable. The net result of these considerations is that for ions Shkarofsky T s 
method offers too little for the effort required to be useful for ionospheric 
applications . 
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Alternate expressions for ion transport coefficients are therefore adopted. 
Although these lack the thoroughness of the electron treatment, the effects on 
practical calculations are probably small, as discussed below. For neutrals 
different complications arise, due to the multiplicity of neutral species. 
However, adequate expressions are available from the literature, and these 
are presented. 

5.4.1 Computation of Electron Transport Tensor Components 

The Shkarofsky method for computing electron transport coefficients is 
described in detail in several works (Shkarofsky, 1961; Shkarofsky, et al., 

1966; Hochstim and Massel, 1969; Mitchner and Kruger, 1973); a lengthy 
description is inappropriate here. Since calculation of the transport 
tensor components is the primary objective here, the assumptions and approach 
used are summarized briefly. Then those equations necessary for making the 
calculations are extracted from the sources above. 

Assumptions made in the derivation of the Shkarofsky method are listed by 
Shkarofsky (1961, Section 2) and, in a slightly different form, by Shkarofsky 
et al., (1966, Chapter 4). These have been placed in the context of the iono- 
sphere by Schunk and Walker (1970), who also demonstrate that ionospheric condi- 
tions are consistent with the assumptions (where appropriate) . In the form 
given by Schunk and Walker (1970) the assumptions and conditions are as follows 

1. Only electron current and heat flow are considered; contributions due 
to ion motions are ignored. 

2. Only elastic collisions are considered (see discussion of fine struc- 
ture excitation of atomic oxygen in subsection 4.3 on the validity of 
this assumption) . 

3. A steady state is assumed - for the ionosphere this is valid because 
time scales for changes in macroscopic properties are long relative 
to electron relaxation times. 

4. The angular dependence of the electron distribution function is 
expanded in Legendre polynomials, truncated after the second term. 
Schunk and Walker (1970) note that this requires the electron flow 
velocity be much less than the mean thermal speed and that the 
electron mean free path be much less than the electron scale height; 
they demonstrate that both hold in the ionosphere. 
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5. The collision frequency is a function of velocity only, not of 
position; i.e. scatterers are uniformly distributed. This is valid 
if electron mean free paths are much smaller than both electron and 
neutral scale heights. 

6. Driving forces (electric fields and gradients) are sufficiently small 
that the equilibrium Maxwellian distribution is not disturbed; i.e. 
the isotropic part of the distrution function is assumed to be 
Maxwellian. This holds if the flow velocity is much smaller than the 
mean thermal speed. 

7. Applied electric and magnetic fields are uniform in space. 

8. The perturbation to the equilibrium Maxwellian distribution is expanded 
in a series of generalized Laguerre polynomials, truncated after four 
terms. (Hochstim and Massel (1969) carry out the expansion to 20 
terms in order to investigate convergence properties. Examination of 
their results indicates that for the ionosphere four terms gives 
sufficient accuracy.) 


These items fairly well indicate the approach as well as the assumptions 
and conditions for the calculation. The electron distribution function is 
expanded first in Legendre polynomials in the form 

f = f +1, (5-28) 

o 1 V 

in the notation of subsection 5.2. This form is then used in the Boltzmann 
equation, together with the previous assumptions. The resulting equations 
relate to f Q . From these relations general expressions for the transport 
tensors are obtained. To solve the equations, f^ is expanded in terms of 
generalized Laguerre polynomials (four terms) . Once is known the expres- 
sions for the transport tensors can be evaluated. The end products of this 
rather lengthy procedure are presented below. 

Since the required equations are taken from Shkarofsky (1961), this work 
is denoted by SI for ease of reference in the remainder of this section. Inso- 
far as feasible, the notation of SI is followed to facilitate reference to the 
complete derivation in that work. The treatment in SI includes a.c. electric 
fields; these effects are suppressed in the results below since they are not 
included in this study. 


5-20 



Consistent with the assumption above, the electron current is related to 
the- electric field and gradients in T and N g by (SI, equation (70)) 

1 - 7 . f - eD . VN - T . VT . (5-29) 

e e e 

(An equation of this form can be obtained from the electron momentum equation 
(4-35) for the steady-state case.) Here a is 'the d.c. electrical conductivity 
tensor, D is the diffusion current tensor, and T is the current flow tensor 
due to thermal gradients at constant density (SI uses t instead of T, but here 
x is used to denote the viscous stress tensor). Correspondingly, the total 
energy flow H due to electrons is related to the electric field and gradients 
in T e and by (SI, equation (72a)) 

H => 7 • £ - eQ • VN e - K • VT g . (5-30) 

Here y is the energy flow tensor due to electric fields, Q is the energy 
diffusion tensor and K is the thermal conductivity tensor at constant density. 
These transport tensors are not independent, being related by (SI, equation 
(76)) 


D = 


kT o 
e 2 N 



(5-31) 


known as Einstein relationships. Calculations are therefore limited to a, T, 
y, and K. 

In Shkarofsky's formulation, the transport tensor components are expressed 
in a compact notation with a deceptively simple form. Only in the limiting 
cases does simplicity remain as the parameters are expanded for evaluation. In 
order to present this expansion in a coherent manner, the approach taken here 
is to begin with the final result required - the transport tensor. These are 
related to other parameters, which are in turn related to still other parameters, 
until all parameters can be related to the fundamental properties of the system: 
electron temperature, electron number density, collision frequencies, magnetic 
field strength and physical constants. 
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All of the transport tensors have the form 





(5-32) 


where the 3-component lies in the direction of the magnetic field. The rela- 
tively simple form is due to the fact that the only anisotropy of the medium 
is due to the magnetic field; off-diagonal terms are due to the Lorentz force. 
Because of the symmetries in the matrix elements, they can be compactly 
expressed for the general case in the following manner (SI, equations (82)): 


N e 
e 


a ll + ia 21 M (<v >g + 'la. h ) 
e g a— b a 


T,, + IT 


N ek 
e 


11 iX 21 M e (<v >^1^) 


Pi, + iP 


5N ek T 
e e 


11 m 21 2M (<v >g +iui,h ) 

e g p— d p 


K,, + iK_, = 


5N k T 
e e 


11 21 M e «v >g K ±i^V 


(5-33a) 


(5-33b) 


(5-33c) 


(5-33d) 


In these equations <v^> the total effective electron collision fre- 
quency, given by (see notation in subsection 5.2.) 

<v > = <v > + <v .> . (5-34) 

g m ei 

The g and h coefficients are correction factors which incorporate the effects 
of the degree of ionization and magnetic field strength into the tensor compo- 
nents. The complex notation is a device for compact expression; it has no 
physical significance. To obtain the components, simply set ui^ - to = 0 
in equations (5-33); these correspond to the zero magnetic field case. 
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In equations (5-33) the only unknown quantities, in the context of the 
fundamental properties noted above, are the g and h correction factors; so it 
remains only to determine them. They are computed from (SI, equations 
(92 a-h)): 



(5-35a) 


(5-35b) 



(5-35c) 

(5-35d) 

(5-35e) 

(5-35f) 

(5-35g) 

(5-35h) 


(The factor 2 in these last two equations was erroneously omitted from equations 
(92g) and (92h) in SI, as determined from SI equations (82) and (91) and veri- 
fied by comparisons of numerical results.) Signs + in the h factors correspond 
to + in equations (5-33) . 
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In equations (5-35) the notation I^qqI denotes the 00 minor of the deter- 
minant | A | , with similar meanings for the other determinants. Matrix A is 

j 

given by (SI, equations (55c), (63a)) 


A - 


<v .> 
ei 


00 


± iw b 

H 01 

CM 

O 

« 

0 

[ 01 

H 11 - 4 ^ 

H 12 

0 

[ 02 

H 12 

"22 ± “h 

0 

I 

0 

0 

H 33± : 


L 16 S 

(5-36) 


The matrix (H_) is an interaction matrix, with contributions due to collisions 
of electrons with ions, electrons, and neutrals: 


H. 


— <v 4 > H < v . > H, . « ^ “44 

e-f xj ei f ij m ij 


hJ® + <vj> H ( 


,en 


(5-37) 


LJ., . - 'V j ' 

ij ei 

(In SI the collision frequencies in this equation are included in the defini- 
tions of H?'!', H??, and H 4 ^.) 

ij ij ij 


For the first two matrices in equation (5-37), the interaction is known 
and fixed; the matrix elements are therefore constants. These are given by 


<v - 


(53), 

(63b, c)) 



1 

3/2 

15/8 

35/16 

3/2 

13/4 

69/16 

165/32 

15/8 

69/16 

433/64 

1077/128 

35/16 

165/32 

1077/128 

2957/256 


<5-38) 


and 


<v 


ee 


0 

1 

3/4 

15/32 \ 

(5-39) 

1° 

3/4 

45/16 

309/128 1 * 

\o 

15/32 

309/128 

5657/1024/ 
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The matrix (H^) en contains effects of electron-neutral interactions; 
its elements vary with the velocity dependence of the momentum transfer cross 
sections for electron-neutral collisions. In general, these matrix elements 
are given by integrals over velocity-weighted products of the momentum transfer 
cross section and electron distribution function (isotropic part) . For the 
case of a simple power law dependence (see equations (5-2) and (5-3)), these 
matrix elements have been evaluated in SI. The purpose of the method derived 
in section 5.2. was to permit use of these results for arbitrary velocity 
dependence. From SI, equations (48) and (63d), the electron-neutral inter- 
action matrix is given by 


/ 


<v 


r(r-2) 

8 


r(r-2)(r-4) 

48 


r(r +2r+20 y 
16 


r(r-2)(r +2r+30) 
96 


r(r-2) 

8 


r(r +2r+20) 
16 


"22 

r(r 4 +4r 3 +68r 2 +8r+420) 

384 


r(r-2)(r-4) 

48 


r(r-2)(r +2r+30) 
96 


\ 


r (r 4 +4r 3 +68r 2 +8r+420) 


384 


33 


/ 


(5 -40a) 


where 


en 4 3 2 

H 22 = (r + 4r + 52r + 96r + 280 >/ 64 


(5 -40b) 


H 33 = (r6 + 6r5 + 142r4 + 528r3 + 3952r2 + 6864r + 15120) /2304.(5-40c) 
In SI this matrix is expressed in terms of "m", which is related to r by 


m = (r+3)/2. 


(5-41) 


The notation used above seems more convenient and compact. 
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Once the power law index r has been determined by the procedure in sub- 
section 5.2., all information necessary for the calculation of the transport 
tensor components is available in terms of the fundamental properties desired. 
These equations hold for the general case of arbitrary degree of ionization 
and magnetic field strength. For the ionosphere, there is a simplification 
associated with the fact that for electrons the magnetic field is always in 
the strong field limit (w^>> <v m > ) • This applies, of course, only to the 
transverse components. 


In the strong magentic field limit, transverse components of the transport 
tensors are expressed in the following manner (cf. equations (5-33) above): 
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CT 11 ± la 21 


N e / +i g <v >\ 
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(5-42b) 


(5-42c) 


(5-42d) 


from SI, equation (98) . The limiting expressions for the correction factors 
are then (SI, equations (99a) - (99d)): 


h^ = 1 for a = a, T, p, K 


g a = 1 

1 [ "-I + (r+2)<u m >/<u ei > 1 

8 T ' 2 L 1 + < V /<V ei > J 


(5-43) 

(5-44a) 

(5-44b) 
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(5-44c) 


'florn 


2 p * ?(^)<v.>/ < v. 1 >1 

>■ ‘ 5 L 1 + <v 7 'vi 7 J 

1/10 + ^2/5 + [ (r+5) (r+4)<v >/20<v .>] 
m • el 

®K 1 + <v >/<v 4 > 

m 61 


(5-44d) 


Since the transverse components can be evaluated from equations (5-42) - 
(5-44 ) % the lengthy procedure involving the interaction matrices can be avoided. 
This is not so for the parallel (3, 3) components; however, for this component 
there is a small simplification. Because this component is obtained by setting 
w to zero, it is seen from equation (5-36) that 

D 


**ij ^ <v ei ># » 


(5-45) 


in computing the parallel component only. With g correction factors given by 
equations (5-35a, c, e, g) , the parallel components are then computed from 
(cf. equations (5-33)):- 
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(5-46a) 

(5-46b) 

(5-46c) 

(5-46d) 


where the subscript and superscript o denotes no magnetic field. 
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To make use of these transport tensors In the context of the conservation 
equations, it is necessary to relate them to the unknown quantities, in parti- 
cular to the heat flux vector q^. In SI this quantity is denoted by and is 
shown to be related to the total energy flux H by (SI, p. 1656) • 

q a H - 5kT ?/2e . (5-47) 

^e e 

In the ionosphere it is frequently assumed that the current parallel to magnetic 
field lines vanishes. Physically this is interpreted as a steady state condi- 
tion in which an electric field has built up in such a way that it cancels 
current contributions due to gradients in temperature and density; this is 
termed the thermoelectric effect (Spitzer and Harm, 1953) . The effect in 
equation (5-47) is to make q^ equal to H. It is shown by Hochstim and Massel 
(1969, p. 160) that for this case the heat flux parallel to the magnetic field 
can be computed from 


and 


q e - 


eff 


VT 


where 


5N k T 

Ceff = K < V < 8 K ) eff 


_ 1 


^ 8 K^eff 


’K 


2g ; 


in notation introduced above. 


(5-48) 


(5-49) 


(5-50) 


In view of the approximations made in Section IV regarding transverse 
components in the electron equations, it is worthwhile to examine briefly the 
transport tensor expressions for consistency with those approximations. Let 
X denote any one of the transport tensors, and consider the ratios of trans- 
verse to parallel components. The component parallel to a transverse driving 
force is X^. From equations (5-42) for the strong field transverse case and 
(5-46) for the parallel (zero field) case, the ratio is 
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<v >' 



2 




(5-51) 


In SI a rather comprehensive set of tables is provided, giving g and h correc- 
tion factors for -3<r<3 with magnetic field strength ranging from zero to 
infinitely strong and degrees of ionization from essentially zero to completely 
ionized. The g correction factors are of order unity, always satisfying 
|g|$3. Disregarding them and using maximum values of <v > consistent with 
those used for order-of -magnitude estimates in Section IV gives 


21 

C 33 


10 
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(5-52) 


Similarly for the off diagonal component, corresponding to transport which is 
transverse to both the magnetic field and a transverse driving force, the ratio 
is 


x„. 


<V > 

21 
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X 33 


% 8 


-2 

< 10 , 


(5-53) 


In both cases the direction of the magnetic field clearly dominates the 
transport for isotropic driving forces. However, there may be no driving 
force in the parallel direction, as is frequently assumed to be the case for 
steady-state ionospheric electrostatic fields. The pertinent ratio is then 
the ratio of the transverse component parallel to the driving force (X^) to 
the component transverse to both the magnetic field and the driving force 
(X 21 ): 


X 11 


<v > 

x„. 


i g l 

21 
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10 


-2 


(5-54) 


In words, transport parallel to the driving force is negligible compared to 
that perpendicular to the driving force in the transverse plane. When the 
driving force is an electric field, the current (conductivity) parallel to 
the electric field (in the transverse plane) is called the Pedersen current 
(conductivity), while that perpendicular to the electric field is the Hall 
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current (conductivity) . It is readily recognized that only the Pedersen 
component represents work done by the electric field on the electrons. Since, 
according to equation (5-54), most electron transport transverse to the mag- 
netic field occurs through Hall motion. Joule heating (<*E • J) of electrons 
will be small. 

These considerations are consistent with the results of Section IV and 
this is evident from examination of equations (4-65a) and (4-66a) . In the 
momentum equations, the transverse and parallel component are uncoupled with 
only the Hall drift remaining in the transverse equation. Joule heating was 
dropped from the electron energy equation by comparison with thermal energy 
exchange terms. The additional conclusion that can now be drawn is that the 
spatial gradients in the electron energy equation (4-66a) need be taken only 
parallel to the magnetic field, since the transverse components of K are small 
compared with K~«. Mathematically, this is written (consistent with equation 
(5-48)) 

?’■<'■ V V “ r ll ( K ef £ -£)• «- 55 > 

where s is arc length along a magnetic field line (see Banks, 1966b). 

Equations presented in this part should be «. iequate for any ionospheric 
investigation of electron transport. 

5.4.2 ION Transport 

The formalism for a complete treatment of the ion transport problem, 
including arbitrary degree of ionization and arbitrary magnetic field strength, 
has been developed by S. T. Wu (1968); however, its form is impractical for 
numerical calculations. In the absence of a complete, usable method, Only 
what is needed for the conservation equations and is readily available in the 
literature is presented. The term "readily" is used advisedly, since Mitchner 
and Kruger (1973) note that with proper precautions, the Chapman-Enskog approach, 
as applied by Hirschf elder , Curtiss, and Bird (1954), can be used in computing 
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both ion and neutral transport coefficients (for the zero-magnetic field case). 
Some indication of the steps required for such a procedure are provided by 
Chapman (1954) in his evaluation of the viscosity and thermal conductivity 
coefficients for a fully ionized gas without a magnetic field. Rather than 
engage in such an effort, the results of Chapman (1954) are presented; 
limitations on their applicability and consequent effects on calculations are 
then examined briefly . 


References to equations (4-65c) , (4-66b) and (4-33) shows that only 
coefficients of viscosity and thermal conductivity are required for the ion 
conservation equations . In the previous section, it was seen that in a 
magnetic field the thermal conductivity becomes a tensor of rank two. 

Shkarofsky, et al. (1966, Chapter 8-8) show that viscosity is affected in a 
similar manner, although the results are included in the stress tensor rather 
than in an explicit tensor viscosity. They note that a scalar viscosity is 
applicable only if << <v i± >; this condition seldom holds in the ionosphere. 
Thus coefficients derived for a fully ionized gas in the absence of a magnetic 
field apply only for directions parallel to the field lines. For these condi- 
tions Chapman (1954) obtains the following expressions for viscosity and thermal 
conductivity coefficients of a single ion species: 


n 
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I El (kT) 5/2 
e*A 2 (2) 


(5-56) 


Kt 


. 15k 
4M ± n i 


(5-57) 


In the work cited, A^(2) is an integral based on a cutoff at the mean 
separation of particles. This has subsequently been modified in Chapman and 
Cowling (1970, p. 178) to use a Debye length cutoff, although the form of the 
result is unchanged. Thus A£(2) can be expressed in terms of the argument of 
the Coulomb logarithm. A, equation (5-20b) , as 

A 2 (2) = 2 [in (1+A 2 ) - A 2 /(l+A 2 )]. (5-58a) 
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Since A»1 (Spitzer, 1956; Banks, 1966a), this can be written to a good approxi- 
mation as 


A 2 (2) * 2 [ 2 £n A-l ] • 

For consistency with the results for electron-ion collisions, this should fur- 
ther be approximated by 

A 2 (2) « 4 An A . (5-58b) 

In the ionosphere lnA>10 (see Spitzer, 1956). It may now be noted that Chapman 
and Cowling (1970) use an expression corresponding to A which differs from 
equations (5-20) by a factor 4/3. Such a difference is inconsequential in 
terms of the approximation leading to equation (5-58b) . 


Equation (5-56) can thus be rewritten as 
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» 


(5-59) 


while equation (5-57) remains unchanged. This result may be compared wd^th a 
similar result obtained by Shkarofsky et al. (1966, equation (8-140)). That 
form can be shown to be the same with a factor 3/ [4 \fl{ . 7326) ] in place of 5/8 
in equation (5-59). These numbers have decimal values of .7239 and .6250 
respectively. However, the expression obtained by Chapman and Cowling (1970) 
is in the first approximation (Chapman - Enskog method of successive approxi- 
mations) . They note that in the second approximation, the expression for the 
first approximation is simply multiplied by 1.15. When this is done, the 

numerical constant is .7188, in very close agreement with the result of 

T 

Shkarofsky et al. (1966). The corresponding correction factor for in the 

second approximation is 1.25. So the recipe for accurate ion transport 

coefficients is to compute n. from equation (5-59). With this value, evaluate 
T 1 

; then multiply the former by 1.15 and the latter by 1.25. 
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As indicated at the outset, these expressions are of limited applicability. 

The question is, how important are the limitations? Consider first the effect 

of the geomagnetic field. Chapman and Cowling (1970, Chapter 19) find that 

transport perpendicular to the magnetic field is reduced from the zero field 

2 

case by a factor 1/ [ 1 + (w.t) ] for direct transport and by a factor (a> t)/ 

2 1 l 

[1 4* ] for transverse transport. The terms direct and transverse, as 

used in this sense, refer to directions parallel to and perpendicular to the 

driving force (external force or gradient) in the plane perpendicular to the 

magnetic field, t is a relaxation time which is given approximately by the 

inverse of the total ion collision frequency. Thus at low altitudes (<130km) 

where direct transport is little affected by the magnetic field, while 

transverse transport varies approximately as to./v. . At high altitudes, on 

— li _ _2 

the other hand, where o)j» direct transport varies as (uj^/v^ and 

transverse transport as Hence, at high altitudes transport is 

confined to directions parallel to the magnetic field (except for strong trans- 
verse electric field driving forces) . 

For thermal conduction, observation and theory agree that ions have 
approximately the neutral temperature up to altitudes of 250km to 300km. So 
below that transition, thermal conduction is unimportant for the ions. Above 
this transition so that ion thermal conduction transverse to the geo- 

magnetic field is unimportant. Thus, the thermal conduction coefficient for 
the zero-magnetic-field case is sufficient for practical purposes. Further- 
more, above the thermal transition region, ion-neutral collision frequencies 
are small compared with collision frequencies of ions with charged particles. 
Hence, the ion thermal conductivity for a fully ionized gas is appropriate. 

In sum, the limitations on the application of the ion thermal conductivity 
coefficient, given by equation (5-52), should have little practical effect 
on ionospheric calculations. 

For viscosity, the physical situation is rather different and not so 
well-defined, since less is known about ionospheric velocity fields than 
about its thermal structure. According to present understanding, signifi- 
cant ion flow velocities have two origins: electric fields and neutral 
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flow (via collisions). As discussed in Sections I and IV, electrostatic fields 
in the ionosphere are predominately transverse to the geomagnetic field. At 

•> 4 - 

high altitudes, where co i >> v^ 9 the Hall drift (EXB) results (e.g. electron 
momentum equation (4-36)). At low altitudes, where v^» ions essen- 
tially move with the neutral flow. A transition between these motions occurs 
at intermediate altitudes where = 1 ^ (~120km - 140km), resulting in vertical 
velocity gradients. Since large electric fields are observed primarily at 
high latitudes where geomagnetic field lines are almost vertical, large ion 
drifts due to electric fields are primarily horizontal. Thus, in this case, 
the flow is mainly horizontal with primarily vertical gradients which are 
largest in the region where This qualitative picture is given quanti- 

tative substance in the next section. 

When neutral flow causes the ion motion, a similar situation results. 
Collisional coupling at low altitudes ( v ^ n >> w^) causes ions to move with 
the (horizontal) neutral wind. With increasing altitude the ion coupling 
undergoes a transition from collision dominated to geomagnetic field dominated, 
until at high altitudes the ions move primarily along magnetic field lines in 
the absence of transverse electric fields. Although this picture is somewhat 
complicated by geomagnetic field geometry variation with latitude, the physical 
results are basically the same as above. Velocity gradients are predominately 
vertical and strongest in the dynamic transition region • 

Physically, viscosity represents momentum transport due to spatial grad- 
ients in velocity. From the discussion above, this corresponds essentially 
to the transport of horizontal momentum in vertical directions. Since the 
physical conditions for maximum velocity gradients are now more precisely 
defined ( v in ~ a) j[) than those used for order-of-magnitude estimates in Section 
IV, a new estimate may be useful. In equation (4-39a) the same values are 

used to obtain* a maximum velocity gradient (~.l s ^) , but the collision fre- 

2 -1 

quency is now restricted to the order of the gyrofrequency (-10 s ). 

Under these conditions the ratio of the viscous stress tensor to scalar 
pressure is 

I -j^^l < 10' 3 . (5-60) 
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Hence viscous effects on ion motion are negligible in the dynamic transition 
region and below. 

At high altitudes, where the total ion collision frequency is much reduced, 
it might appear that viscosity could be important. However, since the primary 
driving force for ion flow at high altitudes is electric fields, this is 
unlikely for the following reasons. The largest measured electric fields occur 
at high latitudes; that is, they are of magnetospheric origin and mapped into 
the ionosphere along geomagnetic field lines (see Figure 1-4) . Because magnetic 
field lines are approximately equipotentials and the magnetic field is almost 
uniform over the altitude range of interest here, the electric field varies 
little along a field line over this range. Correspondingly, there is little 
variation in the ExB drift velocity, resulting in very small velocity gradients. 
Consequently, along with the reduction in collision frequency at high altitudes, 
there is a reduction in the vertical velocity gradient, due to a uniformity in 
altitude of the driving force. So viscosity is unlikely to be important for 
ions at either high or low altitudes. If it is significant, it will probably 
be at high altitudes and for gradients along field lines, conditions for which 
the viscosity coefficient above is appropriate. 

5.4.3 Neutral Transport 

Neutral gas transport is not complicated by considerations of magnetic 
fields or degrees of ionization. Moreover, the kinetic theory for transport in 
single species monatomic gases is well-established (see Chapman and Cowling, 
1970). However, difficulties arise in gas mixtures, particularly mixtures 
involving polyatomic molecules. These difficulties are discussed in some 
detail by Chapman and Cowling (1970, chapters 12-14) and will not be elaborated 
here. Using approximate formulas and some inspired curve-fitting, Banks and 
Kockarts (1973, chapter 14) have devised simple, yet reasonably accurate 
expressions for the viscosity and thermal conductivity coefficients for the 
mixture of neutral gases in the upper atmosphere. These are presented below. 

The temperature dependence of the viscosity coefficients of all primary 
neutral species has been approximated in a common power law form, based on 
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experimental data for all except atomic hydrogen. Thus all the viscosity 
coefficients vary as T'^ for 200 o K<T<2000°K, a range adequate for most upper 
atmospheric conditions. Because of this common temperature dependence. Banks 
and Kockarts (1973) have been able to reduce the viscosity for the mixture to 
the following simple form: 

n = A T* 69 , (5-61a) 

n 

where 

A = y A n fl n . (5-61b) 

L n n L n 
n n 

The summation runs over all neutral species, and the A coefficients are given 

n 

in Table 5-4. 


Table 5-4. COEFFICIENTS FOR COMPUTING THE NEUTRAL 

VISCOSITY COEFFICIENT (BANKS AND KOCKARTS, 
1973) 


NEUTRAL 

SPECIES 

A n (10"®gm cm"^ s"^) 

n 2 

4.03 

°2 

3.43 

0 

3.90 

H e 

3.84 

H 

1.22 


Thermal conductivity is treated in a similar manner. If all upper atmo- 
spheric gases were monatomic, a relation similar to equation (5-57) for ions 
could be used to obtain the thermal conductivity coefficient from the viscosity 
coefficient. However, the internal degrees of freedom associated with the 
diatomic molecules and 0^ also affect the thermal conductivity. Analysis of 
experimental data allows correct inclusion of these effects. Banks and Kockarts 
(1973) are again able to express all temperature variations in a common power 
law form, permitting the thermal conductivity for the mixture to be written 
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K T - BI - 69 
n 

where 

b - 1 b „"»^ % • 
n n 

The B n coefficients are given in Table 5-5. 

Table 5-5. COEFFICIENTS FOR COMPUTING THE NEUTRAL 
THERMAL CONDUCTIVITY COEFFICIENT (BANKS 
AND KOCKARTS, 1973) 


NEUTRAL 

SPECIES 

B n (erg cm - ^ s~^ °K -1 ) 

n 2 

56.0 

CM 

O 

56.0 

0 

75.9 

H e 

299.0 

H 

379.0 


These relations may be used over the temperature range 200°K<T n £2000°K. Although 
the error in 0 ^ due to curve fitting approximations approaches 20 percent at 
2000°K, at altitudes where is important, temperatures and errors are much 
lower. Errors for the other gases are smaller than for 0^. 

Since viscosity and thermal conductivity are the only transport coefficients 
appearing in the conservation equations, none others are required for closure of 
the neutral equations. 


(5-62a) 
(5 -6 2b) 
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Section VI 


NEUTRAL WINDS IN THE AURORAL E REGION 
DURING GEOMAGNETIC DISTURBANCES 


The conservation equations developed in previous sections are somewhat 
general in the sense that only those approximations appropriate to all E and F 
region conditions have been applied. For the particular physical circum- 
stances of any specific study, further approximations may be justified. More- 
over, observational data may be used to specify some variables so that certain 
of the equations need not be solved at all, or may be greatly simplified. 

Also, numerical constraints may require simplifying assumptions which limit 
the applicability of the results to special circumstances. For these reasons, 
results of the previous sections must be regarded as appropriate starting 
points for a wide variety of upper atmospheric investigations, rather than as 
finished products. The present section is an example of one such investigation, 
and it illustrates various methods for bridging the gap between the initial 
set of conservation equations and the final working equations which can be 
solved numerically within computational constraints. 

6.1 INTRODUCTION 

Since the incoherent scatter radar facility at Chatanika, Alaska, became 
operational in 1971, thermal, dynamic, and electrical properties of the iono- 
sphere above that location have been observed under a variety of conditions 
(e.g. Watkins and Banks, 1974; Brekke, et al. , 1973, 1974a, b). Of particular 
interest are winds generated during geomagnetic disturbances. It is well 
established that global disturbances of both the charged and neutral components 
of upper atmospheric gases are associated with geomagnetic storms (Jacchia, et 
al., 1966; DeVries, 1972; Matuura, 1972; Wu, et al. , 1974). The mechanisms 
by which these disturbances occur at low and middle latitudes remain under 
study; but neutral winds originating in the auroral region are proposed in 
some theoretical explanations (e.g., Obayashi and Matuura, 1972). Such winds 
have apparently been observed in the E region at low latitudes following a 
magnetic storm (Smith, 1968), and simple estimates, based on observations at 
high latitudes, support the plausibility of such occurrences (Rees, 1971a). 
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Hence, extended observations of these neutral winds in the vicinity of their 
origin, may provide decisive information in establishing the validity of these 
and other theories. 

In this study, attention is directed toward certain aspects of the neutral 
winds as determined by incoherent scatter radar observations. Because ion 
velocities, rather than neutral velocities, are observed directly by the 
radar, the neutral wind velocities must be derived from the ion velocities. 

One of the approximations required to determine the neutral velocity is that 
this velocity is uniform over an altitude range of 25 km to 50 km in the E 
region. Although Brekke , et al. (1973, 1974a) recognize that this condition 
does not hold, no quantitative estimates of its effects on the derived neutral 
winds are provided. One of the objectives here is to examine the effect of 
the neutral velocity altitude structure on neutral velocities derived from the 
incoherent scatter radar data. 

A second objective concerns the origin of neutral winds during geomagnet- 
ically disturbed conditions. Brekke, et al. (1974a) note that incoherent 
scatter radar data alone provides insufficient information to separate the 
effects of ion drag (collisional coupling) and Joule heating on neutral winds. 
From observations of high latitude rocket-borne chemical releases at E region 
altitudes and ground-based magnetometer data, Rees (1971a) concludes that the 
neutral winds are generated by collisional coupling with electric-field-driven 
ions. This view is shared by Meriwether, et al. (1973), based on similar data 
for altitudes above 200 km. On the other hand, from chemical release data 
below 180 km, Ro.thwell, et al. (1974) concur in the conclusion of Stoff regen 
(1972), based on similar data for higher altitudes, that during geomagnetic 
disturbances, the deviations of observed neutral winds from theoretical models 
are due to high latitude heating and resultant pressure gradients. 

To examine this problem, an approach which is partially direct and par- 
tially indirect is employed. In the direct approach, the effects of collisional 
coupling are computed by numerically solving the ion and neutral momentum 
equations simultaneously, based on observed electric fields as the driving 
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force and omitting pressure gradients. If there is good agreement between the 
calculated and observed velocities, this indicates that collisional coupling 
is adequate to explain the observations and that effects of pressure gradients 
are inconsequential. On the other hand, large discrepancies would indicate 
the opposite. Then, the indirect approach is employed: the discrepancy 

between the observed and calculated velocities is used to estimate the omitted 
driving force (pressure gradient) . 

Data used to carry out the calculations and to provide the observational 
reference velocities are described in, subsection 6.2. In subsection 6.3 
details of the theoretical and numerical methods are presented. Model calcula- 
tions are presented in subsection 6.4 to demonstrate in a simple context the 
effects to be encountered under realistic conditions with observed electric 
fields. In subsection 6.5 a series of calculations, designed to meet the 
objectives outlined above, are discussed in detail. Finally, results of this 
study are discussed and summarized in subsection 6.6. 


6.2 OBSERVATIONAL DATA 
6.2.1 Characteristics of the Data 

Data used in the calculations and comparisons below were obtained by the 
incoherent scatter radar facility at Chatanika, Alaska on 15 May 1974. This 
facility became operational in July, 1971; and early observations have been 
reported recently in the literature (Doupnik, et al. 1972; Brekke, et al. , 
1973, 1974 a, b; Banks, et al. , 1973; Banks, et al. , 1974; among others). 

Leadabrand, et al. (1972) have described the equipment and its operating 
characteristics; and Banks and Doupnik (1974) have recently reviewed the 
experimental methods and operational modes used with the radar. Therefore, 
only those aspects pertinent to the analysis of the present data need be 
summarized here. 

Location and geomagnetic field characteristics of Chatanika, Alaska, are 
given in Table 6-1. Operating at a frequency of 1290 MHz, the Chatanika. inco- 
herent scatter radar has a completely steerable dish antenna with a beam width 
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of .6 degrees. For the experiment on 15 May 1974 the recently developed 
azimuth scan procedure (Banks and Doupnik, 1974) was employed. Elevation of 
the antenna was held constant at 76.5 degrees above horizontal and the antenna 
azimuth was changed at the rate of 1 degree per second. Data were integrated 
over 15 second intervals (15 degrees in azimuth) to obtain estimates of ion 
drift velocity vectors and electron density profiles. The 24 estimates, for 
each 360 degree scan, were then analyzed, using a Fourier analysis technique, 
to obtain the ion velocity vector in 8 separate range gates and a time- 
averaged electron density profile for the 6 minute time period. At azimuth 
299 degrees (geomagnetic west) on each scan, the direction of antenna rotation 
was reversed. 

Table 6-1. CHATANIKA, ALASKA: LOCATION AND GEOMAGNETIC 

FIELD PROPERTIES 


LOCATION: Geographic - 65.1°N, 147.45°W 

Geomagnetic Dipole - 64.75°N, 105°W 

TIME: Alaskan Standard Time (AST) = UT-10 hours 

GEOMAGNETIC FIELD: Dip Angle (I) = 76.5° 

Magnetic Declination (6) = 29° East of 

Geographic North 

Magnetic Field Intensity (B) = .54 Gauss 


A diagram of the azimuth scan method of sampling the ionosphere is pre- 
sented in Figure 6-1. It is seen that in addition to time-averaging over 6 
minutes, this technique averages spatially. For example, the centers of the 
first two range gates correspond to altitudes of 109 km and 167 km; given the 
elevation angle, the radar beam samples over horizontal distances of 55 km and 
85 km at these altitudes. In addition, due to the time length (x) of the 
radar pulse, the radar signal averages ionospheric properties over the spatial 
extent of the pulse (ex). Brekke, et al. (1973) note that although the total 
extent of each range gate is 96 km, the effective length is reduced to 25 km 
to 50 km due to weighting of the radar signal. This weighting has three 
components: (1) a triangular factor which is unity at the center of the range 

gate, decreasing linearly to zero at the ends; (2) a factor which decreases 
as the inverse square of the range; and (3) the electron number density. The 
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Figure 6-1. DIAGRAM OF THE AZIMUTH SCAN MODE OF OBSERVATION FOR THE 
INCOHERENT SCATTER RADAR AT CHATANIKA, ALASKA (BANKS 
AND DOUPNIK, 1974) 


first two of these multiplied together are designated by W(z) and are pre- 
sented graphically by Brekke, et al. (1973); those results are used here with 
appropriate translation to the proper altitude. 

Interspersed with .the long radar pulses, used for velocity measurements, 
were short pulses (with greater range resolution) for electron density deter- 
mination. Beginning at a range of 73.6 km (71.6 km in altitude w due to the 
elevation angle), successive values are spaced at range increases of 4.5 km 
(4.38 km altitude). However, since the signal samples a range of 10 km, 

adjacent electron density values are not totally independent (Doupnik, 1975). 

4 -3 

Doupnik (1975) also notes that for electron densities below 10 cm , true 
densities are underestimated by the radar because the Debye length becomes 
comparable to the radar wavelength. In addition, measured densities must be 
multiplied by a factor of (1 + T e /T\)/2 to obtain true values. For E region 
altitudes, T^/T^ = 1 may be a reasonable approximation, however, not for 
higher altitudes. Since present interest is in the E region (and in the 
absence of temperature measurements), the electron density measurements are 
used without modification. 


Geomagnetic activity during the experiment was relatively high globally; 

15 May 1974 was designated one of the! 5 disturbed days for that month (Lincoln, 

1974). Values of the global three-hour geomagnetic range index K varied from 

P 

3- to 5+, averaging about 4. These, together with local index K for College, 
Alaska, are shown in Table 6-2. 


Table 6-2. 


GLOBAL (K p ) AND LOCAL (K) (COLLEGE, ALASKA) 
THREE-HOUR GEOMAGNETIC RANGE INDICES FOR 
15 MAY 1974 


PERIOD: 

1 

2 

3 

4 

5 

6 

7 

8 

K p 

4+ 

3- 

4- 

5- 

5+ 

4+ 

3 

4+ 

K 

4 

2 

3 

6 

7 

6 

3 

3 


6.2.2 Data Analysis 

From the ion flow velocities measured in different range gates, values of 
the electric field and the altitude-averaged E region neutral wind can be 
determined. The procedure for accomplishing this is described in some detail 
by Brekke, et al. (1973); it will be summarized here. This procedure is 
followed exactly since it provides the currently-accepted observed values from 
the incoherent scatter radar method. Part of the objective of this study is 
to compare theoretical calculations with standard observed values; so in 
treating the data, the established procedures are used without deviation. 


At high altitudes where ok >> , ions move with the Hall drift velocity 

(E x B c/B^) (cf. electron momentum equation (4-36)). The vector equation 
relating the ion velocity to the electric field can then be solved for the 
electric field in terms of the ion velocity at high altitudes. At low alti- 
tudes where w. < v. , collisions with neutrals as well as the electric field 
l in 

affec't the ion motion. With the electric field known from the high altitude 
ion velocity measurements, the ion velocities observed at low altitudes can be 
used to determine the neutral wind velocity. 


These qualitative considerations can now be put on a quantitative basis. 
If gravitational and pressure gradient terms are assumed to be small compared 
with electric and magnetic field terms, the ion momentum equation (4-40) can 
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be written in the form 


in,-*- -+ s . E B 

~ (v n ' V + ¥ C + V i X B = 0 




( 6 - 1 ) 


From this equation it is apparent that if >> the electric field com- 

ponents transverse to the magnetic field can be obtained from ion flow 
velocities from the equation 


E + — x B - 0 . (6-2) 

c 

Brekke, et al. (1973) find that the ion velocity measured in the second range 
gate (centered at 167 km) . shows little, rotation from those at higher range 
gates. They conclude that the electric field can be adequately determined 
from equation (6-2), using, observed ion velocities from the second range gate. 
In the model to be discussed later (similar to that used by Brekke, et al. 
(1973)), it is found that at this altitude v in / (JLj i ~ -05, thus the collision 
term in equation (6-1) is expected to be unimportant, consistent with the 
observations. 


By taking scalar and cross products of S with equation (6-1) , that equa- 


tion can be solved explicitly for v^ with the result 
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where the component of 2 parallel to S has been assumed to be negligible com- 
pared to the components transverse to $. Equation (6-3) holds at each altitude 
however, observed ion velocities are weighted averages over altitude. There- 
fore, similar averaging must be applied to this equation. 


Based on the weighting of the incoherent scatter radar signal, discussed 
in the previous section, Brekke, et al. (1973) determine that the weighted 
averages can be represented as 


v z 2 Z 2 

= [ N g (z) W(z) v ± (z) dz/f N 
Jz. Jz. 


(z) W(z) dz 


(6-4) 
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Limits of integration correspond to the boundaries of the range gate; this is 
accomplished in effect by the vanishing of the weighting function W(z) at 
these boundaries. This weighting function represents the linear decrease 
from the center of the range gate divided by the square of the range; it is 
shown for the first range gate in Figure 6-2. It is a consequence of this 
weighting by W(z) and N^(z) that the effective width of the first range gate 
is reduced from 96 km to 25 km to 50 km. 



Figure 6-2. INCOHERENT SCATTER RADAR WEIGHTING FUNCTION W(Z) 

(BREKKE , ET AL. 1973) 

In order to apply this averaging technique to the right side of equation 
(6-3) and obtain useful results, Brekke, et al. (1973) make the following 
assumptions: (1) the neutral wind velocity has a uniform value ^ throughout 

the first range gate; (2) the electric field determined from the second range 
gate is mapped down magnetic field lines to the bottom of the first range 
gate with no change in either magnitude or direction. Under these conditions, 
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equation (6-3) can be multiplied by N^Cz) W(z) , integrated over altitude and 
normalized by the denominator on the right side of equation (6-4) to obtain 


v 


K V 

K 3 n 


C 1 V 

1 n 


+ K • $)S/B 2 . 

L n 


(6-5) 


The K-coefficients are given by 
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(z) W(z) dz 


(6-6a) 


(6-6b) 


(6-6c) 


(6-6d) 


The ion velocity on the left side of equation (6-5) now represents the velocity 
as observed by the incoherent scatter radar. 



To determine the neutral velocity, equation (6-5) is solved explicitly for 
in the same manner as equation (6-3) was obtained, with the result 


n 




t + t \ 

¥ + ¥ c ) 


(<i - < 2 k 3 ^ 



[f X ! c + (^ 



(6-7) 


(Equation (6-7) corrects a typographical error in equation (11) of Brekke, 
et al. (1973).) From this equation the height-averaged neutral wind is deter- 
mined from incoherent scatter radar data. Ion velocity ^ is obtained from 
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the first range gate; E is obtained from ion velocities in the second range 
gate (through equation (6-2)); and the x -c °ef f icients are evaluated from 
electron number density measurements and a collision frequency model (to 
be discussed later). 

By means of these techniques, Chatanika incoherent scatter radar observa- 
tions made on 15 May 1974 have been analyzed. Resulting electric field values 
are shown in Figure 6-3a; and horizontal neutral velocities, determined from 
equation (6-7), are shown in Figure 6-3b. It should be noted that the x- 
component is defined positive to the south, while the meridional component 
in the figures of Brekke et al. (1973, 1974a) is defined positive to the 
north. Values shown in Figures 6- 3a, b are used in the calculations and com- 
parisons below. 

The electric field components in Figure 6-3a show two distinct periods of 
enhancement, the first extending from about 06 to 11 hours UT and the second 
from about 12 to 18 hours UT. During the calm periods, electric fields have 
values generally less than about 10 mV m \ while during disturbances they 
have values of several tens of millivolts per meteri Rapid variations in 
magnitude are the rule; but reversals of direction which persist with the 
large magnitude electric fields, occur only once, shortly before local midnight 
.(1000 UT) . This reversal is a general feature of this local time period. It 
is associated with the Harang discontinuity (Heppner, 1973), originally observed 
as a reversal of ion flow direction. 

Neutral velocities shown in Figure 6-3b are notable for their rapid and 
sometimes quite large fluctuations in both magnitude and direction. Since 
this indicates large variations during the peirod over which the data are 
averaged, large errors are likely (Brekke et al. 1973, 1974a). Nevertheless, 
the trend of the data is clear: except for occasional large fluctuations, 

magnitudes of the neutral velocity components are generally near or less than 
about 125 m s \ Periods of greatest velocity fluctuations correspond to 
periods of enhanced electric field. These results will be examined more 
closely in the comparisons with calculated velocities. 
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Figure 6-3. ELECTRIC FIELDS (a) AND NEUTRAL WINDS (b) IN GEOGRAPHIC 
COORDINATES, AS DETERMINED FROM INCOHERENT SCATTER RADAR 
OBSERVATIONS AT CHATANIKA, ALASKA ON 15 MAY 1974 
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6.3 THEORETICAL APPROACH 


6.3.1 Physical and Mathematical Description 

As noted previously, the effects of collisional coupling in neutral wind 
generation are to be determined. In simplest terms, an external electric 
field causes the ions to move. Collisions of ions with neutrals both modify 
the ion motion and cause the neutrals to move. It is the resultant motions of 
the neutrals, as functions of altitude and time, that are of particular inter- 
est. Model calculations of this nature have been presented by Fedder and 
Banks (1972) • However, they differ from the present calculations in two 
important respects: theirs neglected coriolis effects and they were only 

model calculations with no explicit comparison with experiment possible. The 
first of these differences is examined in similar model calculations in sub- 
section 6.4 to determine the importance of the coriolis term. The second 
difference, use of observed electric fields and comparison of corresponding 
calculated velocities with observed velocities, allows the importance of this 
collisional coupling in neutral wind generation to be assessed. 

The theoretical framework has been developed in previous sections. Here, 
in the application of this framework, it remains to apply those additional 
approximations which the more restrictive conditions of a particular problem 
permit. And then a method for determining each remaining unknown quantity must 
be specified, either by including an appropriate equation in the system to be 
solved, or by adopting appropriate model values. This translation of the 
physical problem to a mathematical and numerical problem is accomplished 
below. 


The coupled ion and neutral equations of motion (4-65c) and (4-65d) are 
the appropriate starting points. In consonance with discussions in Section V, 
the ion pressure tensor term is neglected in comparison with the scalar pres- 
sure term. Thus, the initial equations are 
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in Cartesian tensor notation (equation (4-33)). 


(6-8b) 


(6-8c) 


In addition to being coupled to each other, these equations are coupled 
to the continuity and energy equations through mass densities (p^, p^) f press- 
ures (p^, p^) and collision frequencies (v^ n , v ni ) . Thus, additional infor- 
mation or equations are required in order to solve equations (6-8) . From the 
incoherent scatter radar observations, electron number density profiles are 
available* Ion composition cannot be obtained from these observations, so an 
ion composition model, based on the rocket-borne mass spectrometer measurement 
of Kopp, et al. (1973) at high latitudes (68°N) , is used. With this information, 
the ion continuity equation is no longer required. Since electric fields and 
velocities are generally large throughout the experiment, the order- 
of-magnitude estimates in Section IV indicate that the pressure gradient and 
gravitational terms can be dropped from the ion equation (6-8a) . In addition, 
at the altitudes of interest (90 km to 200 km) the ion temperature can be 
assumed equal to the neutral temperature (required only for collision fre- 
quency calculations). These approximations and assumptions obviate the use of 
the ion energy equation. 

To supplement the neutral equation of motion a model atmosphere is adopted 
to provide neutral number densities and temperatures. This approach is taken 
partly to keep the computational effort reasonable and partly due to lack of 
information. Since temperature measurements were not made in the experiment 
of 15 May 1974, there is no means for monitoring any theoretical calculations 
of temperature. Also, by using a static model for neutral temperature and 
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densities, the source of calculated velocity variations can be more easily 
identified. For a series of calculations of increasing complexity, this 
approach is a reasonable first step. 

The neutral atmospheric model adopted is the Jacchia 1971 model (Jacchia, 
1971; also cf. CIRA 1972) for an exospheric temperature of 1000°K (see Figure 
1-3). Although this choice of exospheric temperature is arbitrary, it does 
correspond to the exospheric temperature used by Brekke, et al (1973, 1974a) 
in their data analysis (however, they employed a neutral model by Banks and 
Kockarts (1973)). Effects of model differences will be discussed later. It 
may be noted that recent temperature measurements reported by Watkins and 
Banks (1974) indicate that 1000 °K is a reasonable choice for the exospheric 
temperature for a problem such as this. Their observations show no pronounced 
diurnal variations; however, temperature excursions of 100°K to 300°K over 
time periods of an hour, associated with large electric fields, are common. 

As a result of using a model atmosphere, the continuity and energy equa- 
tions are no longer required. However, there is an additional consequence for 
the neutral momentum equation. The Jacchia 1971 model is a static diffusion 
model' in which it is assumed that each neutral constituent is in diffusive 
equilibrium. For present purposes the important feature of this condition is 
that the vertical pressure gradient exactly balances the gravitational accel- 
eration in the momentum equation - there is no vertical motion. In reality, 
small vertical velocities are observed (Rieger, 1974). However, because the 

gravitational and pressure gradient terms are larger than the velocity terms, 

2 3 

by factors of 10 to 10 under these conditions, a relatively small adjustment 
of the pressure gradient can counteract vertical motions, causing them to 
vanish or remain quite small. Consequently, vertical motions are suppressed 
in equations (6-8 b) and (6-8c) : dif fusive equilibrium is assumed. 

Finally, all horizontal gradients are neglected. This is required by the 
fact that all observations are made at a single location: no information on 

the horizontal variation of properties is available. Implicit in this assump- 
tion is horizontal uniformity, which is a bad assumption at high latitudes. 
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However, horizontal gradients affect only certain terms and if these terms ^txe 
small relative to others, the consequences may be negligible. One of the 
objectives of this calculation is to test this assumption. 

• L 

With the approximations and assumptions made, the ion equation of motion 
in equation (6-8a) reduces to 



which is equivalent to equation (6-1). The neutral equation (6-8b) is corre- 
spondingly reduced to 
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where the form of remains as given in (6-8a) , but many of the terms vanish. 


To solve equations (6-9) and (6-10) numerically, they must first be 
decomposed into scalar component equations. Based on physical considerations, 
it is convenient to treat these equations in two different Cartesian coordin- 
ate systems; the computer can easily transform the coordinates numerically as 
required for calculations. The ion equation is treated in a system designated 
here as "plasma" coordinates (1, 2, 3), in which the 3-component is parallel 
to B, the 2- component is toward geomagnetic east in the plane transverse to B, 
and the 1-component lies in the geomagnetic meridian plane; this is shown 
in Figure 6-4a. The neutral equation is decomposed in a local Cartesian coordin- 
ate system (x, y, z) oriented along geographic coordinates with x positive 
to the south, y positive to the east, and z positive vertically upward, as shown 
in Figure 6-4b . Angular relations between the coordinate systems are shown in 
Figure 6-4c. Transformation equations between these coordinates for an arbitrary 
vector q are: 

(i) From plasma coordinates (1, 2, 3) to local geographic coordinates 
(x, y, z) 

q x = q 2 sin tS - cos 6(q^ sin I + q^ cos I) (6-lla) 
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PLASMA COORDINATES LOCAL GEOGRAPHIC COORDINATES 
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(b). 


LOCAL HORIZONTAL 




GEOGRAPHIC 

NORTH 


GEOMAGNETIC 
► NORTH 


(C) 


Figure 6-4. DEFINITIONS OF COORDINATE SYSTEMS AND ANGLES RELATING THEM 
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(6-llb) 


q y “ q 2 008 5 + 8in sin 1 + q 3 cos ^ 


cos I - q^ sin I 


(6-llc) 


(ii) From local geographic coordinates (x, y, z) to plasma coordinates 
(1, 2, 3) 


q- * q cos I - sin I(q cos <5 - q sin 6) 
X z x y 


q 2 - q x sin 6 + qy cos 6 


q, ■ -q sin I - cos I(q cos 6 - q sin 6) 
j z x y 


(6-12a) 


(6-12b) 


(6-12c) 


Consistent with the assumptions of Brekke, et al. (1973, 1974a) and dis- 
cussions in Sections I and IV, the electric field is assumed to be uniform, and 
the component parallel to the geomagnetic field is assumed to be negligible. 

In the plasma coordinate System this is given by 


e 3 - 0 . 


(6-13a) 


If this is used in equation (6-12c) , E can be expressed in the terms of E 

Z X 

and E as 

y 


E - - cot I(E cos 6 - E sin 6) 
z x y 


(6-13b) 


When the ion momentum equation (6-9) is written out in plasma coordinates, 
the directional coupling due to the Lorentz force can be removed algebraically 
giving the explicit results 
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The neutral velocity component equations in local geographic coordinates are 
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(6-15a) 

(6-15b) 

(6-15c) 


where X is geographic latitude. £o obtain the viscous terms in this form, it 
is assumed that the viscosity coefficient varies slowly in altitude compared with 
neutral velocity variations, an assumption which has been found to be valid. 


The set of neutral velocity component equations is seen to be a set of 
coupled (in the coriolis terms) partial differential equations, first order in 
time and second order in space. In addition to being coupled among themselves, 
these equations are coupled to the ion equations through the collision terms. 
Although the ion equations are of a time- independent form, the ion velocities 
are implicit functions of time through the time dependence of the electric 
field and the neutral velocities. Since S and (which varies directly with 

N^) are specified outside the system of equations (from scatter radar observa- 
tions) , numerical solution is the only feasible approach. 

Collision frequencies are computed according to equation (4-14c) and the 
species collision frequencies presented in subsection 5.3.3. Because static 
models of the neutral atmosphere and ion composition are employed, the ion- 
neutral collision frequencies v^ n remain fixed throughout the calculation. 
However, as noted above, \> ni varies with and must be updated from observation. 
The viscosity coefficient is evaluated from equations (5-61a, b) . Since it 
depends only on the neutral model, the viscosity altitude profile is constant 
throughout the calculation. 
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6.3.2 Numerical Methods and Boundary Conditions 

To solve equations (6-61a, b) an explicit finite difference technique is 
employed. For this discussion, let subscript m designate the spatial grid 
index and superscript n the time grid index, as shown in Figure 6-5. The 
partial derivative of a variable V with respect to time is approximated by the 
finite difference expression 


BV 

at 


V n+1 

m 



At 


(6-16) 


where At is the uniform grid spacing in the t direction. Similarly, the 
second partial derivative of V with respect to z is approximated by 

i_/i»\ y? - 2g „“ + y; 

81 l 8 V <4Z) 2 

where Az is the grid spacing in the x-direction. 


To use these expressions, equations (6-15a, b) are first placed in the 


form 


3v 3 2 v 

= Q(z) Tp- + R(v , v , v. , z) 

at 2 nx ny ix* 
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(6-18a) 
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where the functions Q, R, S can be identified by comparison with equations (6- 
15a, b) . When the derivatives in equations (6-18a, b) are replaced according 
to equations (6-16) and (6-17) , the resulting equations can be placed in the 
form 

(v ) n+1 = (v ) n H ^ n Q(z ) [ (v ) n .-| “ 2(v ) n + (v ) n - ] 

nx m nx m , .2 m N nx m+1 nx m nx m-1 

(Az) 


+ At R((v ) n , (v ) n , (v. ) n+1 , z ) 
nx ra ny m * v lx m * m 


(6-19a) 
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Evaluation of equations (6-19a, b) requires the specification of initial 

» 

and boundary conditions. For initial conditions, the neutral velocities are 
taken to be zero. Consequences of this assumption will be discussed later. 
Because a second spatial derivative is involved j two boundary conditions must 
be specified. For finite difference equations of the form in (6-19), it is 
appropriate to specify one condition at each boundary in z. Thus, the lower 
boundary is taken to be 90 km and the neutral velocity is assumed to vanish 
there. At the upper boundary, z = 250 km, the vertical derivative of the 
neutral velocity is assumed to vanish. Tests involving changes of these bound- 
ary conditions have shown that the lower boundary condition has negligible 
effects at 100 km and above, while the upper boundary condition has negligible 
effects on results at and below 200 km. These results are consistent with 
similar results reported by Fedder and Banks (1972) , who find that significant 
boundary condition effects are limited to within About one scale height of the 
boundary (scale heights are approximately 5 km at 90 km and 40 km at 250 km) . 


It has been found that the grid steps At and Az cannot be selected inde- 
pendently, due to numerical instabilities. Although the form of equations 
(6-18a, b) is somewhat more complicated, the stability criterion for the usual 
diffusion equation (equations (6-18a,b) with R and S set to sero) has proved 
to be a reliable indicator. This criterion, as given by Potter (1973), 
requires that the inequality 


At < i (m y 

At - 2 Q(z) 


(6-20) 


be satisfied throughout the domain of the calculation. Accordingly, the grid 
steps employed are Az = 5 km and At = 20 seconds. Inequality (6-20) is satis- 
fied at all altitude from 90 km to 250 km, although the limit is approached at 
the upper boundary. Considerations of requirements of the physical problem, 
computer storage and processing times, and numerical stability were all taken 
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into account to arrive at the grid step sizes and location of the upper 
boundary for the calculation. 

One of the early numerical tests performed with this set of equations was 
the use of the time-dependent ion equations of motion, to determine conditions 
for which the time- independent form is valid. The time-dependent ion equations 
of motion have the finite difference form (taking the 1-component as an 
example) 

(v il ) » +1 ‘ + At |\ B 1 c + "i (v i2 ) «" 

- - ‘’'.I*.”]] • «- 21 > 

With an initial condition of v, =0, it was found that for At > x. , where x. 

_ i in in 

* 1/v, is the mean time between ion-neutral collisions, numerical instabil- 
in 

ities rapidly developed. This is easily understood from the equations. 

Initially, the electric field term is the only nonzero term in both ion and 

neutral equations. If it accelerates the ions for a time At > , the ions 

attain an unrealistically large velocity. Neutrals respond slowly to ion 

collisions because of much larger number densities, so the neutral velocity 

remains close to zero. In the next time step and thereafter, the collision 

term in (6-21) dominates, with the result that the sign of the velocity 

alternates with each succeeding time step, and grows rapidly in magnitude. 

For At < t. , this does not occur, and reasonable results can be obtained, 
in 

For a constant electric field, it was found that the ion velocity reaches 
essentially equilibrium values in 6 to 8 collision times (t^ ) • Since x^ n 
varies from 3 x 10 seconds at 90 km to 1 second at 250 km (for models used 
here) , a wide range of equilibrium times corresponds to the altitude range 
under study. The difficulties in using equation (6-21) under such circum- 
stances are apparent. However, the important point of these results is not the 
range, but the fact that all times are short compared to the time scales of 
interest here. Since the Chatanika incoherent scatter radar* observations are 
averaged over about 6 minutes, that is the time scale for change of the elec- 
tric field. Hence, use of the time- independent ion momentum equations introduces 
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essentially no error into the calculation; and the ion velocity field can be 
considered to be always in equilibrium with the electric field and neutral 
velocity field. This results in a considerable reduction of the computational 
effort. 

6.4 MODEL CALCULATIONS 

Before performing an extended calculation and attempting to interpret the 
results, it is helpful to assess the importance of the various factors involved. 
While order- of -magnitude estimates are useful in judging the immediate effect 
of a variable on the system, total effects over extended periods of time are 
difficult to determine by this means. This is one purpose of the model calcu- 
lations to be presented here. By adopting a simple electric field model and a 
static ion density model, the temporal response, or rise time, of the ion- 
neutral wind system can be determined. Effects of coriolis and viscous forces 
can be ascertained by selectively femoving these terms from the equations and 
re-solving the problem. Once these effects are known, a better interpretation 
of the results for the physically realistic problem can be formulated. 

A model calculation of this nature has been performed by Fedder and Banks 
(1972). They consider a polar cap case in which magnetic field lines are 
vertical; and they retain the time-dependent form for the ion equations. 
Otherwise, their assumptions, approximations and resulting equations corre- 
spond to those above, with one important exception: they neglect the coriolis 

term. Effects of the coriolis force have been found to be significant by 
Heaps (1972) in calculations of high- latitude thermospheric winds over extended 
times; and it is expected that they should be significant here also. Since 
the analysis of Fedder and Banks (1972) treated the effects of collisional 
coupling and viscosity in some detail, the model calculations here are primar- 
ily to supplement those results by including coriolis effects. To accomplish 
this most clearly, the quasi-step function electric field used by Fedder and 
Banks (1972) is employed here: 

E = -E(t)y (6-22a) 

where 
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where 

E(t) =0, t < 0 

E(t) 20 [1 - cos(irt/1800) ] mV m- 1 , 

0 t _< 1800 seconds 

E(t) = 40 mV m \ t > 1800 seconds 


(6-22b) 


For these calculations, geomagnetic field geometry corresponding to 
Chatanika, Alaska is used. Differences in results between this geometry and 
the polar cap geometry of Fedder and Banks (1972) are small except for the 
vertical component of ion velocity, which is not considered in detail here. A 
representative ionization profile is used in conjunction with the ion compo- 
sition and neutral atmospheric models previously discussed to compute required 
collision frequencies. The ionization-density profile and resulting collision- 
frequency and gyrofrequency profiles are shown in Figure 6-6. 

Calculations for the electric field model in equations (6-22a, b) have 
been carried out' for a period of 24 hours, both with and without the coriolis 
force. The ion and neutral horizontal velocity components (geographic coordi- 
nates) thus obtained are shown in Figure 6-7 as functions of time for altitudes 
of 115 km, 125 km, and 150 km. Velocities computed without the coriolis force 
in Figures 6-7 a, c are consistent with the results of Fedder and Banks (1972); 
the more rapid response of their neutral velocities to the ion velocity driving 
force is due to their use of larger collision frequencies. 


In this approximation, it is seen that neutral velocity x-components 
(Hall drift direction) approach the ion- flow velocity over time periods of ten 
hours and more, the longer times corresponding to lower altitudes. This 
altitude variation of response times is understood in terms of neutral densi- 
ties increasing exponentially with decreasing altitude, in effect reducing the 
degree of ionization, even for little change in ionization number density. At 
low altitudes, the neutral atmosphere thus requires a larger total transfer of 
momentum to achieve a given velocity than at higher altitudes. Since the rate 
of momentum transfer varies less with altitude than the neutral number density, 
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Figure 6-7. VELOCITY COMPONENTS AS FUNCTIONS OF TIME AT ALTITUDES 
115 KM, 125 KM, AND 150 KM FROM MODEL CALCULATIONS: 
(a) x-COMPONENTS FOR CORIOLIS FORCE OMITTED; (b) x- 
COMPONENTS FOR CORIOLIS FORCE INCLUDED 
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Figure 6-7. VELOCITY COMPONENTS AS FUNCTIONS OF TIME AT ALTITUDES 
115 KM, 125 KM, AND 150 KM FROM MODEL CALCULATIONS: 
(c) y-COMPONENTS FOR CORIOLIS FORCE OMITTED; (d) y- 
COMPONENTS FOR CORIOLIS FORCE INCLUDED 





more time is required for the neutral velocity at lower altitudes to come into 
equilibrium with the ion (collisional) driving force. From the opposite point 
of view, the effect of the neutrals on ions is clearly evident in the response 
of the ions to the electric field. At 115 km, this response time is very 
large and is closely coupled with the neutral response time because of the 
close collisional coupling of the ions to the neutrals. 

The y-components are more difficult to understand immediately because 
they involve a balance between the geomagnetic field and collisional coupling 
with neutrals in controlling ion motions. However, upon reference to equations 
(6-14a, b) , it is seen that because of ion-neutral collisions at lower alti- 
tudes, both components of both the electrical field and neutral velocity 
affect each component of ion velocity. This directional coupling is the 
origin of the relatively large ion velocities parallel to the electric field 
(Pedersen drift) in the early stages at 115 km and 125 km. The neutral veloc- 
ities again respond to the collisional coupling and approach the ion velocity. 
However, because of the changing neutral velocity components the ion velocity 
y-component begins decreasing and becomes smaller than the neutral velocity. 

At this time the inertia of neutral gas results in a "coasting" motion, now 
becoming a driving force for the ions in maintaining their motion. At the 
same time neutral momentum is being diffused in altitude through viscous 
effects, as discussed in detail by Fedder and Banks (1972). 

These results may now be compared with those of Figure 6- 7b, d, for which 
the coriolis term is included. The effects on the neutral velocities are 
profound, and extend to the ion velocities at 125 km and 115 km where the 
collisional coupling of ions to neutrals is important. The magnitudes of the 
neutral velocity x-components are reduced to about one third of the equilibrium 
values being approached when coroilis forces are neglected, with the sign 
changing at lower altitudes. In the y-components of neutral velocity, the 
differences at different altitudes are particularly pronounced. Since the 
coriolis force basically couples the directions, the oscillatory behavior is a 
manifestation of a rotating velocity vector. The different effects at differ- 
ent altitudes are due to the variation of response times of the neutral velo- 
city field with altitude. 
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For example, at 150 km, where the response time is somewhat smaller than 
the 13 hour period of the coriolis term, the x-component of neutral velocity 
begins to grow rapidly, as in Figure 6-6a. Through coriolis coupling, however, 
part of this motion is diverted into the y-component. When a final steady- 
state is achieved the y-component is much larger than for the previous case 
and the x-component is much smaller. This is not a simple case of just 
rotating the velocity vector however, because the magnitude of the velocity 
vector is smaller when the coriolis force is included. The neutral response 
to collisions with ions is not rapid enough for the neutral velocity to ap- 
proach the ion velocity before coriolis effects significantly alter the 
direction of the former. 

At 115 km, it appears from Figures 6-7b, d that the neutral velocity 
field will oscillate due to the coriolis force for a considerable period (of 
the order of days). This appears to be due to the fact that the response time 
of the neutrals to ion drag exceeds the coriolis period (see Figures 6-7 a,b). 

An immediate consequence of these results is that large ion-neutral 
velocity differences are maintained at all altitudes, even when an equilibrium 
condition is reached. This is in contrast to the case when coriolis forces 
are ignored, in which the equilibrium condition is identical ion and neutral 
velocities. As a result. Joule heating, which varies as the square of the 
ion-neutral velocity difference (see equations (4-52) (4-54)) should remain 
relatively large under equilibrium conditions when the coriolis force is taken 
into account. This is consistent with the conclusions of Heaps (1974), who 
examined the effects of the coriolis force on Joule heating in a single-fluid- 
approximation calculation. 

Effects of coriolis forces are now examined explicitly in terms of alti- 
tude structure. Calculated altitude profiles of the velocity fields are shown 
for t - 5 hours and t - 24 hours for no coriolis force in Figures 6-8a, c and 
with the coriolis force included in Figures 6-8b,d. Figure 6-8a shows the 
extent to which the equilibrium is approached in 5 hours and Figure 6-8c shows 
virtually the final state, when no coriolis coupling of directions is consi- 
dered. At high altitudes the intuitive picture is fulfilled; below 125 km the 
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Figure 6-8. ALTITUDE STRUCTURE OF VELOCITY COMPONENTS FROM MODEL CALCULA- 
TIONS: (c) FOR T = 24 HOURS WITH CORIOLIS FORCE OMITTED; (d) 

FOR T = 24 HOURS WITH CORIOLIS FORCE INCLUDED 
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complexities introduced by the transition from ion-neutral collisional cbuplip^ 
to ion magnetic field coupling is evident, but with a similar end result: ion 

and neutral velocities are almost the same at all altitudes. 

In contrast, altitude plots in Figure 6-8b, d, for which the coriolis 
force is included, show a very different altitude structure for ions and 
neutrals at all altitudes. Effects on the ion velocity are seen to be small 
except at low altitudes, as expected. At high altitudes, the profile in 
Figure 6-8d represents virtually an equilibrium configuration, while at low 
altitudes the structure can be expected to change somewhat, as discussed 
previously. 

Similar calculations have been carried out with the viscous terms omitted. 
Results are not substantially different from those in which viscosity is 
included, however, so those results are not shown. The primary effect of 
viscosity is to diminish maximum values of the velocities and smooth out the 
altitude variations to a small extent. No significant differences in altitude 
structure of the velocity fields are seen. This is not really unexpected, due 
to the nature of viscous forces. It may be noted that Fedder and Banks (1972) 
found viscosity to be effective in transporting horizontal momentum through 
about one scale height. Their calculation was particularly designed to test 
this aspect by beginning with a strong vertical gradient at 155 km and no 
other driving force in the neutral wind equation. Because the collisional 
driving forces maintain the gradients in the present calculations, viscous 
forces are not seen so prominently. Under more realistic conditions, effects 
could well be more pronounced. Nevertheless, viscosity is not expected to 
alter the basic altitude configuration of the velocity fields to any great 
extent. 

The inverse situation is now examined; that is, the decay of the velocity 
fields from an initial configuration is calculated for the case in which the 
electric field is set to zero. It is noted that this should correspond to a 
very unlikely condition, because even in the absence of an external electric 
field, a dynamo electric field would be generated by the E region wind system. 
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However, it was seen in Figure 6-3a that electric fields are observed to "turn 
off" rather abruptly; so this example may not be so unreasonable. In any 
event the purpose here is simply to gain some idea of the inertial properties 
of the wind system, primarily to see how important the initial conditions are 
to the calculated winds some hours later. 

The initial wind configuration is that shown in Figure 6-7d, which is the 
calculated velocity field after 24 hours under the influence of the electric 
field in equation (6-22) . The full set of equations (6-14) and (6-15) is 
employed, as it was to obtain the initial configuration. Temporal variations 
for the velocity components at 115 km, 125 km, and 150 km are shown for 10 
hours in Figures 6-9a, b. Altitude structure of the velocity fields at 3 
hours and 10 hours after decay begins are shown in Figures 6-9c, d. From the 
detailed numerical results it is found that the coriolis term dominates the 
motion, particularly at lower altitudes. This is evident in Figures 6-9a, b. 
Viscosity is found to be more important in the decay mode than when the veloc- 
ity field is building up. 

Perhaps the most important feature of these results, however, is the 
inertia in the 110 km to 120 km region. After 10 hours with no momentum 
input, velocities there are reduced by about one half, while at 200 km the 
velocities have all but vanished. This indicates that initial conditions for 
the calculations with observed electric fields may be significant for an 
extended period after the calculation begins. It also shows that small elec- 
tric fields for the initial 2 or 3 hours of a calculation do not guarantee 
that initial neutral velocities will be small or will become small before the 
velocity fields are significantly affected by the larger electric fields of 
interest. This must be considered in the comparison of calculated and observed 
neutral velocities. 

Results of interest from these model calculations may be summarized as 
follows. In the development of the velocity fields through collisional coup- 
ling alone, 8 to 10 hours are required at high altitudes (>150 km) for the 
neutral velocity to approach the ion velocity, while at low altitudes (<120 km) 
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Figure 6-9. VELOCITY FIELD DECAY (NO ELECTRIC FIELD) FROM INITIAL 
CONFIGURATION GIVEN IN FIGURE 6-8d: (a) x-COMPONENTS 

AS FUNCTIONS OF TIME FOR ALTITUDES 115 KM, 125 KM, AND 
150 KM; (b) SAME AS (a) FOR y-COMPONENTS; 
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Figure 6-9. VELOCITY FIELD DECAY (NO ELECTRIC FIELD) FROM INITIAL 
CONFIGURATION GIVEN IN FIGURE 6-8d: (c) ALTITUDE 

PROFILES AT T = 3 HOURS: (d) ALTITUDE PROFILES AT 
T = 10 HOURS 


more than 24 hours are required. The coriolis force plays a significant 
role in both development and decay of the velocity fields. It reduces the 
magnitude of the maximum neutral velocity attained, and maintains a larger 
ion-neutral velocity difference in a steady state. The coriolis term causes 
an increase in the time required for an equilibrium velocity field to be 
achieved, particularly at low altitudes. In the decay mode from large initital 
velocities in the absence of driving forces, the coriolis term is dominant. 
Finally, large initial velocities can persist without driving forces in excess 
of 10 hours at low altitudes; hence effects of initial conditions can be 
expected for several hours into a calculation at these altitudes. 

6.5 VELOCITY CALCULATIONS FOR 15 MAY 1974 

A series of calculations, based on observations made by the incoherent 
scatter radar facility at Chatanika, Alaska on 15 May 1974, are now presented. 
Since successive calculations are motivated by preceding ones, each set is 
examined and discussed in turn before presenting results of the next. This 
facilitates clarity and allows interpretation to be developed in an orderly 
manner. 

6.5.1 Calculated Ion and Neutral Velocities as Functions of Altitude and Time 

The first set of calculations is the simultaneous numerical solution of 
the ion and neutral momentum equations (6-14) and (6-15) . The ion mean gyro- 
frequency ok and ion-neutral mean effective collision frequency vary with 

altitude but are constant in time under approximations above. From the inco- 
herent scatter observations, electric field components are prescribed as 
functions of time (shown in Figure 6-3a) and assumed uniform throughout the 
region of interest. Remaining parameters in equations (6-14) , speed of light 
c and geomagnetic field strength B, are constants. Neutral wind components 
are computed from equation (6-15) , according to the numerical scheme discussed 
in subsection 6.3.2. Neutral mass density and viscosity coefficient are 
functions of altitude but constant in time. Geographic latitude A and angular 
velocity are constants. The neutral-ion mean effective collision frequency 
varies linearly with ionization density and, hence, varies in both altitude and 
time. To determine these collision frequencies a collision frequency profile 


6-36 



-3 

normalized to 1 ion cm (which varies only with altitude) is multiplied by 
the ionization density profile appropriate to each time. A new electric field 
vector and ionization density profile are provided by the incoherent scatter 
radar at approximately 6 minute intervals; linear interpolation is used to 
determine values at intermediate time points. 

Calculated values of the ion and neutral velocity components in geographic 
coordinates are shown as functions of time in Figure 6-10 for altitudes 110 
km, 125 km, and 150 km. Altitude 110 km is shown in particular because it 
corresponds to the center of the first range gate, data for which is used to 
determine the neutral winds in Figure 6-3b. Several interesting features may 
be noted in these figures. First, the changing relations of the velocity 
components at different altitudes indicates changing altitude structure in the 
velocity fields. This is displayed more clearly in Figures 6-lla, b, c, which 
show the altitude structure of the velocity fields at UT = 1506, 1714, and 
2234 hours. Although these times were selected to display some of the extreme 
velocities, the form of the altitude structure is not atypical, as reference 
to Figures 6-8 and 6-9 will indicate. Vertical gradients in the horizontal 
velocity components can be quite large, particularly at altitudes below 130 
km. It is also of interest to observe in Figures 6-10 and 6-11 that the 
largest velocities for both ions and neutrals do not necessarily occur at the 
higher altitudes, but may also occur in the dynamic transition region, approxi- 
mately 115 km to 130 km for this calculations. 

A second significant feature of Figure 6-10 is the slow and almost sinu- 
soidal variation of the neutral velocity components, in contrast with the 
highly structured temporal variation of the ion-velocity components. This 
lattex structure clearly follows directly from similar structure in the 
electric field components, although the correspondence is not exact. Because 
of large numerical superiority, and therefore greater total inertia, the 
neutral gas averages over the short period fluctuations in the ion motions. 

This response is not unexpected, since it is consistent with the relatively 
long response times seen previously in the model calculations. 
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Figure 6-10. CALCULATED VELOCITIES AS FUNCTIONS OF TIME FOR ALTITUDES 110 KM, 125 KM, AND 150 KM, BASED ON 
OBSERVED ELECTRIC FIELDS AND ELECTRON DENSITIES: (a) x-COMPONENTS 
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Figure 6-10. CALCULATED VELOCITIES AS FUNCTIONS OF TIME FOR ALTITUDES 110 KM, 125 KM, AND 150 KM, BASED ON 
OBSERVED ELECTRIC FIELDS AND ELECTRON DENSITIES: (b) y-COMPONENTS 


VELOCITY (M/S) 




TIME (UT) 


Figure 6-10. CALCULATED VELOCITIES AS FUNCTIONS OF TIME FOR ALTITUDES 
110 KM, 125 KM, AND 150 KM, BASED ON OBSERVED ELECTRIC 
FIELDS AND ELECTRON DENSITIES: (c) x-COMPONENT, ION 

ONLY; (d) x-COMPONENT, NEUTRAL ONLY 
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Figure 6-10. CALCULATED VELOCITIES AS FUNCTIONS OF TIME FOR ALTITUDES 
110 KM, 125 KM, AND 150 KM, BASED ON OBSERVED ELECTRIC 
FIELDS AND ELECTRON DENSITIES; (e) y-COMPONENT, ION 
ONLY; (f) y-COMPONENT, NEUTRAL ONLY 
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Figure 6-11. EXAMPLES OF ALTITUDE STRUCTURE FOR VELOCITIES OF FIGURE 
6-10: (b) AT UT = 1714 HOURS 
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6-10: (c) AT UT = 2234 HOURS 





A major difference between this calculation and the model calculations, 
however, in addition to the large electric field fluctuations in this case, is 
large fluctuations in ionization density. This is illustrated in Figure 6-12 
where the temporal variations of ionization density at 100 km, 110 km, 120 km, 
and 130 km are plotted. Large fluctuations at a given altitude in short times 
are frequent during the times of greatest electric field fluctuation. It 
might be suspected that resulting fluctuations in the neutral-ion collision 
frequency v ni could induce similar variations in the neutral velocities. Upon 
examination of equations (6-15) however, it is recognized that an order of 
magnitude increase in the collision frequency has about the same effect as 
a similar increase in the ion velocity (electric field) : an order of magnitude 

increase in the collision term. Over short times, neither produces a large 
variation in the neutral velocity, again because 1^/N^ <<: 1 at all altitudes. 

Considerable attention is directed to the temporal variation in the 
neutral velocity components because the neutral velocity components inferred 
from the incoherent scatter radar observations, shown in Figure 6- 3b, have 
large fluctuations in short times. Thus, the temporal structure of the observed 
neutral winds more closely resembles that of the calculated ion velocities 
than that of the calculated neutral velocities. Even without considering in 
detail the magnitudes of these velocities, the nature of this discrepancy 
indicates a problem. If the winds inferred from the radar observations are 
consistent with the assumptions made in the data anlaysis and are approximately 
correct (indicative of the actual neutral winds at that location) , then from 
the discussion above, ion-neutral collision coupling alone cannot account for 
them: some important mechanism has been omitted. On the other hand, to the 
extent that collisional coupling is important, Figure 6-11 shows that one of 
the important assumptions in the radar data analysis for neutral wind determin- 
ation is called into question. The data analysis procedure assumes a neutral 
wind which is uniform throughout the first range gate, and it is such a wind 
that is determined. Figure 6-10 shows that considerable altitude structure 
does exist in this region, unless collisional coupling is unimportant, in 
which case the data analysis procedure is invalid anyway (Brekke, et al. 

1973). 
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Figure 6-12. ELECTRON NUMBER DENSITIES AS FUNCTIONS OF TIME FOR ALTITUDES 100. KM, 110 KM, 120 KM, 130 KM 


With respect to the effect of altitude structure on the analysis of radaF 
data, two possibilities can be seen. The radar is assumed to average the 
observation velocity over altitude as in equation (6-4) . Weighting the 
average are the electron (ion) density ^(k) and the radar weighting function 
W(z). In Figure 6-12 it was seen that electron densities can change rapidly* 
in time. Examples of how rapidly the altitude structure of N e (z) can change 
are shown by the altitude profiles in Figure 6-13. It is evident that even if 
the velocity structure remains constant, variation of the electron density 
profile can cause a shift in the altitudes most strongly weighted. From 
Figure 6-11, it is apparent that small shifts in altitude can correspond to 
large changes in velocity. Thus, the scatter radar observations of apparent 
temporal variations in ion velocity can result solely from variations in the 
electron density profiles, due to large vertical gradients in the velocity 
fields. Brekke, et al. (1974a) note this problem; although the mode of oper- 
ation for the present experiment improves the situation somewhat, due to 
better time resolution, it nevertheless remains a difficulty. 

A second possible effect is due to variation in the ion velocity altitude 
structure itself. Since the ion velocities are immediately affected by elec- 
tric field variations, over a period of a few minutes the ion velocity altitude 
structure can change considerably, while that of the neutral velocity is 
virtually unaltered. If an effective altitude of observation is assumed to be 
defined as that altitude at which the weighted average velocity and the actual 
velocity agree, it can be seen that this effective altitude can shift by a 
change of ion velocity altitude structure, even for fixed electron density 
profiles. 

It appears possible that some of these changes may be adequately compen- 
sated by the data analysis procedure, equations (6-6) and (6-7). For example, 
effects of variations in electron density profiles are taken into account to 
some extent by different values of K-coef f icients which essentially weight 
different altitudes (see Figure 3 of Brekke, et al. (1973); note, however, 
that their and labels are reversed). Also, since both the electric 

field and ion velocity appear in equation (6-7), changes in the one may 
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compensate changes in the other, with little effect on the neutral velocity 
determined by that equation. 

Nevertheless, the question remains, how are the calculated velocities of 
Figure 6-9 to be compared with the results of observation in Figure 6-3? 
Evidently, the answer is not to choose a fixed altitude and compare the calcu- 
lated velocities as functions of time at that altitude with observed velocities. 
A promising approach, however, appears to lie in averaging the calculated 
velocities in altitude after the manner of the radar; that is, numerically 
simulate the radar averaging. This has the advantage of providing a common 
form for the observed and calculated velocities, making a valid comparison 
possible. It also offers prospects for resolving some of the difficulties 
noted above. 

6.5.2 Weighted Altitude Averages of Calculated Velocities 

Computation of weighted altitude averages of the calculated velocities 
presents no difficulties. Equation (6-4) defines the manner in which the 
average is weighted, and equations of this type are used to compute weighted 
altitude averages for both ion and neutral velocities. Numerically, the 
procedure is the same as that used in computing the K-coef f icients by equations 
(6-6a through c) . For ion velocities this result should simulate the inco- 
herent radar observation, to the extent that the calculated altitude structure 
corresponds to the actual altitude structure. For neutral velocities, the 
result represents what the radar would obtain if it were sensitive to the 
neutral gas in the way it .is to ions. 

This calculation also provides an opportunity for assessing the accuracy 
of the data analysis procedure for determining the neutral wind. On the one 
hand, the presumably accurate weighted altitude average of the neutral wind is 
obtained by direct numerical integration. On the other hand, the same is done 
for the ions, which should accurately represent the radar result. By treating 
this result as an observation, that is, by using it as an input to the data 
analysis procedure (equations (6-6) and (6-7)), neutral^ winds are obtained 
reflecting all systematic errors and inconsistencies inherent in the procedure. 
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Direct comparison with the accurate results then provides a test of the 
accuracy of the procedure under realistic conditions. It is emphasized that 
this is not a comparison of theory with observation; rather it is a test of 
the data analysis procedure (different theory) to see how well it represents 
what it is supposed to represent. 

Before presenting results of these calculations, it is helpful to intro- 
duce some notational guides. Used consistently, these will help reduce con- 
fusion as different velocities are examined and discussed through the remainder 
of this section. Calculated velocities will be denoted by lower case v, with 
subscript i or n indicating ion or neutral as before. A bar over the velocity, 
v, will indicate a weighted altitude average. Observed^velocities , or those 
obtained from observations, will be denoted by capital V, where the bar indi- 
cates that it is a weighted altitude average. Within this framework, veloc- 
ities will be distinguished by whatever means is appropriate, consistent with 
previous notation. 

In terms of this notation, the weighted averages of the ion and neutral 
velocities^obtained directly from equation (6-4) (<jr its analog) are denoted 
by v^ and v^. The neutral velocity obtained f^om v^ by means of equation 
(6-7), as discussed previously, is denoted by v R f . Horizontal components 
(geographic coordinates) of v^ and v^’ are presented in Figure 6-14. 

— 

Examination of v , with reference to v at constant altitudes in Figure 
n n 

6-10, shows that the process of taking weighted altitude averages introduces 
considerable structure into the time variation of the velocity. It can be 
inferred that fluctuation of the electron density profile is the cause, since 
W(z) does not vary in time and the neutral wind is seen in Figure 6-10 to vary 
relatively slowly at all altitudes. If the general trend of v^ x is taken as a 
reference in Figure 6-14, it appears that perturbations of up to 50 percent 
can result from such shifts in the altitude weighting. The absence of similar 
large perturbations in the y-component appears to be due to the different 
altitude structure and smaller magnitude of this component during periods of 
large electron density fluctuations (see Figures 6-11 and 6-12). Finally, 
even though the weighting function W(z) is centered at 109 km, reference to 
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Figure 6-14. COMPARISON OF WEIGHTED ALTITUDE AVERAGES OF CALCULATED NEUTRAL VELOCITY COMPONENTS OBTAINED 
DIRECTLY BY NUMERICAL INTEGRATION (v ) AND OBTAINED INDIRECTLY THROUGH THE DATA ANALYSIS 
PROCEDURES (v ' ) " 



Figure 6-10 shows that the magnitudes of and correspond to altitudes 
several kilometers higher, the exact altitudes varying with the electron 
density profiles. 

The comparison of v^ with v^' in Figure 6-14 demonstrates that the^data 
analysis procedure consistently underestimates the actual magnitude of v^ 

(and therefore V^) by amounts ranging up to 40 percent. This appears to be a 
systematic effect associated with the altitude structure -of the neutral wind, 
since the more transient effects of electron density fluctuations are evident 

in both v and v ' . The notable feature of this comparison, however, is how 

n n + + 

well the temporal variations of v^ are reproduced by v^ 1 . With allowance for 
the systematic bias in magnitude, it is clear that the data analysis procedure 
does provide a reasonable representation of the weighted altitude average of 
the neutral velocity. However, it is equally clear that the time variation of 
this velocity is not the same as a straight altitude average of the neutral 
wind, because the weighting due to N^fz) has been seen to introduce some of 
this variation. 

With regard to the magnitude of the velocities determined by equation 
(6-7) , a brief digression is instructive. To assure that the data analysis 
procedure was being executed exactly as prescribed by Brekke, et al. (1973) , a 
portion of the data for 15 May 1974 was processed through the computer program 
used for that study by J.. R. Doupnik and P. M. Banks, and used to verify 
comparable calculations in this investigation. Proper use of the procedure 
has been comfirmed. However, in the course of this verification it has been 
seen how sensitive the magnitude of the velocity is to model values entering 
into the collision frequency v^. To compute ion-neutral collision frequencies, 
Brekke, et al. (1973, 1974a) have used the expression 

v. = 7.5 x 10 -10 N , (6-23) 

m n 

with the total neutral number density obtained from a neutral atmospheric 
model given by Banks and Kockarts (1973), corresponding to an exospheric 
temperature of 1000°K. Collision frequencies used here are computed according 
to equation (4-14c) from the individual species collision frequencies given in 
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subsection 5. 3, 3. They also depend on the ion composition model and neutral 
atmosphere model discussed in subsection 6, 3 When these different approaches 
are taken and the final values are computed and compared , the collision fre- 
quencies, used by Brekke, et al (1973, 1974) are larger than those used here 
by a factor of about 1.5 at 90 km, increasing with altitude to a factor of 2 
at 160 km. 

Neutral winds calculated according to equations (6-6) and (6-7) by 
Doupnik and^Banks are denoted by V 1 , while those calculated in this study are 
denoted by V. This notation indicates that both are computed from_^the obser- 
vational data. Horizontal components (geographic coordinates) of V' and V 
for 0600 hours UT to 1800 hours UT on 15 May 1974 are shown in Figure 6-15. 

It is emphasized that both of these represent "observed' 1 neutral wind veloc- 
ities, in the sense that they are obtained from observational data. The 
procedure for their calculation is the same; only the collision frequencies 
differ significantly. This comparison shows basically the same temporal 
variation, but magnitudes are frequently substantially different. These 
differences are not systematic in an obvious manner. 

The purpose of this comparison is to emphasize the sensitivity of the 

neutral winds derived from observation to the collision frequency model employed 

to deduce them. Because of uncertainties involved in the calculation of ion- 

* 

neutral collision frequencies, from the pertinent interactions to the numbers 
of each species of ion and neutral particle actually present, it is difficult 
to support one set of collision frequencies against another. However, to the 
extent that a degree of arbitrariness is involved, its effects on the final 
results - neutral wind velocities - should be recognized and acknowledged. 

In a study of this nature, such considerations* do not undermine the 
comparison between observation and theory, so long as the same collision 
frequencies are used consistently in both the data analysis and the theoretical 
calculations. There is at least a rough scaling with collision frequency of 
neutral velocities obtained from both theory and data analysis, as indicated 
by Figure 6-14. One calculation which could be adversely affected, however. 
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Figure 6-15. NEUTRAL VELOCITIES DETERMINED FROM INCOHERENT SCATTER RADAR 
OBSERVATIONS USING THE COLLISION FREQUENCIES OF THIS STUDY 
(V„ - SAME VELOCITIES AS FIGURE 6-3b) ) AND USING THE 
COLLISION FREQUENCIES OF BREKKE ET AL. (1973) (V '): 

(a) x-COMPONENTS; (b) y-COMPONENTS 
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is that of Joule heating, which depends on (V^ - V^) . Since is determined 
directly by observation, its accuracy is affected only by the experimental 
error. However, because is also subject to this additional bias in the 
data analysis, a Joule heating calculation is susceptible to similar error. 

This is merely noted for future studies because no such calculation is carried 
out here. 

With that digression concluded, the primary comparison can be made between 
the theoretical weighted altitude average of the neutral velocity, v^, and 
that determined from observation, V^. These are presented in Figure 6-16. 

The initial impression is one of substantial disagreement. First, despite the 
temporal structure introduced into v^ by taking the weighted altitude average, 
the temporal variation in is larger in both frequency and magnitude of 
fluctuations. From the comparison in Figure 6-14, it is clear that is a 
reasonably accurate reflection of the weighted altitude average of neutral 
velocity in terms of the observational data. However, Brekke , et al. (1974a) 
note that the data represent time averages - for this experiment over a period 
of 6 minutes. If velocities or electron densities are changing substantially 
during this time, ^he results will be inaccurate. It is evident that this 
happening when large differences are seen between consecutive measurements; 
and Figures 6-3a and 6-12a show this to occur frequently in this experiment. 

So these large rapid fluctuations must be taken as a qualitative indication of 
variability, rather than a quantitative measure of reality. The more stable 
quantities are a more accurate measure of the physical variables. 

it — 

Even from this point of view, the discrepancy between v^ and remains 

substantial. Two characteristics are apparent. First, the general magnitudes 

of V and V are considerably smaller than the broad maximums attained by 
_ nx _ ny ± it 

v and v . Secondly, the general directions of V and v bear little coher- 

nx ny . J n n 

ent relation to one another. Consequently it can be concluded that the theoret- 
ical basis for the calculated results is inadequate for describing the iono- 
spheric conditions actually present when the observations were made. However, 
since the velocities calculated theoretically are somewhat larger in magnitude 
than those observed, the collisional coupling is evidently important. It 
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Figure 6-16. COMPARISON OF OBSERVED NEUTRAL VELOCITIES (V) AND WEIGHTED ALTITUDE AVERAGES OF THE CALCU 
LATED NEUTRAL VELOCITIES (v n ) (a) . x-COMPONENTS; 



appears that other mechanisms of comparable importance are operating in at 
least partial opposition to the collisional force so that the build-up of 
neutral velocities is inhibited. 

It is helpful to recall briefly the basis of the theoretical calculation 
in order to see the implications of these results. The physical model is one 
of large electric fields of magnetospheric origin causing ions to move with 
large flow velocities. Some of this momentum is transferred to the neutral 
gas through collisions, causing it to move. The neutral gas motion is then 
modified somewhat by viscous and coriolis forces, these modifications being 
coupled back to the ions by collisions. From the comparisons in Figure 
6-16, these processes are important, but insufficient to describe the actual 
physical behavior of the sytem. However, more information is now available 
than initially. Not only are velocities due to collisional coupling known, 
but the size of the discrepancy between these velocities and the observation 
velocities is known. This information can now be used to estimate the terms 
omitted form the theoretical analysis, that is, those terms required to resolve 
the discrepancy. 

6.5.3 Estimates of Pressure Gradients 

In order to carry out a calculation of this nature the original assump- 
tions in subsection 6.3.1 must be reviewed. Those involving the ion momentum 
equation are on a sound physical basis and appear to require no modification. 
However, in the neutral momentum equation, the neglect of horizontal gradients 
was recognized as a mathematical expendient at the outset; this must now be 
reconsidered. The uncoupling of the neutral vertical component and use of 
static equilibrium model in its place remains an approximation, but probably 
not a major contributor to the discrepancies in the horizontal component. 

So, the primary requirement for more realistic horizontal neutral equations 
of motion appears to be the inclusion of horizontal gradients. 

Thus, it is now assumed that the neutral velocities determined from 
incoherent scatter radar observations obey the following equations of motion: 
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In order to make use of the observational data’, weighted altitude averages are 
required. This is accomplished formally by multiplying equations (6-24a, b) 
by N^(z) W(z) , integrating over the first range gate, and normalizing to 
(equation ( 6 - 6 d)). When this is done the results may be written 
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where v ni is the weighted altitude average of v n ^, after the manner of equation 
(6-4) . An analogous equation may be written for V . 


Except for the last term in this equation, all terms in which the integral 

symbol is not written explicitly must be considered approximations, since they 

do not follow directly from equation (6-4). The left side in particular 

requires that electron densities vary slowly during the interval in which the 

time derivatives are evaluated. Since successive radar observations are used, 

this is consistent with the similar requirement for accurate observational 

data. It must be noted that only V and V are known from observational 

nx ny 

data. So those terms remaining as integrals are to be regarded as unknown 
quantities. 
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To make a comparison with the calculated results* a similar operation 
must be performed on the original equations (6- 15a, b). This gives 
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Since v and v are known as functions of altitude, the integrals in these 
equations can be evaluated numerically without difficulty; that is, each term 
is a known quantity. 


The difference between equations (6-25) and (6-26a) is now taken; this 
results in 
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with a similar equation for the y-component. 
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Several of these terms are probably small and can be omitted. The first 

term on the right side is negligibly small if the altitude structures of 

and v are similar, or if both terms are small. This term for v has been 
nx nx 

found in numerical calculations to be always at least one to two orders of 
magnitude smaller than the collision terms at these altitudes. Moreover, it 
seems unlikely that has larger vertical velocity gradients, since the 
magnitude of V^ x is generally smaller and vertical gradients in v rx are them- 
selves large. Hence, it appears reasonable to omit this term. As another 
approximation, the other viscous terms (second and third terms on the right 
side) are assumed negligible in comparison with the nonlinear term (fifth tern 
on the right side). This seems reasonable from order of magnitude considera- 
tions, which can be specified more precisely than in Section IV because condi- 
tions are more restricted. The ratio of viscous to nonlinear terms is given 
by (order-of-magnitude only) 
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typical horizontal scale length exceeding 10 to 10 cm. 


With these approximations, equation (6-27) is simplified and rearranged 
so that unknown quantities appear on the left side while terms which can be 
evaluated are on the right. The result is 
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The corresponding equation for the y-component ia 
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The left side is seen to include both the pressure gradient and nonlinear 
terms. Their relative contributions can be estimated, based on observed 
velocities, if it is assumed that the horizontal scale length L is approximately 
the same for both pressure and velocity variations. With that provision, the 
ratio of nonlinear to pressure gradient terms is (order-of-magnitude) 


($ • v)v 

n nx 

1_ 92 

p n 9x 


M 

V z 

n 

n 

k 

T 

n 


(6- 30a) 


From Figure 6-16, observation values of are of the order of 100 m 6 

With assumed values of m - 25 AMU and T - 400 °K, consistent with the neutral 

n n 

model employed, the ratio is 


• V)V 

1 8p 

p n 3x 


nx 


< 10 


-1 


(6-30b) 


Actually, velocities in Figure 6-16 sometimes exceed 100 m s ^ by factors of 2 
to 4; and, on the other hand, T^ very likely exceeds 400°K by substantial 
amounts, on occasion, due to local heating. For this reason, the nonlinear 
term is retained in equations (6-29a, b). However, from (6-30b) , it is clear 
that it is the pressure gradient which is generally being evaluated by these 
equations. 


Numerical evaluation of terms on the right side of equations (6-29a, b) 
gives the accelerations shown in Figure 6-17a, c. These are based on the 
calculation and observation values shown in previous figures, and at time 
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Figure 6-17. PRESSURE GRADIENT ACCELERATIONS AS DETERMINED FROM 

EQUATIONS (6-29a, b): (a) x-COMPONENT; (b) y-COMPONENT 
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AVERAGED BETWEEN CONSECUTIVE OBSERVATION TIMES 
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AVERAGED BETWEEN CONSECUTIVE OBSERVATION TIMES 


points corresponding to the observations. The rather large fluctuations which 
frequently occur between consecutive points tend to obscure the general trend. 
Since they may result from inaccuracies in the data, due to rapidly changing 
conditions, and from contributions from the nonlinear term for the larger 
velocities, these fluctuations do not necessarily represent good values of 
the pressure gradient. To see the general magnitude of this term more clearly, 
these fluctuations are smoothed somewhat by taking simple averages between 
consecutive points. These results are shown in Figure 6-17b, d. 

Despite this averaging of results, considerable temporal structure remains. 
Nevertheless , trends in the magnitude and direction of the pressure gradient 
are now more apparent, particularly during the second electric field enhance- 
ment. In the x-component, the first two hours indicate no definite direction. 
During the first enhancement (0600 to 1100 hours UT; see Figure 6-3a) this 
component fluctuates considerably, but appears to be somewhat more negative 
than positive, indicating generally higher pressures (heat sources) to the 
north. Between enhancements, the pressure gradient is relatively small, and 
perhaps slightly more negative (northward) than positive. For the duration of 
the second electric field enhancement (1200 to 1800 hours UT) , the meridional 
pressure gradient is predominately positive, indicating heat sources primarily 
to the south. Thereafter, this component becomes small, fluctuating about 
zero for about two hours, and then remains small but primarily negative (heat 
sources to the north) for the remainder of the experiment. 

A similar examination of the y-component of the pressure gradient reveals 
both similarities and differences. The initial two hour period is again 
characterized by small magnitudes with indefinite direction. During the first 
electric field enhancement, the sign is more definite than for the x-component 
and the fluctuations smaller in amplitude. The predominately positive pressure 
gradient (heat source to the east) appears to become primarily negative 
shortly before 1000 hours UT and remains moderately small until about 1230 
hours. For the duration of the second enhancement this component is primarily 
negative (heat sources to the west) , although large rapid fluctuations occur 
between 1500 hours and 1800 hours UT. In contrast with the x-component. 
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however, the y-component does not become appreciably smaller in magnitude 
after termination of the second electric field enhancement. Rather it becomes 
more definite in its sign, still indicating heat sources to the west. This 
situation remains until about 2200 hours UT, when the sign becomes indefinite. 

Magnitudes of these accelerations are seen to be a few tens of centimeters 
per second per second during electric field enhancements and a few centimeters 
per second per second otherwise. This distinction is less clear for the y- 
component, however, than for the x-component. These magnitudes may be compared 
with pressure gradient accelerations used in global theoretical models of the 
neutral wind system. Pressure gradients in the models of Kohl and King (1967) 
and Blum and Harris (1975) are both based on one of the neutral model atmos- 
pheres of Jacchia. Amplitudes of these pressure gradients for E region alti- 
-2 

tudes are < 1 cm s . Thus, pressure gradients during disturbed conditions at 
high latitudes are typically at least an order of magnitude larger than those 
driving the global neutral wind models, which are due primarily to solar 
heating at low latitudes. It is therefore not surprising that observations of 
neutral winds at high latitudes under geomagnet ically disturbed conditions 
frequently disagree with these models (e.g. Stoff regen, 1972; Meriwether, 
et al. 1973; Brekke, et al. 1974a; Rothwell, et al. 1974). For these cond- 
itions no other calculated or experimentally derived estimates of horizontal 
pressure gradients are known to be available for comparison. Accelerations 
due to ion drag are of the same order as these due to pressure gradients. 

To interpret these results, it is important to note a significant differ- 
ence between ion drag and pressure gradients as driving forces for neutral 
winds. For ion drag the direction of the electric field as well as its distri- 
bution and magnitude must be considered. Pressure gradients, however, depend 
only on the distribution (geometry) and magnitude of pressure enhancements 
(primarily heat sources). Since these results involve only pressure gradients, 
interpretation requires an explanation of the implied distribution of heat 
sources. Local heating, for example Joule heating associated with ion-neutral 
velocity differences at the point of observation, need not be considered, 
since it would represent a local maximum for which the local gradient would 
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vanish. Rather it is heat sources in surrounding regions, which are not 
observed by the radar, that must be considered. 

Comparison of the electric fields in Figure 6- 3a and the electron number 
densities at constant altitude in Figure 6-12 shows that periods of enhanced 
electric field occur during the same periods as auroral E layers. This is not 
meant on a (time) point by point basis, because large electric field magnitudes 
are often anticorrelated with E region electron densities; 'but rather periods 
of large rapid fluctuations in electric field magnitude coincide with similar 
fluctuations in E region electron densities. Since a developed auroral E 
layer is a reliable indicator of auroral activity, it can be concluded from 
Figure 6-12 that the ionosphere over Chatanika was in the auroral oval, that 
is, in the midst of auroral activity, from shortly after 0600 hours to about 
1800 hours UT, with a possible exception during a brief period between 1200 
hours and 1300 hours UT. Further, auroral heating in the form of both energetic 
particle precipitation and Joule heating was occurring concurrently in the 
same locations. Based on statistical models of the auroral oval (Feldstein 
and Starkov, 1967), it can be assumed that before and after this period, all 
auroral activity in the vicinity of Chatanika was to the north. 

In this context the pressure gradients in Figures 6-17c, d can be examined. 
During the initial two hour period, no particular directionality is evident in 
the pressure gradient. After Chatanika enters the auroral, oval near 0600 
hours UT (2000 hours AST), the pressure gradient indicates that heating lies 
primarily to the north and east, which is consistent with the general picture 
of auroral activity being concentrated toward the local midnight and early 
morning sector (Akasofu, 1968). From 0900 hours to 1100 hours UT (1000 hours 
UT is AST midnight), no consistent directions are indicated; this can be 
viewed as a consequence of a somewhat symmetric distribution of heat sources 
about Chatanika. The period 1100 hours to 1300 hours is somewhat ambiguous; 
a more detailed discussion is deferred. During the second active period, from 
1300 hours to 1800 hours UT (0300 hours to 0800 hours AST) the pressure gradi- 
ents indicate predominant heating to the south and west. Again the heating 
toward the west is consistent with primary activity in the local midnight and 
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early morning sector of the auroral oval. The definite southward direction of 
the heat sources shows that Chatanika is well within the auroral oval. Soon 
after 1800 hours UT the electric field decreases abruptly to very small values, 
concurrent with the disappearance of the auroral E layer. Thus, Chatanika has 
left the auroral oval, or vice versa. 

After 1800 hours UT, local activity is at a low level. The y-component 
of the pressure gradient, however, remains relatively large in magnitude with 
a definite westward direction, decaying slowly over the next 4 hours. It 
appears that heating on the night side has built up and maintained pressure 
differentials which are affecting the dayside region. In contrast, the x- 
component drops rapidly to very small magnitudes by 1830 hours UT, and remains 
small and essentially undirected for about an hour. Then it becomes definitely 
negative, indicating heat sources to the north where the dayside auroral-oval- 
cusp region is located. 

Separation of spatial and temporal effects in observations at a single 
location is difficult, particularly when it involves a dynamic region like the 
auroral oval, which has a relatively fixed orientation with respect to the sun 
but can expand or contract due to magnetospheric electrodynamics. Nevertheless, 
some plausible inferences can be drawn. The period 1100 hours to 1300 hours 
is a good example* Between 1100 hours and 1200 hours UT the electric field 

(Figure 6-3a) has become quite small, while a small auroral E layer (Figure 6- 

* 

12) persists. Between about 1215 hours and 1245 hours UT the auroral E layer 
vanishes and a large electric field develops. Several possible explanations 
for these observations exist: (1) the auroral oval could have contracted so 

that Chatanika was to the south of its equatorward border; (2) Chatanika could 
have encountered an extended region of low activity (first a region of low 
electric field, then low auroral precipitation), while activity continued in 
surrounding regions ; this could occur by Chatanika rotating through a more or 
less fixed region of low activity, or the region could drift by Chatanika due 
to electrodynamic forces; or (3) activity could simply cease in time over a 
large spatial region due to source region temporal variations. 
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The first possibility does not seem to be supported by the available 
information. It could not have occurred except during the brief period when 
the auroral E layer vanished. Calculations of Jones and Rees (1973) show 
that an auroral E layer will decay to background levels in a minute or less 
after removal of the ionization source. This appears to happen near 1215 
hours UT, but not prior to that time. However, during the period of no auroral 
E layer, a large electric field develops. For all other times in these obser- 
vations, large electric fields occur in the immediate vicinity of auroral 
activity. Furthermore, during this same time a substantial pressure gradient 
develops directed first toward the south. This could not happen if Chatanika 
were to the south of the auroral oval at this time. For this same reason, the 
second explanation appears more likely than the third. It seems probable that 
Chatanika rotated into a slot between bands of activity oriented along the 
auroral oval. However, data from other locations would be required to substan- 
tiate this. 

Similar considerations can be applied to the time period after 1800 hours 
UT. Here the evidence supports the view that the position of Chatanika rela- 
tive to the auroral oval has shifted to the south. Electric fields and E 

/T*' 

layer number densities are consistent with this interpretation. Also the x- 
component of the pressure gradient favors this view. Heating to the south of 
Chatanika is indicated for several hours before. Shortly before 1730 hours 
UT, the pressure gradient begins to decline rapidly. This would be expected 
if the auroral oval were contracting northward so that pressures previously 
built up in the south were becoming balanced by pressures developing in the 
north. Such a balancing is more likely than a decay of the pressure in the 
south because these decay times are of the order of hours (Thomas and Ching, 
1969; Rees, 1971b). From Figure 6-17c, it appears that approximate balance is 
achieved about 1830 hours UT, lasting for approximately an hour. Continued 
heating to the north then prevails, shifting the dominant pressure to the 
north of Chatanika. 

While interpretations such as these cannot be proved without information 
covering an extended spatial region, they are consistent with the observations 
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at Chatanika and the general morphology of auroral activity. The results, as 
interpreted above, support the schematic model suggested by Stoffregen (1972) 
for explaining discrepancies between observed high latitude neutral winds and 
winds determined from the global theoretical wind models mentioned previously. 
Figure 6-18 shows his proposed schematic modification of global temperature 
distribution. This proposal is based on an analysis of 25 chemical release 
observations of neutral winds, restricted in local time to twilight hours, and 
for altitudes generally above 200 km in altitude. Nevertheless, the required 
heating would take place primarily in the E region and the general features of 
the diagram are consistent with the discussion of pressure gradients above. 

In particular the westward pressure gradient after 1800 hours UT (0800 hours 
AST) is satisfactorily explained, whereas the pressure gradient produced by 
solar heating alone should be eastward and of much smaller magnitude. 

A possible period of inconsistency is in the first two hours of the 
experiment, when auroral heating would perhaps be expected to produce an 
eastward directed pressure gradient. However, effects of this heating may be 
small, so that some form of balance is achieved with the small westward pres- 
sure gradient associated with solar heating and heating in the afternoon 
sector of the auroral oval. Also because of uncertainties involving initial 
conditions, the earliest calculated results are not necessarily reliable; this 
could skew the pressure gradient calculation during that time period. 

The broad features of the pressure gradient results have been discussed 
at some length because it is felt that they are more reliable than the small 
time scale features, due both to noise in the observational data, and to 
approximations required in deriving equations (6-29a, b) . However, the tem- 
poral structure evident in Figure 6-17 requires comment. This structure is 
almost certainly not all real. Some is due simply to experimental uncertain- 
ties in the observed ion velocities. Likewise, some is due to altitude struc- 
ture of the velocity fields and the weighted altitude averaging, as was seen 
to occur in v^. It should be reemphasized that the pressure gradients calcu- 
lated are weighted altitude averages. The altitude structure of these hori- 
zontal pressure gradients is not known. Conceivably, close to rapidly varying 
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WIND DIRECTIONS (NO ACTIVITY) 

( HIGH ACTIVITY ) 


Figure 6-18. SCHEMATIC MODIFICATION OF EXOSPHERIC TEMPERATURE DISTRI- 
BUTION PROPOSED BY STOFFREGEN (1972) TO INCORPORATE 
AURORAL HEAT SOURCES: (a) ORIGINAL DISTRIBUTION (JACCHIA 

1965); (b) STOFFREGEN'S (1972) MODIFIED DISTRIBUTION 





localized sources of Joule and particle heating, it could be complicated, but 
vertical motions would likely diminish such complexities rapidly. 

Despite these disclaimers, however, the temporal structure should not be 
dismissed entirely because there are physical reasons for expecting it to 
exist, based on discussions above. If the interpretations made previously are 
basically correct, Chatanika was more or less surrounded by particle and Joule 
heating while in the auroral oval. From the variability of the electron 
densities and electric fields, this heating is evidently highly variable in 
space and time. Thus, while heating may, for example, occur predominately to 
the south of Chatanika, sudden localized heating to the north may also occur. 
Depending on the strength and location of the northern heat source relative to 
the observation point, the southward pressure gradient could be diminished. or 
even reversed. Hence, rapid changes of neutral wind magnitudes and directions 
could be a consequence of the sudden development of heat sources in various 
nearby locations and the resultant fluctuations in pressure gradients. 

In this context back pressures must be considered. When a moving gas 
flows into a stationary gas (or one moving more slowly), density will increase. 
The localized pressure build up will result in a pressure gradient (back 
pressure) which opposes the motion of the fast gas but also transmits momentum 
to the stationary gas causing it to move. As a result, the original fast 
moving gas will be slowed and probably diverted. The precise role of back 
pressures in the auroral oval itself is difficult to assess. It is possible 
that much of the fluctuation in magnitude of the pressure gradients is associ- 
ated with back pressures. However, conservation of momentum permits a back 
pressure only to diminish the original pressure gradient, not to reverse its 
sign. Thus, those frequent fluctuations which involve a reversal of component 
direction require a different explanation, such as offered previously. 

Finally, some idea of the horizontal scale lengths (L) associated with 
these pressure gradients can be gained from order- of -magnitude estimates. 

Let A denote values of the acceleration due to the weighted altitude average 
of the pressure gradient, as determined from equations (6-29a, b). The order- 
of-magnitude is then estimated from 
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When this relation is solved for L, the result is 
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Based on T = 500 °K and m = 25 AMU, values for L are given in Table 6-3. 

Thus, scale lengths corresponding to pressure gradients in Figure 6-17 are 
seen to range generally from a few hundred to a few thousand kilometers, which 
are not unreasonable. 


Table 6-3. HORIZONTAL SCALES ASSOCIATED WITH PRESSURE 
GRADIENT ESTIMATES 


1- 2 
P n v p (cm s ) 

1 

10 

100 


L (km) 

2 x 10 4 

2 x 10 3 

2 x 10 2 . 



6.6 DISCUSSION AND CONCLUSIONS 

In recent years several analyses of neutral wind observational data have 
examined the problem of discrepancies between theoretical winds predicted by 
global models and those observed under magnetically disturbed conditons. 

These have considered altitudes above 150 km; however, results obtained here 
are pertinent to those discussions. Basically the issue comes down to whether 
ion drag (collisional coupling) or auroral heating (Joule and charged particle 
precipitation) is primarily responsible for deviations from the models. 
Stoffregen (1972) and Rothwell, et al. (1974) contend that ion drag cannot 
account for the observations and that pressure gradients due to heating are 
the primary driving forces. On the other hand, Rees (1971a) and Meriwether, 
et al. (1973) favor ion drag, particularly in the evening and midnight sector; 
and the latter specifically argue against pressure gradients as accounting for 
discrepancies at other times. 
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The observations of groups on both sides of the question are based on 
chemical releases, which are somewhat brief (lasting a few minutes) and are 
restricted in local time coverage. The arguments are semi-quantitative in 
nature, relying heavily on the relative directions of the ion and neutral 
velocities. The more recent papers on both sides cite the results of Fedder 
and Banks (1972) as indicating that the equilibrium directions for ion and 
neutral winds are coincident above about 200 km. In this regard the neglect 
of coriolis effects by Fedder and Banks (1972) is significant. Figure 6-8d 
shows that a directional difference between ion and neutral velocities is 
maintained in a steady state, even at higher altitudes, when the coriolis 
force is taken into account. From the detailed numerical results, at 250 km 
v n is rotated in a clockwise direction through an angle of A7 degrees from v^ 
in this calculation. This angle develops over a period of time, but even at 3 
hours into the calculation this angle has reached a value of about 30 degrees. 
The time for development and the final equilibrium angle depend on the model 
values used. At lower altitudes the angles can be somewhat larger (e.g. about 
56 degrees for 150 km at 24 hours). Such effects should be considered in 
future discussions of ion drag. 

An additional consideration is the variability of ion velocity. Figures 
6-3a, b show that the electric field (hence ion velocity) can change substan- 
tially in magnitude and direction while neutral winds at high altitudes are 
presumably less variable. Thus, observations over relatively short periods of 
time may be subject to considerable deviation from the general trend over a 
period of one or more hours. Comparison of ion and neutral velocity direc- 
tions for brief periods might then be quite misleading. 

Results of these calculations indicate that both ion drag and pressure 
gradients due to auroral heating are important. While heating appears to 
predominate during this experiment, the effects of ion drag are large enough 
to modify the winds significantly, even though not decisively. Moreover, it 
is not obvious that these results apply generally. Both Joule heating and ion 
drag result from the same physical phenomenon — ion collisions with neutrals. 
Joule heating varies with |v^ - | ^ while ion drag varies with (v A - v n > • 
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It therefore appears that as)^ - v^J increases, a transition will occur, 
below which ion drag will dominate the bulk motion and above which pressure 
gradients due to Joule heating will assume that role. An important factor in 
this consideration, however, is the spatial distribution of electric fields, 
which is the driving function for the ions. This precludes any simple esti- 
mate of such a threshold. Furthermore, if particle precipitation heating is 
larger than Joule heating the question is moot. The point is that under 
different conditions, also to be found in the high latitude ionosphere, ion 
drag may predominate. 

In summary, the following conclusions are drawn the calculations. 

Collisional coupling of neutral particles with ions driven by large electro- 
static fields can cause neutral winds of several hundred meters per second, 
having considerable altitude structure below about 130 km. Coriolis forces 
have the effect of maintaining ion-neutral velocity differences, even in a 
steady state for constant electric fields. Attainment of a steady state at 
low altitudes requires longer times when coriolis forces are taken into account 
than in their absence. 

When observed electric fields and electron number densities are used in a 
simple collisional coupling calculation, short time variations in these prop- 
erties are not reflected in the resulting neutral wind at a given altitude. 

The neutral gas inertia is so large relative to the rate of collisional momen- 
tum transfer, that such variations are averaged out. However, when weighted 
altitude averages of the neutral velocity are computed, short time scale 
perturbations of up to 50 percent are found to occur, due to fluctuations in 
the electron density weighting factor. 

Comparison of weighted altitude averages of neutral velocities computed 
by direct numerical integration with those obtained indirectly from similar 
averages of ion velocities by means of the incoherent scatter radar data 
analysis procedure reveals: (1) magnitudes are underestimated by up to 40 

percent; but (2) temporal structure and general directions are reasonably well 
reproduced. Magnitudes of velocities determined from the radar data are also 
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somewhat sensitive to the collision frequencies used, which includes the 
interaction model, the ion composition model, and the neutral atmosphere 
model. Some of the temporal structure is associated with variability in 
ionization density profiles due to the fact that weighted altitude averages of 
ion velocities are observed. With due recognition of these factors, however, 
it appears that neutral velocities derived from the incoherent scatter radar 
observations give a reasonable indication of the time evolution of neutral 
winds . 

From a direct comparison of neutral velocities derived from the incoherent 
scatter radar observations with the weighted altitude averages of the calcu- 
lated velocities, it is found that ion drag is inadequate to explain the 
observations. An equation is derived to compute the pressure gradient required 
to remove the discrepancies between observed and calculated results. The 
magnitude of the acceleration associated with this pressure gradient is found 
to vary from a few centimeters per second per second in quieter periods to a 
few tens of centimeters per second per second in more active periods. Accel- 
erations due to ion drag are of the same order. These are generally at least 
an order-of-magnitude larger than those associated with solar heating. Direc- 
tions of the pressure gradients are consistent with the location of the radar 
relative to the auroral oval and expected auroral heat sources. Fluctuations 
in the pressure gradients, to the extent that they are real, are thought to 
result from shifting balances in highly variable heat sources and from back 
pressures . 

It appears that ion drag and pressure gradients are both required to 
explain the neutral winds generated in the auroral E region during geomagnetic 
disturbances . 
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Section VII 
SUMMATION 


7.1 SUMMARY AND CONCLUSIONS 

This work is comprised of two basic parts: a formal development of the 

Conservation equations for ionospheric gases, and an application of these 

# 

equations to ionospheric dynamics, in which the generation of neutral winds 
'during geomagnetic disturbances is investigated. Governing equations for 
ionospheric gases are developed from kinetic theory, beginning with the Boltz- 
mann equation; and a particle viewpoint is maintained throughout. The final 
formulation is based on the three-fluid approximation, which requires that the 
electron, ion, and neutral gases be treated separately, with different temper- 
atures and flow velocities for each. The domain of applicability of the final 
equations is approximately 90 km to 800 km in altitude for all geographic 
locations. 

Starting from the Boltzmann equation a formal equation of transport for 
an arbitrary particle property is derived in Section II. This property is 
successively specified to be mass, momentum, and kinetic energy, resulting in 
the species conservation equations for these properties. Collision terms are 
treated separately in Section III in the binary collision formalism of the 
Boltzmann collision integral. Formal results are extended by approximating 
the species distribution functions as independent, displaced Maxwellian distri- 
butions. In this context a new formalism is developed for the collision terms 
in the momentum and energy equations; it simplifies calculations for specific 
interactions, as compared with similar expressions presently in the literature. 
Explicit evaluations of collision terms are carried out in the Appendix for a 
sufficient variety of interactions to meet most ionospheric requirements. In 
addition, proofs of some of the limiting properties of the general expressions 
are provided there. 

In Section IV, the formalism of the species conservation equations is 
applied to the ionosphere in terms of the three-fluid approximation. This 
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results in separate sets of conservation equations for electrons, ions, and 
neutral particles, coupled through the collision terms. Extreme values of 
ionospheric properties and parameters are considered for the altitude region 
90 km Ito 800 km, for which the three-fluid approximation is assumed to hold. 
Based oh "these values, order-of-magnitude estimates are made to determine the 
range of values possible for each term in the complete set of conservation 
equations. Those terms which can be clearly identified as negligible under 
all ionospheric conditions for this altitude region are omitted. The resulting 
set of equations provides a reasonable starting point for a wide variety of E 
and F region ionospheric investigations, from purely theoretical model studies 
to the data analysis of ground-based or in situ observations. 

In order to close the set of conservation equations , the higher order 
velocity moments introduced into the equations by the derivation method have 
to be related to lower order velocity moments, essentially independently of 
the set of conservation equations. This is accomplished in Section V by means 
of transport tensors and coefficients. Electron transport is treated in a 
comprehensive manner by the method of Shkarofsky (1961). In order to take 
advantage of simplifications in the use of this method when electon-neutral 
collision frequencies vary with relative velocity as a power law, a new method 
is derived for approximating an arbitrary collision frequency velocity depen- 
dence as a power law. The equations presented represent a condensation of 
Shkarofsky 1 s procedure, for application to ionospheric calculations. Ion and 
neutral transport coefficients available from the literature are also provided. 
Explicit sets of collision frequencies for electron-ion, electron-neutral, and 
ion-neutral collisions are presented, as computed from results in Section III 
and the Appendix, or as taken from the literature. Those sets of collision 
frequencies incorporating the effects of flow velocity differences have not 
previously been available. 

The second basic portion of this work applies the results of the first 
five Sections to a specific problem in Section VI. In the process, it illus- 
trates the different methods available for modifying the set of initial equa- 
tions (the products of the first part) so that numerical solutions can be 
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obtained. These methods include additional approximations based on the physical 
conditions of the particular problem, use .of empirical or analytical models as 
substitutions for equations or variables, use of observational data in the 
same manner, and mathematical assumptions required for consistency with these 
other techniques . 

The problem investigated involves neutral winds observed during magnetic- 
ally disturbed conditions in the auroral E region. One objective of the 
calculations is to determne the effect of altitude structure in the ion and 
neutral winds on the determination of neutral winds from incoherent scatter 
radar observations. A second objective is to determine the relative importance 
of collisional coupling (ion drag) and auroral heating in the generation of 
these winds. Observations made by the incoherent scatter radar at Chatanika, 
Alaska, on 15 May 1974 are used in the calculations and analysis. 

Theoretical calculations are based on large auroral electric fields 
causing ions to move, with collisions between ions and neutrals both modifying 
the ion flow and causing neutral gas flow. Viscous and coriolis terms modify 
the neutral gas motion; pressure terms are omitted from the initial phase of 
the. model. Calculations are restricted to altitudes between 90 km and 250 km. 

To become familiar with the effects of the coriolis term and inertial proper- 
ties of the neutral gas, model calculations are first carried out with a 
simple electric field model. From these calculations it is determined that 
the neutral atmosphere responds much more slowly to ion drag at lower altitudes 
(<125 km) than at higher altitudes (>150 km). The coriolis term causes a 
substantial ion-neutral velocity difference to be maintained in a steady- 
state, resulting from both a smaller neutral velocity magnitude (compared with 
the ion velocity) and a difference in direction of about 45 degrees. Due to 
the large inertia of the neutral gas relative to the ion gas at low altitudes, 
effects of initial conditions can last for many hours at these altitudes, when 
ion drag is the only driving force. 

For the primary calculations, observed electric fields and electron 
density profiles are employed, providing 6-minute time resolution. Despite 
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the considerable time variation in these quantities, the neutral gas effec- 
tively filters out such variability due -to its large inertia relative to the 
ion gas. Large vertical gradients in horizontal ion and neutral velocity 
components are found to result from the ion drag model. To make comparisons 
with observed velocities, weighted altitude averages of the calculated veloc- 
ities must be taken numerically. By this means it has been possible to assess 
semiquantitatively the technique for deriving neutral velocities from the 
incoherent scatter radar observations. It has been found that this data 
reduction technique systematically underestimates the velocity magnitude by up 
to 40 percent , but provides a fairly good indication of temporal variation in 
the weighted altitude average of the neutral velocity. A sensitivity to the 
collision frequencies used in the data reduction has also been noted. In 
addition, it has been determined that the process of taking weighted altitude 
averages itself can introduce fluctuations in magnitude of up to 50 percent in 
resultant velocities. This is due to variation in the electron density profile 
part of the weighting factor, together with the altitude structure in the 
velocity fields. 

Comparison of neutral velocities derived from observation with those 
calculated for E region altitudes shows that the ion drag model used for the 
calculation is inadequate by itself to account for the observed velocities. 

An equation is derived to estimate the pressure gradients required to resolve 
the discrepancy between the calculated and observed velocities. The acceler- 
ations resulting from these pressure gradients are generally comparable to 
those associated with ion drag, but at least an order of magnitude larger than 
those due to solar heating. Directions of the pressure gradients are consis- 
tent with expected locations of auroral heating, relative to the radar location. 
Since accelerations due to pressure gradients and ion drag are of the same 
order, it is concluded that both are important in generating and modifying 
neutral winds in the auroral E region during geomagnetic disturbances. 

7.2 RECOMMENDATIONS FOR FUTURE RESEARCH 

In the course of preparing a work such as this, certain weak areas in 
foundations on which sometimes elaborate calculations are based become apparent. 
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One of these is ion-neutral collision frequencies. The lack of knowledge of 
ion-neutral interactions at intermediate and close ranges leaves large uncer- 
tainties in ion-neutral collision frequencies at temperatures above about 
300°K. One consequence of this in the present work was seen to be a signifi- 
cant uncertainty in the absolute magnitude of the neutral velocities derived 
from observations. 

Another weak area, loosely related to the first, is transport properties 
of ions in a partially ionized gas permeated by a magnetic field of arbitrary 
strength. A complete theoretical treatment, similar to that presented for 
electrons, is needed. Because the ion gas undergoes a transition from large 
ion-neutral collision frequencies to &mall collision frequencies relative to 
the ion gyrof requency , it is clear that transport transverse to the magnetic 
field is more complex and more significant for ions than for electrons (for 
which it is largely inhibited). Large ion gas inertia, as compared with the 
electron gas, and closer collisional coupling to the neutrals result in ions 
being less mobile than electrons. However, charge neutrality assures that 
macroscopic densities of ions and electrons are the same. Thus, ionization 
transport due to mechanical forces are constrained by the ion motion, although 
electrons play the primary role in responding to electromagnetic forces. At 
present it is difficult to separate the limitations imposed by lack of a 
comprehensive theoretical transport treatment from those due to uncertainties 
in collision frequencies; both require improvement. 

With respect to the development of conservation equations for the E and F 
regions, several modifications could be made to extend the range of applica- 
tions. These are, for the most part, straightforward, and require only the 
addition of certain equations. First, problems of ion or neutral composition 
can be treated by adding the appropriate continuity equations for each ion or 
neutral species. Second, the set of conservation equations can be augmented 
by Maxwell's equations in a form appropriate to whatever problem is under study. 
This would be necessary, for example, in studies involving hydromagnetic waves 
or self-consistent electric fields resulting from dynamo action. Third, for 
higher altitudes (>800 km), where light ions are dominant, the three-fiuid 
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approximation may not be reliable. In this case separate sets of conservation 
equations would be required for each species. 

In the study of auroral E region neutral winds, several extensions and 
improvements are suggested. First, data for a less disturbed case could be 
used for a case study similar to this one. This could be done rather easily 
by changing the data input to the computer programs. The primary purpose of 
such a calculation would be to determine if ion drag is more important than 
auroral heating under these conditions. Published observations under such 
conditions indicate that this may be the case. A second addition would be the 
calculation of electric currents and resulting magnetic field perturbations, 
and calculation of Joule heating. Both could be accomplished by extending the 
calculations already made. 

Significant improvements in the present neutral wind calculations can be 
made, but at a significant increase in computational effort. It was demon- 
strated that pressure gradients due to auroral heating play a dominant role in 
the generation and modification of neutral winds in the auroral E region 
during magnetic disturbances. In order to include these pressure gradients in 
a calculation, horizontal variations would have to be taken into account, 
requiring at least two spatial dimensions, ideally three. In addition, the 
neutral continuity and energy equations would be required to replace the 
static neutral model atmosphere. This would have the added desirable effects 
of allowing the calculation of vertical velocities (with consequent effects on 
horizontal velocities) and of permitting inclusion of neutral composition 
changes, which have been observed during magnetic storms. To accomplish this 
in the manner of the present study would require observational data from a 
considerably expanded spatial region. However, initially some of the effects 
could be profitably explored in more restricted model calculations. Much 
remains to be done. 
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Appendix 


ANALYTIC EVALUATION OF COLLISION TERMS 

A.1 LIMITING FORMS FOR THE COLLISION FREQUENCY AND ITS DERIVATIVE AT 
COMMON FLOW VELOCITIES 

In order for the generalized collision frequency, derived in Section 
III, to be consistent with previous work, it should reduce to the standard 
expression in the limit of common flow velocities. In this part of the 
Appendix this is demonstrated to be the case. It is also shown that in 
this limit the derivative of the collision frequency with respect to the 
relative flow velocity vanishes. 


From equation (3-35) the effective collision frequency for momentum 
transfer is given by 
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where 
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Clearly, simply setting v q to zero in equation (A-l) is meaningless. 
Instead, a series expansion is used for the quantity in brackets: 
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Making use of the chain rule and equation 
to be rewritten 
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In the limit v q -*■ 0, only the first term of the series contributes, 
giving 
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This is the expression derived by Banks (1966a) for two static gases at 
different temperatures. It is also equivalent to the expression desig- 
nated by Itikawa (1971) as < ^ e ff > when applied to electrons. 


To evaluate 9v rs /9v Q , l et t * ie integral in equation (A-l) be denoted 


by J : 
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Carrying out the differentiation of v , using equation (A-l), gives 
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Since the limit of the sum equals the sum of the limits, each term could 
be examined independently. However, this turns out to be useful only 
for the first term; the others must be treated collectively because 
important cancellations take place which must be taken into account 
before the limiting process. 


For notational convenience let the successive terms in equation 
(A-6) be designated T^ , T^ » T^, and T^. From equations (A-l) and (A-5) 
the first term may be written 


T = - 2K v v (v ). 
1 o rs o' 
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Taking the limit as v q vanishes gives 

lim T- = 0 

v -K) 
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since, by equation (A-4) , is finite in this limit. From equation 
(A-5) the fourth term can be written 
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The second and fourth terms can be combined immediately resulting -in 
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Making use of the series expansion (A-2) in J gives 
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The hyperbolic cosine in the third term of (A-6) is similarly 
expanded with the result 
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Comparison of equations (A-ll) and (A-12) shows that the series are both 
multiplied by the same factors so that they can be combined directly 
term by term. The result is 
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where the first terms cancelled completely. From this equation it is 
clear that 


lim (T, + T- + T 4 ) * 0 . (A-14) 
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Combining equations (A-6) , (A-8) , and (A-14) gives the final result 


lim ° = lim (T, + T. + <T„ + T, ) = 0 . (A-15) 
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A.2 ANALYTIC EVALUATION OF COLLISION FREQUENCIES FOR SPECIFIC FORMS 
OF INTERACTION 

Forms of interaction among charged and neutral particles of the 
upper atmosphere are varied, and, in some cases, poorly known. In this 
section, collision frequencies and their relative velocity derivatives 
are determined analytically from equation (3-35) for those forms likely 
to be required for most ionospheric applications. Some of these are 
well known theoretically (e.g. , Coulomb and hard sphere interactions), 
while others are useful for the curve fitting of experimental data (e.g. 
polynomials or power laws in relative velocity). Specifically, the 
interactions considered have the forms: 

Q rs (g) * Aj g j (A-16a) 

Q rs (g) “ \ g V (A-16b) 

Q rs (g) = Ag exp(-Bg) , (A- 16c) 

where j is an integer, v is nonintegral, and 6 is arbitrary. Since 
closed- form expressions can be obtained for integral powers, but not 
for arbitrary nonintegral powers, the distinction between j and v is 
useful. 


The more familiar interaction laws are usually written as functions 
of relative distance. When the force between two particles varies 
inversely as a power of the distance between the particles. 
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Chapman and Cowling (1970, Chapter 10) demonstrate from classical two- 
body interaction theory that the momentum transfer cross section (not 
their terminology) varies with relative velocity as 


Q rs (8) 


« 8 


-4/(y-l) 


With this expression, the velocity variation of some familiar power law 
interactions can be identified. Thus, for the inverse square Coulomb 
force, y * 2 and j * -4. For ’’Maxwellian molecules", so called because 
Maxwell first recognized the simplifications resulting from this type 
of interaction (Chapman and Cowling, 1970, p 173), the force law y - 5 
corresponds to J a -1. For the case of hard sphere interactions, y 00 
and j - 0 results, so that «.he momentum transfer cross section is inde- 
pendent of velocity. 


For present purposes, the constants of proportionality in equations 
(A-16) are not of interest; they may be obtained from theory (analytically 
or numerically) or from experiment. The object here is to evaluate the 
appropriate integrals. This is done for -4 £ j _< 4 and for arbitrary v 
and 0 within the mathematical constraints of the integrals. 


A.2.1 Integral Power Law Interactions 

From equation (3-35) the generalized momentum transfer collision 
frequency for species r with species s is given by 
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where v^ = l v r “ v s l is t * ie relative flow velocity. Let J denote the 
integral in this equation. With the momentum transfer cross section 
represented by equation (A-16a) , J may be written 
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Integrating by parts gives 
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The first term vanishes for j > -4. However, the integral is qualita- 
tively different for j < -2 than for j _> -2, so the cases j = -3, j - -4 
will be treated separately. For j < -4, J_. will not be finite. 

First the case of j _> -2 is treated. With the notation a 
and the change of variable z = J* g, equation (A-19) becomes 
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Since both terms in bracket's are even functions of z the integrand has 

definite parity. The integrals to be evaluated all have the form 
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Similarly, if 8 is an even integer, the corresponding result is 
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where 


sinh(2az) exp(-z ) dz 


(A-24) 


Thus, once I and L q are evaluated analytically, all results desired for 

different values of j (8) can be obtained by differentiation with respect 

to a. It remains only to evaluate L and I . 

J o o 
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Consider L q first. If the hyperbolic cosine is expanded in terms of 
exponential functions and the square is completed in the resulting com- 
bined exponents, equation (A-22) can be rewritten as 


L Q “ 2 ex P( a2 ) 


f exp[-(z - a) 2 ]dz + f exp[-(z + a) 2 ] dz 
'o ' o 


(A-25)* 

Changing variables to x * z - a in the first integral and x - z + a in 
the second integral gives 


1 ( 2 jin 

2 exp (a )|j 


2 f°° 2 

exp(-x )dx + exp (-x )dx 

'a 


]• 


Since the integrand exp (-x ) is, an even function of x 


o 

L 


2 f a 2 

exp(-x )dx = exp(-x ) dx 
^ o 


(A-26) 


Equation (A-25) may then be written 


POO 

2 | 2 

L * exp (a ) exp(-x )dx 
* o 


= y ^ exp(a ) 

from standard integral tables. 


(A-27) 


Following a similar procedure for I q leads to the equation, analogous 


to (A-25), 


I Q ■ y exp(o 2 ) 


PCD pCD 

exp (-x )dx - exp (-x )dx 

J —fit * n 


(A-28) 


The symmetry of the integrand about zero, noted above, allows this equa- 
tion to be rewritten 


2 1 2 

I * exp (a ) exp (-x )dx 
° * o 


(A-29) 


This integral is recognized as the error function, which is conventionally 
defined as (Gautschi, 1964) 



(A-30) 


erf(z) = — f exp(-t 2 )dt 

JTJo 

With this definition, equation CA-27) is written 

I 0 -f exp(a 2 ) erf (a) . (A-31) 


Equation (A-27) can now be used to evaluate equation (A-21) for 
appropriate (odd) values of g. From equation (A-20) it is seen that if 
j is in the range -2 £ j _< 4, g will lie in the range 0 £ B £ 8. Carrying 
out the required differentiation gives the following results: 


g * 1: 


z sinh (2az) exp(~z ;dz 


= ^ a exp (a 2 ) 


(A-32a) 


g-3: 


3 2 

z sinh (2az) exp(-z )dz 


- a (3 + 2a 2 ) exp (a 2 ) 


(A- 3 2b) 


5 2 

3 = 5:1 z sinh (2az) exp(-z )dz 


■r 

J o 


a (15 + 20a 2 + 4a 4 ) exp (a 2 ) 


(A-32c) 


3 = 7: 


7 2 

z sinh (2«z) exp(-z )dz 


-iff 

16 


a (105 + 210a 2 + 84 a 4 + 8a 6 ) exp (a 2 ) 


(A-32d) 


For even g s equation (A- 31) must be used in conjunction with 
equation (A-23) . This requires the derivative of the error function, 
which is given by (Gautschi, 1964) 


erf (z) = — exp(-z 2 ) 

/T 


(A-33) 
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Differentiation then gives the following results: 


»r 

J a 


3*0: slnh (2az) «xp(-s )dr 

° 

- exp(o 2 ) erf (a) 


(A- 34a) 


3-2 


* A 


slnh (2az) exp(-z )d* 


• f + exp (a 2 ) erf (a) (1 + 2o 2 ) 


(A- 34b) 


-f 

* A 


i 4 2 

3 * 4:1 z slnh (2az) exp(-z )dz 


- | (5 + 2a 2 ) + exp(a 2 ) erf (a) (3 + 12a 2 + 4a 4 ) (A-34c) 

3-6 :J~ z 6 sinh (2az) exp(-z 2 )dz - | (33 + 28a 2 + 4a 4 ) 


+ exp(a 2 ) erf (a) (15 + 90a 2 + 60a 4 + 8a**) 


(A-34d) 


3 - 8: J* z 8 sinh(2az) exp(-z 2 )dz - (279 + 370a 2 + 108a 4 + 8a 6 ) 

+ § exp (a 2 ) erf (a) (105 + 840a 2 + 840a 4 + 224a 6 + 16a 8 ) 


(A-34e) 


These results, equations (A-32) and (A-34), aay now be used as 
required to evaluate equation (A-20) for different values of j. Beginning 
with j - -2, the following results are obtained straightforwardly. 


J - -2: J 


^ J** (z 2 - 1) sinh (2az) exp(-z 2 )dz 


- [2a +y7exp(a 2 ) erf (a) (2a 2 - 1)] (A-35a) 
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j ■ -1: (Maxwellian Moleculea) 

A 

-i 

1/2 


J -1 


— f*(2a 3 - 3a) aiah(2aa) exp(-a 2 )da 
2a% 1 J o 


A . ( 


2K 


3/2 


exp(a j 


(A 


j - 0: (Hard Spheres) 


A 

3 o m ~ S 2 
° oK 2 >o 


J (a* - 2a 2 ) sinh (2a a) exp(-a 2 ) da 

* O 


— ^ (2a(l ♦ 2a 2 ) 
8aK 


+ ^jTexp(a 2 ) erf (a) (-1 ♦ 4a 2 ♦ 4aSl 


(A 


j “ 1 


: Jj - j (2a 5 - 5a 3 ) sinh (2aa) exp (-a 2 ) da 


•\y 2 ' /r exp(a 2 ) (5 + 2a 2 ) 


(A 


A «co 

j « 2: J_ - — | J (a® - 3a 4 ) sinh (2aa) exp(-a 2 ) da 

aK J 


■— [2a(3 + 16a 2 + 4a 4 ) 


16 aK 


+ y7exp(a 2 ) erf (a) (-3 + 18a 2 ♦ 36a 4 ♦ 8a*) J (A 


i " 3: J 3 “ 7/2 

J 2aK ,,Z 7 o 
n 

A- a 


E (2z 2 - 7a 5 ) t|Lnh (2aa) exp (-a 2 ) da 


" ~^T72 exp(a 2 ) (35 + 28a 2 + 4a 4 ) 


8K 


(A 


-35b) 


-35c) 


-35d) 


35e) 


35f) 
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j ■ 4: J, ■ -~r [ (z 8 - 4z 6 ) Binh (2az) exp(a 2 )dz 

4 aK 4 J o 

■.* ; ■' - ■- ■'■■■ • ■' -t . 

- [2a(15 + 146o 2 + 76a 4 + 8a 6 ) 

32aK 4 

-} ■ 

+ JIT exp (a 2 ) erf (a) (-15 + 120a 2 + '360a 4 
+ 160a 6 + 16a 8 )] . , (A-35g) 


Clearly this procedure could be used to extend results to higher 

{ * ' * '• ‘ * 1 ’’ , 

values of j* However, for the present, these should be sufficient for 
most curve-fit applications. Lower values of j appear to be limited 
by properties of the integrand In equation <A-18) to j - -3 and -4; these 
must be treated individually. For j - -3 equation (A-18) becomes, with 
the change of variables and notation introduced above. 


-3 


— r 

fK>< 


z exp(-z ) d 


[ sinh (2az) I 

. j 


-=L .i& (2..) 2) 

JT 1 2ot 



slnh (2az) 
2az 


(1 - 2z 2 ) exp(-z 2 )dz 


9 


(A- 36) 


where the second equality follows from an integration by parts. 


The first term in equation (A- 36) vanishes in the limits, leaving 
two integrals. One may be evaluated from equation (A- 32a) ; the other 
is 


r 4 -3 p sinh (2az) 

2a /IT Jo 2 


exp(-z z )dz 


-3 


2a /T 


I_ 3 (a) 


(A-37) 


which defines I_^(a). Differentiating I_ 3 with respect to the parameter 
a gives 
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r 

* n 


cosh (2az) exp (-z )dz 


= 2 L q • fT exp(a^) 


» 


(A- 38) 


where equations (A-22) and (A-27) have been used. Since I_^(0) - 0, 
the differential equation (A- 38) can be immediately integrated to obtain 


I_ 3 (o) 



= -5T erf(ia) 


(A- 39) 


(Gautschi, 1964). The final result is then 


-3 


-3 

2 y/iT 


jes 


■ 2 > 


exp(o ) + ^ erf (ia) 


3 


(A-40) 


For j ■ -4 a similar procedure is followed. In the new notation 
equation (A-18) becomes 


J -4 * *-4 [ 

- ‘-4 [ 

<00 

+ ^ f sinh (2az) exp(-z^) dz 


(A-41) 


where an integration by parts has been performed. The first term 
vanishes at the upper limit, but has a limiting value of 1 at z • 0. 
The integral in the second term is just I q , equation (A-24) . So the 
desired result can be written immediately 






2 ^ ex P (<* ) erf (a) 


] 


from equation (A-31). 


(A-42) 
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With equations (A- 35) , (A-40) and (A-42) , collision frequencies 
can be evaluated from equation (A-17) • For notational convenience, the 
j value is given as a superscript on the collision frequency. Final 
expressions for the collision frequencies are, in the original notation: 

j - "4:v;*(v o ) - * - 3 -“ 4 - £erf (/K v q ) - 2v q J^ expC-Kv, 2 ) J (A-43a) 


j « -3:v^(v o ) 


j * -2:\> 2 (v ) 
J rs o' 


„ 8 '3" k + ± /i eXp( “ KV o 2) v o >] 

AA[" 2v Ji 

2K v 3 L <W * 


o' 

(A-43b) 


exp(-Kv Q ) 


+ (2Kv 2 - 1) erf ( Jk v ) 
o o 


] 


j = -lsv 3 (v ) = n A , 
J rs o s -1 


(A-43c) 

(A-43d) 


4K v 


[■*% 


(1 + 2K v 2 ) exp (-Kv 2 ) 
o r . o 


+ (-1 + 4K V 2 + 4K 2 V 4 ) erf ( 
o o 




1 n s A 1 2 

j = lsv _r(v ) = (5 + 2K v ) 

rs o zK o 


(A-43e) 


(A-43f) 


, „ 2, „ s A 2 

j = 2:v (v ) - — r — r 
rs o 8 r 3 v 3 




+ 16Kv 2 + 4K 2 v 4 )exp(-Kv 2 ) 
o o r o 


+ erf( /Fv )(-3 + 18Kv 2 + 36K 2 v 4 + 8K 3 v 6 ) (A-43g) 

f o o o o I 


n A 

j = 3;v 3 (v ) - ~ 8 3 (35 + 28Kv 2 + 4K 2 v 4 ) 

rs o o o 


(A-43h) 
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4 n s A 4 

i - ^ > -Ti 

8 ° 16K v J 


fyr 


v (15 + 146Kv 2 + 76 K 2 v 4 
o o o 


+ 8K 3 v 6 ) exp(-Kv 2 ) + erf ( 7 k*v ) (-15 + 120Kv 

O O w O C 


+ 360K 2 v 4 + 160K 3 v 6 + 16K 4 v 8 ) 
o o o 


] 


(A-43i) 


For the momentum equation (3-52) , collision frequencies as given 
by these equations represent the entire collisional interaction. How- 
ever, in the energy equation the derivative of the collision frequency 
with respect to the relative flow velocity is also required. These 
expressions are given in the next equations. 


j 

j 

3 

j 

i 

j 


-4: 


-4 

3v n A , 
rs _ s -4 

3v 4 


[”./F 


" exp(-Kv 2 ) (3 + 2Kv 2 ) 
^ o o 


- 3 erf ( 
-3 


fv] 

. . iAa r 

2v 4 L 

° -I 

. erf(i yiTv o )(3 + 2Kv q 2 ) I 

h/r 


6v 0 + exp(_Kv o fc) 


„ -2 

3v n A n 

rs _ s -2 

9v 0 4 

o 2K v 

o 


— exp (-Kv ) 
o 


+ erf ( yF v ) (3 - 2Kv 
-1 


■>] 


3v 

-1:—^ = 0 


8v 

< 

av 


0: 


n A 
rs s 


9v sir 2 

o 4K 


^ [* JT-. 


exp (-Kv 2 ) (2Kv 2 - 3) 
o r o' o 


+ erf ( TFv ) (3 - 4Kv 2 + 4K 2 V 4 ) 

¥ O O ° I 


av 


1: 


rs 


= 2n A v 
3v s 1 o 


(A-44a) 


(A- 4 4b) 


(A-44c) 

(A-44d) 


(A-44e) 
(A-44f ) 
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j = 2: 


3v 3n A 0 

rs s 2 


3v qtt-3 4 

o 8K v 

o 


[•/? 


v exp(-Kv 2 ) (-3 + 4Kv 2 + 4K 2 v 4 ) 
o o o o 


+ erf ( JF v ) (3 - 6Kv 2 + 12K 2 v 4 + 8Kv 6 ) 
■ o o o o 

3 


] 


3v ~ 2 n A 0 v 0 

j - 3 :t— ~ = 3 3 ° (7 + 2Kv 2 ) 

J 3v K o 


(A-44g) 

(A-44h) 


j - 4: 


3v 5n A. 
rs s 4 


3v A 4 

o 16K v 

o 


[•/? 


v exp(-K v 2 ) (-9 + 18Kv 2 
O’ o o 


+ 44K 2 v 4 + 8K 3 v b ) + erf ( /Fv )(9 - 24Kv Z + 72K Z v 

O O W O O ( 


r3 .. 6 


2.. 4 

o 


+ 96K 3 v 6 + 16K 4 -v 
o 


.*> 


(A-441) 


With these relations and the collision frequencies in equations (A-43) , 
the collisional transfer of energy and momentum can be evaluated for 
arbitrary temperatures and flow velocities for a variety of interactions. 


A.2.2 Nonintegral Power Law Interactions 

Experimental data of momentum transfer cross sections or results of 
numerical calculations may be fit to power law curves of arbitrary power. 
The appropriate integrals for this case are evaluated here. Because of 
similarity with the previous case, the analysis is the same through 
equation (A-20) . So here, the starting point is 


J = ^ (K)- (v+4)/2 
v Za j 


sinh (2az) exp (-z ) 


b V+2 [2z 2 - (v + 4) ]dz 


(A-45) 


where v is nonintegral. As previously, this equation holds only for 
v > -4. 


Several approaches can be taken for treating this integral; however, 
all appear to end up in infinite series or in one or more of the less 
common special functions. For this reason, the most expeditious route 
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to a series solution is taken, by expanding the hyperbolic sine in a 
power series: 


sinh (2az) = 


(2az) 


2n+l 


n-0 < 2 ” + » 1 


Equation (A-45) can then be rewritten 


(A- 4 6) 


J 

v 


K ; f2a) 2n+1 L r 

2»K (v+ « /2 „-0 < 2n+1 >' Jc 


exp (-z 2 ) z ( 2n+v+ ^dz 


- (v + 4) | exp(-z^) z^ n+v+ ^^dz | , (A-47) 

where now \) > -3 is required for all integrals to be finite. From one 
of the integral representations of the gamma function, the following 
result is obtained: 



(A-48) 


which has the form of the integrals in equation (A-47). Use of this 
equation in (A-47) gives 


J 

v 


A v r (2q) 2n / 2n + v + 4 

k (v+4)/2 ^ (2n + 1)! n r V 2 


where the gamma function property 


(A-49) 


T(z + 1) = z r(z) 
has been used. 


(A- 50) 


For the simple case v = -1, it can be verified that equation (A-49) 
is equivalent to equation (A- 35b). However, the path is sufficiently 
devious and tedious that the approach used previously for treating 
integral power laws appears shorter and less susceptible to error than 
working back through the power series from the general result given in 
equation (A-49) for arbitrary v (restricted only by v > -3). 
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The collision frequency is found from equation (A-17) and (A-49) to 
be given by 


v 


V 

rs 


<V 


2n A exp(-Kv 2 ) « (4Kv 2 ) n 

s v _ o y o 

f kW* » o 2 Jl < 2 " + »' ° 


/ 2n + v + 
\ 2 



(A-51) 


Differentiation of equation (A-51) with respect to v q has the result, 
after some manipulation of series. 


3v V 2n A exp (-Kv 2 ) ~ (4Kv 2 ) n (n - 1) 

rs s v r o v o _ 

9v o ' n-2 <2 » + 1)! 

* o 

• [n(v + 3) + 1] r ( 2n + - 2 - v + 2 ) . (A-52) 


With equations (A-51) and (A-52) the collision terms in the momentum 
and energy equations can be evaluated for the nonintegral power law 
interaction. 


A. 2.3 Exponential Law Interaction 

The final case to be treated is that in which the momentum transfer 
cross section varies with velocity according to an exponential law, 
equation (A-16c) . Then the integral in equation (A-17) becomes 

r 2 4 T 8lnh <2K8 Vl 

J 6 - Ag exp C-Kg 2 - Sgjg 4 d |_ 2kgv o J • < A ' 53) 


Expanding the differential and changing variables permits equation (A-53) 
to be rewritten 


2oK 


r 2 r 3 

exp(-z - yz) I z cosh (2az) 
)Jdz , 


- ^ sinh (2az) 


(A- 5 4) 


where z - Kg 2 and y = 3/ /k 7 The* hyperbolic sine and cosine functions 
can be rewritten in terms of exponentials to obtain 
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J 3 


4aK 


exp[-z - (y - 2a) z] 


? / o 2 

exp[-z Z - (y + 2a) z] (z J + -| 
o \ a 


y-iy 

(z 3 + fe) dz J . 


(A-55) 


These integrals all have the form 


fOO 

* t 2 X n , 

[ = exp (-z - pz)z dz 

n 

' n 


(A-56) 


As previously, the integrals for arbitrary n > 0 can be obtained 
from the integral corresponding to n = 0 by differentiating with respect 
to the parameter p: 


n * 

* a 1 

I = (-) n — - 

n 9p n 


where 




exp(-z - pz)dz 


(A-57) 


(A-58) 


Equation (A-58) can be treated by completing the square in the exponential 
and changing variables; the result is 


h = exp (i p2 )[ x - erf (!)] 

Differentiation with respect to p gives 
* 


(A-59) 


(A-60a) 


if = £ { _2 + F 9 exp (i p2 ) [ x ■ erf (f)] } 

- 2p + {7 exp £l - erf + P 2 J | (A-60b) 


2 * 

3 1 , 

o 1 


3P 

.3 * 


^i 2 = i6 - 2(4 + p2) + \F p exp - erf (f)][ 6 + p2 ]| ' 


(A-60c) 
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These equations are used together with equations (A-56) and (A-57) 
to evaluate (A-55) with p - Y - 2a in the first integral ,and p - a + 2« 
fin the second integral, giving the result 


J 


e 


Sj <4a[2 + y 2 + 4a 2 ] 

32alT ( 

- J7exp [jX^_2«) J [i - erf 
. [2 oy - 8a 2 + Y 2 + a(Y - 2a) 3 + 2] 


- f 1 exp J l 1 - erf (p-^)] 

• [2oy + 8a 2 - Y 2 + a(Y + 2a) 3 - 2] I . 


(A-61) 


Then in the original notation, the collision frequency, from equation 
(A-17) , is 



Since the limit of 3 * 0 corresponds to a momentum transfer cross section 

which is constant (hard sphere) , equation (A-62) should reduce to 
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equation (A~43e) in this limit. A quick check shows that it does, thus 
verifying all terms except those containing B* Finally the derivative 


of the collision frequency with respect to v q is given by 


3v j? n A exp(-Kv 2 ) 
rs _ s p o 

3v 


16K 2 v 4 
o 


3 - 2Kv 


«>[( 

[(f 


[‘ 

a--< 


4v 2 (k + e 2 ) 

o 


- 3 (£ +2 )] 


3 - 2Kv 

I 

2 /F 


■)] 


(2Kv - 3) 3 - 24K 2 v 4 + 323 Kv 3 
o o o 


+ 8v 2 (K - 3 2 ) - 3v 
o o 


- exp 


( 6 + r)- 3 ( 2 + r)] 

+ 2Kv o\ ~l 

2/T )j 


[(^) I 

Ih'-*) 

(£-*) + , (r +2 )]| 


(2Kv + 3) 3 + 24 K 2 v 4 + 163 Kv 3 
o o o 


+ 4v 2 (3 2 - 2K) + 3v 
o o 


(A-63) 


Comparisons of this equation in the limit 3 -* 0 with equation (A-44e) 
verifies the terms not containing 3- 


This concludes the analytic calculation of collision frequencies 
for nonequilibrium gases moving relative to one another. Results pre- 
sented here should be adequate for most ionospheric applications. 
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